VECMAtk: a scalable verification, validation and uncertainty quantification toolkit for scientific simulations.

We present the VECMA toolkit (VECMAtk), a flexible software environment for single and multiscale simulations that introduces directly applicable and reusable procedures for verification, validation (V&V), sensitivity analysis (SA) and uncertainty quantication (UQ). It enables users to verify key aspects of their applications, systematically compare and validate the simulation outputs against observational or benchmark data, and run simulations conveniently on any platform from the desktop to current multi-petascale computers. In this sequel to our paper on VECMAtk which we presented last year [1] we focus on a range of functional and performance improvements that we have introduced, cover newly introduced components, and applications examples from seven different domains such as conflict modelling and environmental sciences. We also present several implemented patterns for UQ/SA and V&V, and guide the reader through one example concerning COVID-19 modelling in detail. This article is part of the theme issue 'Reliability and reproducibility in computational science: implementing verification, validation and uncertainty quantification in silico'.

Hereditary angioedema with F12 mutation: factors modifying the clinical phenotype.

BACKGROUND: Hereditary angioedema (HAE) with normal C1 inhibitor (C1Inh) associated with the c.983C>A and c.983C>G mutations of the F12 gene (FXII-HAE) is a rare condition, and presents with highly variable clinical expression. On the basis of data gathered from a large carrier cohort, we assessed the modifiers affecting the clinical phenotype. METHODS: We analyzed clinical and biological data recorded from 118 mutation carriers (80 symptomatic and 38 asymptomatic), 58 noncarrier relatives from 40 families, and 200 healthy donors. Disease severity was scored in relation to frequency and location of edema, as well as age at disease onset. To predict FXII-HAE disease severity, we analyzed the biological phenotype [C1Inh, C4, spontaneous amidase, angiotensin-I-converting enzyme (ACE), aminopeptidase P (APP), and carboxypeptidase N/M (CPN)] by means of logistic regression (Akaike information criterion) and odds ratio (OR). RESULTS: Meaningful variables contributed to FXII-HAE, with the kinin catabolism enzymes ACE and CPN exhibiting a significant inverse relationship with disease severity (OR = 0.36, 95% CI 0.23-0.59, P < 0.001; OR = 0.58, 95% CI 0.36-0.91, P < 0.05, respectively). CPN activities were 37.5 (28.5-41.3) nmol/ml/min and 38.5 (32.8-45.6) for FXII-HAE asymptomatic and symptomatic carriers, respectively, and 37.9 (30.5-43.7) nmol/ml/min for noncarriers. Angiotensin-I-converting enzyme activities were 58 (44-76) and 49 (35-59) nmol/ml/min for FXII-HAE asymptomatic and symptomatic carriers, respectively, and 56 (49-66) nmol/ml/min for noncarriers. CONCLUSIONS: The FXII-HAE is associated with modifiers, for example kinin catabolism enzymes, ACE and CPN, different from those recognized in HAE with C1Inh deficiency.

Congenital myopathy with focal loss of cross-striations revisited.

In 1977 Wijngaarden et al. reported a Dutch family with a congenital myopathy characterized by external ophthalmoplegia and a remarkable histological feature, focal loss of cross-striations. A small number of other families with similar clinical and pathological features led to the consideration of this congenital myopathy as a distinct entity. Here we present more than 30years of follow-up from the Dutch family and report recently identified compound heterozygous mutations in the skeletal muscle ryanodine receptor (RYR1) gene, c.10627-2A>G and p.Arg3539His (c.10616G>A). Focal loss of cross-striations on muscle biopsy is another histopathological feature that should raise the possibility of RYR1 involvement.

Can bubble columns be an alternative to fibrous filters for nanoparticles collection?

The most effective and widely used dedusting techniques to separate nanoparticles of a carrier fluid are fibrous media. The main problem is the clogging of the filter that induces a pressure drop increase over time and thus requires a regular cleaning of the media (or its replacement). Following these observations, this study proposes to investigate the potential of bubble columns for nanoparticles collection. Despite collection efficiencies lower than those of fibrous filters, experimental results show that bubble columns present likely performances for the collection of nanoparticles and have collection efficiency even more important when the liquid height is high and bubbling orifices have low diameters. Experiments have also revealed the presence of a most penetrating particle size for a particle diameter range between 10 and 30 nm. The model developed in this article highlights a good agreement between the theoretical collection efficiency by Brownian diffusion and experimental collection efficiencies for particles lower than 20 nm. Nevertheless, the modelling may be extended to other collection mechanisms in order to explain the collection efficiency increase for particles higher than 20 nm and to confirm or infirm that electrostatic effects can be the cause of this efficiency increase.

Mutational spectrum and phenotypes in Danish families with hereditary angioedema because of C1 inhibitor deficiency.

BACKGROUND: Hereditary angioedema (HAE), type I and II, is an autosomal dominant disease with deficiency of functional C1 inhibitor protein causing episodic swellings of skin, mucosa and viscera. HAE is a genetically heterogeneous disease with more than 200 different mutations in the SERPING1 gene. A genotype-phenotype relationship does not seem to exist in HAE, although the polymorphism c.-21T>C of exon 2 has been reported to be associated with a more severe phenotype. We aimed to establish the mutational spectrum of C1 inhibitor deficiency in Denmark and investigate the possible disease-aggravating effect of the c.-21T>C polymorphism. METHODS: Hereditary angioedema was diagnosed based on clinical features and C1 inhibitor deficiency. A general severity score ranging from 0 to 10 was developed based on age at disease onset, clinical manifestations and treatment experiences. SERPING1 gene investigation was performed by exon sequencing followed by multiplex ligation-dependent probe amplification genomic rearrangement analysis in all known Danish HAE families. RESULTS: Fifty-nine patients with HAE from 26 families were included in this study. The mean disease severity score was 7.12 [1-10], and the mean C1 inhibitor function was 26% [20-46%]. The sensitivity of the mutational screening was 96%, and 13 new mutations were found in this Danish patient cohort. Nine patients (15%) carried the c.-21T>C polymorphism, but they didn't have a more severe phenotype. CONCLUSION: Thirteen new mutations were identified in the Danish HAE population. No correlation between the c.-21T>C polymorphism, the biochemical values of C1 inhibitor function and the clinical severity score was found.

Recessive RYR1 mutations cause unusual congenital myopathy with prominent nuclear internalization and large areas of myofibrillar disorganization.

AIMS: To report the clinical, pathological and genetic findings in a group of patients with a previously not described phenotype of congenital myopathy due to recessive mutations in the gene encoding the type 1 muscle ryanodine receptor channel (RYR1). METHODS: Seven unrelated patients shared a predominant axial and proximal weakness of varying severity, with onset during the neonatal period, associated with bilateral ptosis and ophthalmoparesis, and unusual muscle biopsy features at light and electron microscopic levels. RESULTS: Muscle biopsy histochemistry revealed a peculiar morphological pattern characterized by numerous internalized myonuclei in up to 51% of fibres and large areas of myofibrillar disorganization with undefined borders. Ultrastructurally, such areas frequently occupied the whole myofibre cross section and extended to a moderate number of sarcomeres in length. Molecular genetic investigations identified recessive mutations in the ryanodine receptor (RYR1) gene in six compound heterozygous patients and one homozygous patient. Nine mutations are novel and four have already been reported either as pathogenic recessive mutations or as changes affecting a residue associated with dominant malignant hyperthermia susceptibility. Only two mutations were located in the C-terminal transmembrane domain whereas the others were distributed throughout the cytoplasmic region of RyR1. CONCLUSION: Our data enlarge the spectrum of RYR1 mutations and highlight their clinical and morphological heterogeneity. A congenital myopathy featuring ptosis and external ophthalmoplegia, concomitant with the novel histopathological phenotype showing fibres with large, poorly delimited areas of myofibrillar disorganization and internal nuclei, is highly suggestive of an RYR1-related congenital myopathy.

Early onset of hypokalaemic periodic paralysis caused by a novel mutation of the CACNA1S gene.

We report a precocious and atypical form of hypokalaemic periodic paralysis, with clinical manifestations at birth and first episodes of paralysis occurring as early as 1 year of age, although onset of this disease usually occurs between 5-35 years. Extensive molecular analysis showed that the disease was caused by a novel de novo p.Arg897Ser mutation in the CACNA1S gene. The mutation mapped to a new region of the protein, the S4 voltage sensing segment of domain III, at odds with previously reported mutations that exclusively affected domains II and IV.

Muscle imaging in dominant core myopathies linked or unlinked to the ryanodine receptor 1 gene.

OBJECTIVE: To characterize the muscle involvement of patients with central core disease (CCD) caused by mutations in the ryanodine receptor 1 gene (RYR1) and to compare these findings with those from patients with core myopathies unlinked to the RYR1 gene. METHODS: We performed a systematic muscular imaging assessment in 11 patients with an RYR1 gene mutation and compared these findings with those of 5 patients from two unrelated families with autosomal dominant core myopathies not linked to RYR1, ACTA1, or MYH7 gene loci. RESULTS: All patients with RYR1 CCD had a characteristic pattern with predominant involvement of the gluteus maximus, adductor magnus, sartorius, vastus intermediolateralis, soleus, and lateral gastrocnemius muscles. In contrast, muscle CT in the first family not linked to RYR1 showed predominant affection of the gluteus minimus and hamstring muscles, whereas the second family presented with predominant involvement of the gluteus minimus, vastus intermediolateralis, tibialis anterior, and medial gastrocnemius muscles. In addition to muscle imaging data, we present detailed information on the clinical and pathologic findings of these novel phenotypes of core myopathies not linked to RYR1. CONCLUSIONS: Our data suggest genetic heterogeneity in autosomal dominant core myopathies and the existence of additional unidentified genes.

[Oculo-cerebro-renal Lowe syndrome: clinical, biochemical and molecular studies in a Moroccan patient]. / Syndrome oculo-cérébro-rénal de Lowe: étude clinique, biochimique et moléculaire chez un patient marocain.

The oculo-cerebro-renal syndrome of Lowe is a rare X-linked disorder, caused by the inositol biphosphate 5-phosphatase deficiency, localized to the Golgi complex. Several mutations were reported in patient's OCRL gene leading to enzyme deficiency. We report a Moroccan case of OCRL syndrome of Lowe with a neo mutation in exon 10. The patient aged of 19 months was referred to our medical centre because of a psychomotor retardation. He had a medical history of eye abnormalities including cataract and bilateral glaucoma, diagnosed when he was 5 weeks old. Cataract has been treated after chirurgical therapy but ocular hypertonia persisted. Physical examination revealed an axial hypotonia and walking difficulties. Laboratory tests revealed a moderate acidosis (20 mmol/L), a slight decrease of serum phosphate level (24 mg/L) and an increased serum phosphatase activity. Further studies showed mild proteinuria, urinary bicarbonates loosing and generalised hyperaminoaciduria. Based on both clinical and biological data, Lowe syndrome has been suggested. In this context, molecular investigation has been performed using dHPLC/sequencing techniques which allow identifying an original mutation c.776T>C (p.Phe259Ser), localized on the exon 10 of the OCRL gene. The mutation was not found in the probant's mother suggesting a neo mutation. Lowe syndrome is a rare hereditary X-linked disorder resulting from a variety of heterogeneous mutations of OCRL gene. Indeed, numerous mutations have been reported, variations were noted concerning their localization as well as their type. To our knowledge, this is the first report of the neo mutation c.776T>C of OCRL gene and the first published case report of the Lowe syndrome in a Moroccan patient.

Clinical and histopathological aspects of central core disease associated and non-associated with RYR1 locus.

We analysed the clinical, histochemical, ultrastructural and genetic data of patients affected by central core disease (CCD) studied during the last 20 years. From a total series of 86 CCD-families, we have identified 46 CCD families with RYR1 mutations (16 autosomal dominant, 8 autosomal recessive, 17 sporadic cases and 5 de novo mutations). Out of the other 40 CCD families, the RyR1 gene was entirely excluded in 7 families, by cDNA sequencing or linkage analysis, indicating a genetic heterogeneity of CCD.

Prenatal diagnosis of DMD in a female foetus affected by Turner syndrome.

OBJECTIVES: We report on a prenatal diagnosis of DMD complicated by a 45,X karyotype that was revealed only in the chorionic villus long-term culture. METHODS: Cytogenetic investigations were performed on both short-term (STC) and long-term cultures (LTC) of the chorionic villus sample. Familial segregation was performed using a panel of intragenic polymorphic markers, and multiplex PCR was used to characterize exonic deletion. RESULTS: Investigations performed for sex determination after STC of the chorionic villus sample showed a normal karyotype 46,XX, while the karyotype performed after LTC revealed a homogeneous monosomy X. Cytogenetic analysis performed on amniotic fluid cells showed 45,X/46,XX mosaicism. Familial segregation analysis for DMD showed loss of heterozygosity for the STR49 marker in the DNA of the proband, her mother and the foetus. Dystrophin gene analysis on the 45,X cells led to the identification of a deletion of exon 50. CONCLUSIONS: The report described a rare situation of monosomy X associated with a DMD genotype. The data confirmed the DMD carrier status of the proband and her mother and indicated that the foetus had a high risk to combine a Turner phenotype and DMD. This study illustrated the potential risk of using short-term culture of villi as the only source of biological material for prenatal diagnosis.

Comparative analysis of in vitro contracture tests with ryanodine and a combination of ryanodine with either halothane or caffeine: a comparative investigation in malignant hyperthermia.

BACKGROUND: The diagnosis of susceptibility to malignant hyperthermia (MH) is currently performed on muscle biopsies subjected to halothane-caffeine in vitro contracture tests (IVCTs). There is a consensus on our need to improve the diagnostic potential of IVCTs if we are to maximize the information available for research and diagnosis in MH. This study was designed as a pilot comparative study and we aimed at comparing the ryanodine test and new tests using a combination of ryanodine, halothane and caffeine. METHODS: One hundred and thirty-two subjects (52 MHS and 80 MHN) were included in this study and new IVCTs were performed in additional muscle biopsy specimens. The contracture time-course was compared considering the onset time of contracture (OT) and the time to reach a 10 mN contracture (10T). Cut-off values were determined using ROC analyses. RESULTS: For the ryanodine test, sensitivity and specificity calculated for OT were 84.6% and 90.4%, respectively, and were better than those obtained using 10T. Combined tests using either caffeine and ryanodine or halothane and ryanodine did provide higher sensitivities (from 85.3 to 93.9%). A better specificity was only observed for the IVC tests combining halothane (cumulated) and caffeine both with ryanodine (93.9% for both). The largest sensitivity was observed when halothane was used as a bolus and combined with ryanodine. The specificity was always larger with the combined tests as compared to the test using ryanodine alone (from 79.1 to 90.9%). This superiority was confirmed, at least in part, when comparing genetic investigations and the results of new tests in a subgroup of subjects. CONCLUSIONS: This pilot study showed a clear diagnostic potential for new IVC tests combining halothane, the triggering agent of MH, and ryanodine acting at the calcium release channel, and should be considered as a first step in the investigation of combined tests.

[Genetic of diseases by abnormal functioning of the skeletal muscle-calcium releasing complex]. / Génétique des pathologies associées à un dysfonctionnement du complexe de mobilisation calcique du muscle squelettique.

Myoplasmic calcium homeostasis is an essential feature of skeletal muscle contraction. The calcium mobilisation complex (CMC) located at the level of the triadic junction plays a major role for the regulation of calcium fluxes between extra-cellular, cytoplasmic and intra-cellular compartments. The ryanodine receptor type I (RYR1), which is located at the level of the terminal cisternae of the sarcoplasmic reticulum is a key component of the CMC. RYR1 allow the release into the myoplasm of the intralumenal stores of calcium. RYR1 interacts with other proteins: DiHydroPyridine Receptor, triadin, calsequestrin, FKBP12, calmodulin. Malignant hyperthermia (MHS) and congenital core myopathies have been associated with a dysfunction of the CMC. MHS is an autosomic dominant pharmacogenetic disease. The MH crisis is induced by exposure of the predisposed patients to halogenated volatile anaesthetics. MHS is characterised by a genetic heterogeneity and two genes, RYR1 and CACNA1S, have been associated so far with the disease. Mutations in the RYR1 gene have been recently associated with heat stroke, a related syndrome. Central Core Disease (CCD) and Multi minicore Disease (MmD) are congenital myopathies presenting with clinical variability and characterized by the presence of specific although heterogeneous muscle histological features: the cores. Clinical boundaries between the two diseases may overlap and the specific diagnosis is often based on the nature of the cores. These diseases show genetic heterogeneity with both autosomic dominant and recessive mode of inheritance and mutations in the SEPN1, RYR1, ACTA1, TPM3 genes have been reported. Mutations associated with MHS were mainly identified into 2 regions of the N-terminal part of RYR1. Functional role of these two domains is still unclear. Mutations responsible for congenital myopathies mainly mapped to the C terminal region of RYR1 that form the transmembrane calcium channel. Functional studies of the RYR1 mutations have shown that MHS mutations were mainly associated with an alteration of the calcium fluxes in response to caffeine or halothane while CCD mutations would result in a leaky RYR1 channel or would alter the Excitation-Contraction coupling at the level of the CMC.

Recent advances in the diagnosis of malignant hyperthermia susceptibility: how confident can we be of genetic testing?

Malignant hyperthermia (MH) is a condition that manifests in susceptible individuals only on exposure to certain anaesthetic agents. Although genetically heterogeneous, mutations in the RYR1 gene (19q13.1) are associated with the majority of reported MH cases. Guidelines for the genetic diagnosis for MH susceptibility have recently been introduced by the European MH Group (EMHG). These are designed to supplement the muscle biopsy testing procedure, the in vitro contracture test (IVCT), which has been the only means of patient screening for the last 30 years and which remains the method for definitive diagnosis in suspected probands. Discordance observed in some families between IVCT phenotype and susceptibility locus genotype could limit the confidence in genetic diagnosis. We have therefore assessed the prevalence of 15 RYR1 mutations currently used in the genetic diagnosis of MH in a sample of over 500 unrelated European MH susceptible individuals and have recorded the frequency of RYR1 genotype/IVCT phenotype discordance. RYR1 mutations were detected in up to approximately 30% of families investigated. Phenotype/genotype discordance in a single individual was observed in 10 out of 196 mutation-positive families. In five families a mutation-positive/IVCT-negative individual was observed, and in the other five families a mutation-negative/IVCT-positive individual was observed. These data represent the most comprehensive assessment of RYR1 mutation prevalence and genotype/phenotype correlation analysis and highlight the possible limitations of MH screening methods. The implications for genetic diagnosis are discussed.

Evaluation and clinical interest of mannan binding lectin function in human plasma.

The mannan binding lectin (MBL) plays a major role in innate immunity through its ability to activate complement upon binding to carbohydrate arrays on the surface of various microorganisms. The question of a possible association of the MBL structural gene polymorphism and the oligomeric state of MBL was poorly documented. For these reasons, it appears difficult to evaluate MBL in blood patients on the only basis of protein contents, even in combination with MBL genotyping. This study reports a method to calculate a specific activity for circulating MBL, that relies on: (i) the availability of purified MBL; and (ii) a simplified MBL activity assay based on complement activation. The three-step MBL purification from human plasma reported here is characterized by a highly purified MBL, that occurs in two different oligomeric forms. The results on the specific activity of these forms show that the higher oligomeric forms of MBL have the ability to induce C4 cleavage more efficiently than the corresponding lower oligomers. The usefulness of this approach is illustrated by its potential interest in the biological exploration of certain pathology, for example in the follow-up of chronic hepatitis C. Further investigation is needed to establish whether MBL specific activity (MBLsa) is correlated to the polymorphic state of the molecule. The relative simplicity of the test described here allows better investigation on the relationship between MBL biological activity and its genotype.

Predicting the risk of cystic fibrosis with abnormal ultrasound signs of fetal bowel: results of a French molecular collaborative study based on 641 prospective cases.

Hyperechogenic fetal bowel is prenatally detected by ultrasound during the second trimester of pregnancy in 0.1-1.8% of fetuses. It has been described as a normal variant but has often been associated with severe diseases, notably cystic fibrosis (CF). The aim of our study was to determine the risk of CF in a prospective study of 641 fetuses with ultrasonographically abnormal fetal bowel and the residual risk when only one mutation is detected in the fetus. Fetal cells and/or parental blood cells were screened for CFTR mutations. Two screening steps were used, the first covering the mutations most frequently observed in French CF patients (mutation detection rate of 70-90%) and, when a CF mutation was detected, a DGGE-sequencing strategy. We observed a 3.1% risk of CF when a digestive tract anomaly was prenatally observed at routine ultrasound examination. The risk was higher when hyperechogenicity was associated with bowel dilatation (5/29; 17%) or with the absence of gall bladder (2/8; 25%). The residual risk of CF was 11% when only one CF mutation was detected by the first screening step, thereby justifying in-depth screening. Mutations associated with severe CF (DeltaF508 mutation) were more frequently observed in these ultrasonographically and prenatally detected CF cases. However, the frequency of heterozygous cases was that observed in the normal population, which demonstrates that heterozygous carriers of CF mutations are not at increased risk for hyperechogenic bowel. In conclusion, fetal bowel anomalies indicate a risk of severe cystic fibrosis and justify careful CFTR molecular analysis.

Familial and sporadic forms of central core disease are associated with mutations in the C-terminal domain of the skeletal muscle ryanodine receptor.

Central core disease (CCD) is an autosomal dominant congenital myopathy. Diagnosis is based on the presence of cores in skeletal muscles. CCD has been linked to the gene encoding the ryanodine receptor (RYR1) and is considered to be an allelic disease of malignant hyperthermia susceptibility. However, the report of a recessive form of transmission together with a variable clinical presentation has raised the question of the genetic heterogeneity of the disease. Analyzing a panel of 34 families exclusively recruited on the basis of both clinically and morphologically expressed CCD, 12 different mutations of the C-terminal domain of RYR1 have been identified in 16 unrelated families. Morphological analysis of the patients' muscles showed different aspects of cores, all of them associated with mutations in the C-terminal region of RYR1. Furthermore, we characterized the presence of neomutations in the RyR1 gene in four families. This indicates that neomutations into the RyR1 gene are not a rare event and must be taken into account for genetic studies of families that present with congenital myopathies type 'central core disease'. Three mutations led to the deletion in frame of amino acids. This is the first report of amino acid deletions in RYR1 associated with CCD. According to a four-transmembrane domain model, the mutations concentrated mostly in the myoplasmic and luminal loops linking, respectively, transmembrane domains T1 and T2 or T3 and T4 of RYR1.