Uniform random generation of large acyclic digraphs

Directed acyclic graphs are the basic representation of the structure underlying Bayesian networks, which represent multivariate probability distributions. In many practical applications, such as the reverse engineering of gene regulatory networks, not only the estimation of model parameters but the reconstruction of the structure itself is of great interest. As well as for the assessment of different structure learning algorithms in simulation studies, a uniform sample from the space of directed acyclic graphs is required to evaluate the prevalence of certain structural features. Here we analyse how to sample acyclic digraphs uniformly at random through recursive enumeration, an approach previously thought too computationally involved. Based on complexity considerations, we discuss in particular how the enumeration directly provides an exact method, which avoids the convergence issues of the alternative Markov chain methods and is actually computationally much faster. The limiting behaviour of the distribution of acyclic digraphs then allows us to sample arbitrarily large graphs. Building on the ideas of recursive enumeration based sampling we also introduce a novel hybrid Markov chain with much faster convergence than current alternatives while still being easy to adapt to various restrictions. Finally we discuss how to include such restrictions in the combinatorial enumeration and the new hybrid Markov chain method for efficient uniform sampling of the corresponding graphs.Comment: 15 pages, 2 figures. To appear in Statistics and Computin

Expanding the phenotype associated to KMT2A variants: overlapping clinical signs between Wiedemann–Steiner and Rubinstein–Taybi syndromes

Lysine-specific methyltransferase 2A (KMT2A) is responsible for methylation of histone H3 (K4H3me) and contributes to chromatin remodeling, acting as \u201cwriter\u201d of the epigenetic machinery. Mutations in KMT2A were first reported in Wiedemann\u2013Steiner syndrome (WDSTS). More recently, KMT2A variants have been described in probands with a specific clinical diagnosis comprised in the so-called chromatinopathies. Such conditions, including WDSTS, are a group of overlapping disorders caused by mutations in genes coding for the epigenetic machinery. Among them, Rubinstein\u2013Taybi syndrome (RSTS) is mainly caused by heterozygous pathogenic variants in CREBBP or EP300. In this work, we used next generation sequencing (either by custom-made panel or by whole exome) to identify alternative causative genes in individuals with a RSTS-like phenotype negative to CREBBP and EP300 mutational screening. In six patients we identified different novel unreported variants in KMT2A gene. The identified variants are de novo in at least four out of six tested individuals and all of them display some typical RSTS phenotypic features but also WDSTS specific signs. This study reinforces the concept that germline variants affecting the epigenetic machinery lead to a shared molecular effect (alteration of the chromatin state) determining superimposable clinical conditions

Mast cells and acute coronary syndromes : Relationship between serum tryptase, clinical outcome and severity of coronary artery disease

Objective: To assess the relationship between serum tryptase and the occurrence of major cardiovascular and cerebrovascular events (MACCE) at 2-year followup in patients admitted with acute coronary syndrome (ACS). To compare serum tryptase to other validated prognostic markers (maximum high-sensitivity troponin (hs-Tn), C reactive protein (CRP) levels at admission, Synergy between percutaneous coronary intervention with Taxus and Cardiac Surgery (SYNTAX) score). Methods: We measured serum tryptase at admission in 140 consecutive patients with ACS and in 50 healthy controls. The patients' follow-up was maintained for 2 years after discharge. The predictive accuracy of serum tryptase for 2-year MACCE was assessed and compared with hs-Tn, CRP and SYNTAX score. Results: Serum tryptase levels at admission were significantly higher in patients with ACS compared with the control group (p=0.0351). 2 years after discharge, 28/140 patients (20%) experienced MACCE. Serum tryptase levels, maximum hs-Tn measurements and SYNTAX score were higher in patients who experienced MACCE compared with those without (p<0.0001). Conversely, we found no significant association between MACCE and CRP. The predictive accuracy of serum tryptase for MACCE was set at the cut-off point of 6.7 ng/mL (sensitivity 46%, specificity 84%). Conclusions: In patients with ACS, serum tryptase measured during index admission is significantly correlated to the development of MACCE up to 2 years, demonstrating a possible long-term prognostic role of this biomarker

Phenotypic overlap of Roberts and Baller-Gerold syndromes in two patients with craniosynostosis, limb reductions, and ESCO2 mutations

Baller-Gerold (BGS, MIM#218600) and Roberts (RBS, MIM#268300) syndromes are rare autosomal recessive disorders caused, respectively, by biallelic alterations in RECQL4 (MIM*603780) and ESCO2 (MIM*609353) genes. Common features are severe growth retardation, limbs shortening and craniofacial abnormalities which may include craniosynostosis. We aimed at unveiling the genetic lesions underpinning the phenotype of two unrelated children with a presumptive BGS diagnosis: patient 1 is a Turkish girl with short stature, microcephaly, craniosynostosis, seizures, intellectual disability, midface hemangioma, bilateral radial and thumb aplasia, tibial hypoplasia, and pes equinovarus. Patient 2 is an Iranian girl born to consanguineous parents with craniosynostosis, micrognathism, bilateral radial aplasia, thumbs, and foot deformity in the context of developmental delay. Upon negative RECQL4 test, whole exome sequencing (WES) analysis performed on the two trios led to the identification of two different ESCO2 homozygous inactivating variants: a previously described c.1131+1G&gt;A transition in patient 1 and an unreported deletion, c.417del, in patient 2, thus turning the diagnosis into Roberts syndrome. The occurrence of a Baller-Gerold phenotype in two unrelated patients that were ultimately diagnosed with RBS demonstrates the strength of WES in redefining the nosological landscape of rare congenital malformation syndromes, a premise to yield optimized patients management and family counseling

Rothmund-Thomson syndrome: Insights from new patients on the genetic variability underpinning clinical presentation and cancer outcome

Biallelic mutations in RECQL4 gene, a caretaker of the genome, cause Rothmund-Thomson type-II syndrome (RTS-II) and confer increased cancer risk if they damage the helicase domain. We describe five families exemplifying clinical and allelic heterogeneity of RTS-II, and report the effect of pathogenic RECQL4 variants by in silico predictions and transcripts analyses. Complete phenotype of patients #39 and #42 whose affected siblings developed osteosarcoma correlates with their c.[1048_1049del], c.[1878+32_1878+55del] and c.[1568G>C;1573delT], c.[3021_3022del] variants which damage the helicase domain. Literature survey highlights enrichment of these variants affecting the helicase domain in patients with cancer outcome raising the issue of strict oncological surveillance. Conversely, patients #29 and #19 have a mild phenotype and carry, respectively, the unreported homozygous c.3265G>T and c.3054A>G variants, both sparing the helicase domain. Finally, despite matching several criteria for RTS clinical diagnosis, patient #38 is heterozygous for c.2412_2414del, no pathogenic CNVs out of those evidenced by high-resolution CGH-array, emerged as contributors to her phenotype

Genetic characterization of cassava (Manihot esculenta) landraces in Brazil assessed with simple sequence repeats

Based on nine microsatellite loci, the aim of this study was to appraise the genetic diversity of 42 cassava (Manihot esculenta) landraces from selected regions in Brazil, and examine how this variety is distributed according to origin in several municipalities in the states of Minas Gerais, São Paulo, Mato Grosso do Sul, Amazonas and Mato Grosso. High diversity values were found among the five above-mentioned regions, with 3.3 alleles per locus on an average, a high percentage of polymorphic loci varying from 88.8% to 100%, an average of 0.265 for observed heterozygosity and 0.570 for gene diversity. Most genetic diversity was concentrated within the regions themselves (HS = 0.52). Cluster analysis and principal component based scatter plotting showed greater similarity among landraces from São Paulo, Mato Grosso do Sul and Amazonas, whereas those from Minas Gerais were clustered into a sub-group within this group. The plants from Mato Grosso, mostly collected in the municipality of General Carneiro, provided the highest differentiation. The migration of human populations is one among the possible reasons for this closer resemblance or greater disparity among plants from the various regions

Differential signature of the centrosomal MARK4 isoforms in glioma

Background: MAP/microtubule affinity-regulating kinase 4 (MARK4) is a serine-threonine kinase expressed in two spliced isoforms, MARK4L and MARK4S, of which MARK4L is a candidate for a role in neoplastic transformation. Methods: We performed mutation analysis to identify sequence alterations possibly affecting MARK4 expression. We then investigated the MARK4L and MARK4S expression profile in 21 glioma cell lines and 36 tissues of different malignancy grades, glioblastoma-derived cancer stem cells (GBM CSCs) and mouse neural stem cells (NSCs) by real-time PCR, immunoblotting and immunohistochemistry. We also analyzed the sub-cellular localisation of MARK4 isoforms in glioma and normal cell lines by immunofluorescence. Results: Mutation analysis rules out sequence variations as the cause of the altered MARK4 expression in glioma. Expression profiling confirms that MARK4L is the predominant isoform, whereas MARK4S levels are significantly decreased in comparison and show an inverse correlation with tumour grade. A high MARK4L/MARK4S ratio also characterizes undifferentiated cells, such as GBM CSCs and NSCs. Accordingly, only MARK4L is expressed in brain neurogenic regions. Moreover, while both MARK4 isoforms are localised to the centrosome and midbody in glioma and normal cells, the L isoform exhibits an additional nucleolar localisation in tumour cells. Conclusions: The observed switch towards MARK4L suggests that the balance between the MARK4 isoforms is carefully guarded during neural differentiation but may be subverted in gliomagenesis. Moreover, the MARK4L nucleolar localisation in tumour cells features this MARK4 isoform as a nucleolus-associated tumour marker

Search for direct production of charginos and neutralinos in events with three leptons and missing transverse momentum in √s = 7 TeV pp collisions with the ATLAS detector

A search for the direct production of charginos and neutralinos in final states with three electrons or muons and missing transverse momentum is presented. The analysis is based on 4.7 fb−1 of proton–proton collision data delivered by the Large Hadron Collider and recorded with the ATLAS detector. Observations are consistent with Standard Model expectations in three signal regions that are either depleted or enriched in Z-boson decays. Upper limits at 95% confidence level are set in R-parity conserving phenomenological minimal supersymmetric models and in simplified models, significantly extending previous results

Jet size dependence of single jet suppression in lead-lead collisions at sqrt(s(NN)) = 2.76 TeV with the ATLAS detector at the LHC

Measurements of inclusive jet suppression in heavy ion collisions at the LHC provide direct sensitivity to the physics of jet quenching. In a sample of lead-lead collisions at sqrt(s) = 2.76 TeV corresponding to an integrated luminosity of approximately 7 inverse microbarns, ATLAS has measured jets with a calorimeter over the pseudorapidity interval |eta| < 2.1 and over the transverse momentum range 38 < pT < 210 GeV. Jets were reconstructed using the anti-kt algorithm with values for the distance parameter that determines the nominal jet radius of R = 0.2, 0.3, 0.4 and 0.5. The centrality dependence of the jet yield is characterized by the jet "central-to-peripheral ratio," Rcp. Jet production is found to be suppressed by approximately a factor of two in the 10% most central collisions relative to peripheral collisions. Rcp varies smoothly with centrality as characterized by the number of participating nucleons. The observed suppression is only weakly dependent on jet radius and transverse momentum. These results provide the first direct measurement of inclusive jet suppression in heavy ion collisions and complement previous measurements of dijet transverse energy imbalance at the LHC.Comment: 15 pages plus author list (30 pages total), 8 figures, 2 tables, submitted to Physics Letters B. All figures including auxiliary figures are available at http://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/HION-2011-02