Bayesian models for syndrome- and gene-specific probabilities of novel variant pathogenicity

BACKGROUND: With the advent of affordable and comprehensive sequencing technologies, access to molecular genetics for clinical diagnostics and research applications is increasing. However, variant interpretation remains challenging, and tools that close the gap between data generation and data interpretation are urgently required. Here we present a transferable approach to help address the limitations in variant annotation. METHODS: We develop a network of Bayesian logistic regression models that integrate multiple lines of evidence to evaluate the probability that a rare variant is the cause of an individual's disease. We present models for genes causing inherited cardiac conditions, though the framework is transferable to other genes and syndromes. RESULTS: Our models report a probability of pathogenicity, rather than a categorisation into pathogenic or benign, which captures the inherent uncertainty of the prediction. We find that gene- and syndrome-specific models outperform genome-wide approaches, and that the integration of multiple lines of evidence performs better than individual predictors. The models are adaptable to incorporate new lines of evidence, and results can be combined with familial segregation data in a transparent and quantitative manner to further enhance predictions. Though the probability scale is continuous, and innately interpretable, performance summaries based on thresholds are useful for comparisons. Using a threshold probability of pathogenicity of 0.9, we obtain a positive predictive value of 0.999 and sensitivity of 0.76 for the classification of variants known to cause long QT syndrome over the three most important genes, which represents sufficient accuracy to inform clinical decision-making. A web tool APPRAISE [http://www.cardiodb.org/APPRAISE] provides access to these models and predictions. CONCLUSIONS: Our Bayesian framework provides a transparent, flexible and robust framework for the analysis and interpretation of rare genetic variants. Models tailored to specific genes outperform genome-wide approaches, and can be sufficiently accurate to inform clinical decision-making

Both doublecortin and doublecortin-like kinase play a role in cortical interneuron migration

Type I lissencephaly, a genetic disease characterized by disorganized cortical layers and gyral abnormalities, is associated with severe cognitive impairment and epilepsy. Two genes, LIS1 and doublecortin (DCX), have been shown to be responsible for a large proportion of cases of type I lissencephaly. Both genes encode microtubule-associated proteins that have been shown to be important for radial migration of cortical pyramidal neurons. To investigate whether DCX also plays a role in cortical interneuron migration, we inactivated DCX in the ganglionic eminence of rat embryonic day 17 brain slices using short hairpin RNA. We found that, when DCX expression was blocked, the migration of interneurons from the ganglionic eminence to the cerebral cortex was slowed but not absent, similar to what had previously been reported for radial neuronal migration. In addition, the processes of DCX-deficient migrating interneurons were more branched than their counterparts in control experiments. These effects were rescued by DCX overexpression, confirming the specificity to DCX inactivation. A similar delay in interneuron migration was observed when Doublecortin-like kinase (DCLK), a microtubule-associated protein related to DCX, was inactivated, although the morphology of the cells was not affected. The importance of these genes in interneuron migration was confirmed by our finding that the cortices of Dcx, Dclk, and Dcx/Dclk mutant mice contained a reduced number of such cells in the cortex and their distribution was different compared with wild-type controls. However, the defect was different for each group of mutant animals, suggesting that DCX and DCLK have distinct roles in cortical interneuron migration

Normal bone physiology, remodelling and its hormonal regulation

The skeleton has structural and locomotor functions, and is a mineral reservoir. Bone turnover by osteoclasts and osteoblasts is a lifelong process, incorporating growth, modelling and remodelling to repair microdamage and access the mineral reservoir. Signalling between bone cells is essential for the co-ordination of these processes. Osteoblasts regulate osteoclast activity through the RANK/RANK ligand/OPG system, and osteocytes regulate osteoblast activity through sclerostin secretion. If resorption and formation are balanced there is no net change in bone mass after each remodelling cycle, but with ageing and some disease states resorption exceeds formation, leading to negative bone balance, decreased bone mass and loss of microstructural integrity. The rate of remodelling is determined by factors including mechanical loading and endocrine influences. The most important endocrine regulator of bone turnover is oestrogen, but other hormones regulating bone metabolism include IGF-1, PTH and gut and adipocyte hormones

Role of electron-phonon coupling and thermal expansion on band gaps, carrier mobility, and interfacial offsets in kesterite thin-film solar cells

The efficiencies of solar cells based on kesterite Cu2ZnSnS4 (CZTS) and Cu2ZnSnSe4 (CZTSe) are limited by a low open-circuit voltage due to high rates of non-radiative electron-hole recombination. To probe the origin of this bottleneck, we calculate the band offset of CZTS(Se) with CdS, confirming a weak spike of 0.1 eV for CZTS/wurtzite-CdS and a strong spike of 0.4 eV for CZTSe/wurtzite-CdS. We also consider the effects of temperature on the band alignment, finding that increasing temperature significantly enhances the spike-type offset. We further resolve an outstanding discrepancy between the measured and calculated phonon frequencies for the kesterites, and use these to estimate the upper limit of electron and hole mobilities based on optic phonon Fröhlich scattering, which uncovers an intrinsic asymmetry with faster (minority carrier) electron mobility.</jats:p

A search for High Mass Stars Forming in Isolation using CORNISH & ATLASGAL

Theoretical models of high mass star formation lie between two extreme scenarios. At one extreme, all the mass comes from an initially gravitationally-bound core. At the other extreme, the majority of the mass comes from cluster scale gas, which lies far outside the initial core boundary. One way to unambiguously show high mass stars can assemble their gas through the former route would be to find a high mass star forming in isolation. Making use of recently available CORNISH and ATLASGAL Galactic plane survey data, we develop sample selection criteria to try and find such an object. From an initial list of approximately 200 sources, we identify the high mass star forming region G13.384+0.064 as the most promising candidate. The region contains a strong radio continuum source, that is powered by an early B-type star. The bolometric luminosity, derived from infrared measurements, is consistent with this. However, sub-millimetre continuum emission, measured in ATLASGAL, as well as dense gas tracers, such as HCO+(3-2) and N2H+(3-2) indicate that there is less than 100 M$_{\odot}$ of material surrounding this star. We conclude that this region is indeed a promising candidate for a high mass star forming in isolation, but that deeper near-IR observations are required to put a stronger constraint on the upper mass limit of young, lower mass stars in the region. Finally, we discuss the challenges facing future studies in proving a given high mass star is forming in isolation

Histiocytoid cardiomyopathy and microphthalmia with linear skin defects syndrome: phenotypes linked by truncating variants in NDUFB11

Variants in NDUFB11, which encodes a structural component of complex I of the mitochondrial respiratory chain (MRC), were recently independently reported to cause histiocytoid cardiomyopathy (histiocytoid CM) and microphthalmia with linear skin defects syndrome (MLS syndrome). Here we report an additional case of histiocytoid CM, which carries a de novo nonsense variant in NDUFB11 (ENST00000276062.8: c.262C > T; p.[Arg88*]) identified using whole-exome sequencing (WES) of a family trio. An identical variant has been previously reported in association with MLS syndrome. The case we describe here lacked the diagnostic features of MLS syndrome, but a detailed clinical comparison of the two cases revealed significant phenotypic overlap. Heterozygous variants in HCCS (which encodes an important mitochondrially targeted protein) and COX7B, which, like NDUFB11, encodes a protein of the MRC, have also previously been identified in MLS syndrome including a case with features of both MLS syndrome and histiocytoid CM. However, a systematic review of WES data from previously published histiocytoid CM cases, alongside four additional cases presented here for the first time, did not identify any variants in these genes. We conclude that NDUFB11 variants play a role in the pathogenesis of both histiocytoid CM and MLS and that these disorders are allelic (genetically related)

Unified force law for granular impact cratering

Experiments on the low-speed impact of solid objects into granular media have been used both to mimic geophysical events and to probe the unusual nature of the granular state of matter. Observations have been interpreted in terms of conflicting stopping forces: product of powers of projectile depth and speed; linear in speed; constant, proportional to the initial impact speed; and proportional to depth. This is reminiscent of high-speed ballistics impact in the 19th and 20th centuries, when a plethora of empirical rules were proposed. To make progress, we developed a means to measure projectile dynamics with 100 nm and 20 us precision. For a 1-inch diameter steel sphere dropped from a wide range of heights into non-cohesive glass beads, we reproduce prior observations either as reasonable approximations or as limiting behaviours. Furthermore, we demonstrate that the interaction between projectile and medium can be decomposed into the sum of velocity-dependent inertial drag plus depth-dependent friction. Thus we achieve a unified description of low-speed impact phenomena and show that the complex response of granular materials to impact, while fundamentally different from that of liquids and solids, can be simply understood

The cometary composition of a protoplanetary disk as revealed by complex cyanides

Observations of comets and asteroids show that the Solar Nebula that spawned our planetary system was rich in water and organic molecules. Bombardment brought these organics to the young Earth's surface, seeding its early chemistry. Unlike asteroids, comets preserve a nearly pristine record of the Solar Nebula composition. The presence of cyanides in comets, including 0.01% of methyl cyanide (CH3CN) with respect to water, is of special interest because of the importance of C-N bonds for abiotic amino acid synthesis. Comet-like compositions of simple and complex volatiles are found in protostars, and can be readily explained by a combination of gas-phase chemistry to form e.g. HCN and an active ice-phase chemistry on grain surfaces that advances complexity[3]. Simple volatiles, including water and HCN, have been detected previously in Solar Nebula analogues - protoplanetary disks around young stars - indicating that they survive disk formation or are reformed in situ. It has been hitherto unclear whether the same holds for more complex organic molecules outside of the Solar Nebula, since recent observations show a dramatic change in the chemistry at the boundary between nascent envelopes and young disks due to accretion shocks[8]. Here we report the detection of CH3CN (and HCN and HC3N) in the protoplanetary disk around the young star MWC 480. We find abundance ratios of these N-bearing organics in the gas-phase similar to comets, which suggests an even higher relative abundance of complex cyanides in the disk ice. This implies that complex organics accompany simpler volatiles in protoplanetary disks, and that the rich organic chemistry of the Solar Nebula was not unique.Comment: Definitive version of the manuscript is published in Nature, 520, 7546, 198, 2015. This is the author's versio

Validation of methods for converting the original Disease Activity Score (DAS) to the DAS28

© The Author(s) 2018.The Disease Activity Score (DAS) is integral in tailoring the clinical management of rheumatoid arthritis (RA) patients and is an important measure in clinical research. Different versions have been developed over the years to improve reliability and ease of use. Combining the original DAS and the newer DAS28 data in both contemporary and historical studies is important for both primary and secondary data analyses. As such, a methodologically robust means of converting the old DAS to the new DAS28 measure would be invaluable. Using data from The Early RA Study (ERAS), a sub-sample of patients with both DAS and DAS28 data were used to develop new regression imputation formulas using the total DAS score (univariate), and using the separate components of the DAS score (multivariate). DAS were transformed to DAS28 using an existing formula quoted in the literature, and the newly developed formulas. Bland and Altman plots were used to compare the transformed DAS with the recorded DAS28 to ascertain levels of agreement. The current transformation formula tended to overestimate the true DAS28 score, particularly at the higher end of the scale. A formula which uses all separate components of the DAS was found to estimate the scores with a higher level of precision. A new formula is proposed that can be used by other early RA cohorts to convert the original DAS to DAS28.Peer reviewedFinal Published versio