Statistical Guidance for Experimental Design and Data Analysis of Mutation Detection in Rare Monogenic Mendelian Diseases by Exome Sequencing

Recently, whole-genome sequencing, especially exome sequencing, has successfully led to the identification of causal mutations for rare monogenic Mendelian diseases. However, it is unclear whether this approach can be generalized and effectively applied to other Mendelian diseases with high locus heterogeneity. Moreover, the current exome sequencing approach has limitations such as false positive and false negative rates of mutation detection due to sequencing errors and other artifacts, but the impact of these limitations on experimental design has not been systematically analyzed. To address these questions, we present a statistical modeling framework to calculate the power, the probability of identifying truly disease-causing genes, under various inheritance models and experimental conditions, providing guidance for both proper experimental design and data analysis. Based on our model, we found that the exome sequencing approach is well-powered for mutation detection in recessive, but not dominant, Mendelian diseases with high locus heterogeneity. A disease gene responsible for as low as 5% of the disease population can be readily identified by sequencing just 200 unrelated patients. Based on these results, for identifying rare Mendelian disease genes, we propose that a viable approach is to combine, sequence, and analyze patients with the same disease together, leveraging the statistical framework presented in this work

Loss of the sphingolipid desaturase DEGS1 causes hypomyelinating leukodystrophy.

Sphingolipid imbalance is the culprit in a variety of neurological diseases, some affecting the myelin sheath. We have used whole-exome sequencing in patients with undetermined leukoencephalopathies to uncover the endoplasmic reticulum lipid desaturase DEGS1 as the causative gene in 19 patients from 13 unrelated families. Shared features among the cases include severe motor arrest, early nystagmus, dystonia, spasticity, and profound failure to thrive. MRI showed hypomyelination, thinning of the corpus callosum, and progressive thalamic and cerebellar atrophy, suggesting a critical role of DEGS1 in myelin development and maintenance. This enzyme converts dihydroceramide (DhCer) into ceramide (Cer) in the final step of the de novo biosynthesis pathway. We detected a marked increase of the substrate DhCer and DhCer/Cer ratios in patients' fibroblasts and muscle. Further, we used a knockdown approach for disease modeling in Danio rerio, followed by a preclinical test with the first-line treatment for multiple sclerosis, fingolimod (FTY720, Gilenya). The enzymatic inhibition of Cer synthase by fingolimod, 1 step prior to DEGS1 in the pathway, reduced the critical DhCer/Cer imbalance and the severe locomotor disability, increasing the number of myelinating oligodendrocytes in a zebrafish model. These proof-of-concept results pave the way to clinical translation

Inference of Relationships in Population Data Using Identity-by-Descent and Identity-by-State

It is an assumption of large, population-based datasets that samples are annotated accurately whether they correspond to known relationships or unrelated individuals. These annotations are key for a broad range of genetics applications. While many methods are available to assess relatedness that involve estimates of identity-by-descent (IBD) and/or identity-by-state (IBS) allele-sharing proportions, we developed a novel approach that estimates IBD0, 1, and 2 based on observed IBS within windows. When combined with genome-wide IBS information, it provides an intuitive and practical graphical approach with the capacity to analyze datasets with thousands of samples without prior information about relatedness between individuals or haplotypes. We applied the method to a commonly used Human Variation Panel consisting of 400 nominally unrelated individuals. Surprisingly, we identified identical, parent-child, and full-sibling relationships and reconstructed pedigrees. In two instances non-sibling pairs of individuals in these pedigrees had unexpected IBD2 levels, as well as multiple regions of homozygosity, implying inbreeding. This combined method allowed us to distinguish related individuals from those having atypical heterozygosity rates and determine which individuals were outliers with respect to their designated population. Additionally, it becomes increasingly difficult to identify distant relatedness using genome-wide IBS methods alone. However, our IBD method further identified distant relatedness between individuals within populations, supported by the presence of megabase-scale regions lacking IBS0 across individual chromosomes. We benchmarked our approach against the hidden Markov model of a leading software package (PLINK), showing improved calling of distantly related individuals, and we validated it using a known pedigree from a clinical study. The application of this approach could improve genome-wide association, linkage, heterozygosity, and other population genomics studies that rely on SNP genotype data

A defect in myoblast fusion underlies Carey-Fineman-Ziter syndrome

Multinucleate cellular syncytial formation is a hallmark of skeletal muscle differentiation. Myomaker, encoded by Mymk (Tmem8c), is a well-conserved plasma membrane protein required for myoblast fusion to form multinucleated myotubes in mouse, chick, and zebrafish. Here, we report that autosomal recessive mutations in MYMK (OMIM 615345) cause Carey-Fineman-Ziter syndrome in humans (CFZS; OMIM 254940) by reducing but not eliminating MYMK function. We characterize MYMK-CFZS as a congenital myopathy with marked facial weakness and additional clinical and pathologic features that distinguish it from other congenital neuromuscular syndromes. We show that a heterologous cell fusion assay in vitro and allelic complementation experiments in mymk knockdown and mymk insT/insT zebrafish in vivo can differentiate between MYMK wild type, hypomorphic and null alleles. Collectively, these data establish that MYMK activity is necessary for normal muscle development and maintenance in humans, and expand the spectrum of congenital myopathies to include cell-cell fusion deficits

Recessive mutations in muscle-specific isoforms of FXR1 cause congenital multi-minicore myopathy

FXR1 is an alternatively spliced gene that encodes RNA binding proteins (FXR1P) involved in muscle development. In contrast to other tissues, cardiac and skeletal muscle express two FXR1P isoforms that incorporate an additional exon-15. We report that recessive mutations in this particular exon of FXR1 cause congenital multi-minicore myopathy in humans and mice. Additionally, we show that while Myf5-dependent depletion of all FXR1P isoforms is neonatal lethal, mice carrying mutations in exon-15 display non-lethal myopathies which vary in severity depending on the specific effect of each mutation on the protein

Exome Sequencing Identifies ZNF644 Mutations in High Myopia

Myopia is the most common ocular disorder worldwide, and high myopia in particular is one of the leading causes of blindness. Genetic factors play a critical role in the development of myopia, especially high myopia. Recently, the exome sequencing approach has been successfully used for the disease gene identification of Mendelian disorders. Here we show a successful application of exome sequencing to identify a gene for an autosomal dominant disorder, and we have identified a gene potentially responsible for high myopia in a monogenic form. We captured exomes of two affected individuals from a Han Chinese family with high myopia and performed sequencing analysis by a second-generation sequencer with a mean coverage of 30× and sufficient depth to call variants at ∼97% of each targeted exome. The shared genetic variants of these two affected individuals in the family being studied were filtered against the 1000 Genomes Project and the dbSNP131 database. A mutation A672G in zinc finger protein 644 isoform 1 (ZNF644) was identified as being related to the phenotype of this family. After we performed sequencing analysis of the exons in the ZNF644 gene in 300 sporadic cases of high myopia, we identified an additional five mutations (I587V, R680G, C699Y, 3′UTR+12 C>G, and 3′UTR+592 G>A) in 11 different patients. All these mutations were absent in 600 normal controls. The ZNF644 gene was expressed in human retinal and retinal pigment epithelium (RPE). Given that ZNF644 is predicted to be a transcription factor that may regulate genes involved in eye development, mutation may cause the axial elongation of eyeball found in high myopia patients. Our results suggest that ZNF644 might be a causal gene for high myopia in a monogenic form

Multi-messenger observations of a binary neutron star merger

On 2017 August 17 a binary neutron star coalescence candidate (later designated GW170817) with merger time 12:41:04 UTC was observed through gravitational waves by the Advanced LIGO and Advanced Virgo detectors. The Fermi Gamma-ray Burst Monitor independently detected a gamma-ray burst (GRB 170817A) with a time delay of ~1.7 s with respect to the merger time. From the gravitational-wave signal, the source was initially localized to a sky region of 31 deg2 at a luminosity distance of 40+8-8 Mpc and with component masses consistent with neutron stars. The component masses were later measured to be in the range 0.86 to 2.26 Mo. An extensive observing campaign was launched across the electromagnetic spectrum leading to the discovery of a bright optical transient (SSS17a, now with the IAU identification of AT 2017gfo) in NGC 4993 (at ~40 Mpc) less than 11 hours after the merger by the One- Meter, Two Hemisphere (1M2H) team using the 1 m Swope Telescope. The optical transient was independently detected by multiple teams within an hour. Subsequent observations targeted the object and its environment. Early ultraviolet observations revealed a blue transient that faded within 48 hours. Optical and infrared observations showed a redward evolution over ~10 days. Following early non-detections, X-ray and radio emission were discovered at the transient’s position ~9 and ~16 days, respectively, after the merger. Both the X-ray and radio emission likely arise from a physical process that is distinct from the one that generates the UV/optical/near-infrared emission. No ultra-high-energy gamma-rays and no neutrino candidates consistent with the source were found in follow-up searches. These observations support the hypothesis that GW170817 was produced by the merger of two neutron stars in NGC4993 followed by a short gamma-ray burst (GRB 170817A) and a kilonova/macronova powered by the radioactive decay of r-process nuclei synthesized in the ejecta

Dark Energy Survey Year 1 Results:tomographic cross-correlations between DES galaxies and CMB lensing from SPT+Planck

International audienceWe measure the cross-correlation between redMaGiC galaxies selected from the Dark Energy Survey (DES) year 1 data and gravitational lensing of the cosmic microwave background (CMB) reconstructed from South Pole Telescope (SPT) and Planck data over 1289  deg2. When combining measurements across multiple galaxy redshift bins spanning the redshift range of 0.15<z<0.90, we reject the hypothesis of no correlation at 19.9σ significance. When removing small-scale data points where thermal Sunyaev-Zel’dovich signal and nonlinear galaxy bias could potentially bias our results, the detection significance is reduced to 9.9σ. We perform a joint analysis of galaxy-CMB lensing cross-correlations and galaxy clustering to constrain cosmology, finding Ωm=0.276-0.030+0.029 and S8=σ8Ωm/0.3=0.800-0.094+0.090. We also perform two alternate analyses aimed at constraining only the growth rate of cosmic structure as a function of redshift, finding consistency with predictions from the concordance ΛCDM model. The measurements presented here are part of a joint cosmological analysis that combines galaxy clustering, galaxy lensing and CMB lensing using data from DES, SPT and Planck

Cosmological lensing ratios with DES Y1, SPT and Planck

International audienceCorrelations between tracers of the matter density field and gravitational lensing are sensitive to the evolution of the matter power spectrum and the expansion rate across cosmic time. Appropriately defined ratios of such correlation functions, on the other hand, depend only on the angular diameter distances to the tracer objects and to the gravitational lensing source planes. Because of their simple cosmological dependence, such ratios can exploit available signal-to-noise ratio down to small angular scales, even where directly modelling the correlation functions is difficult. We present a measurement of lensing ratios using galaxy position and lensing data from the Dark Energy Survey, and CMB lensing data from the South Pole Telescope and Planck, obtaining the highest precision lensing ratio measurements to date. Relative to the concordance ΛCDM model, we find a best-fitting lensing ratio amplitude of A = 1.1 ± 0.1. We use the ratio measurements to generate cosmological constraints, focusing on the curvature parameter. We demonstrate that photometrically selected galaxies can be used to measure lensing ratios, and argue that future lensing ratio measurements with data from a combination of LSST and Stage-4 CMB experiments can be used to place interesting cosmological constraints, even after considering the systematic uncertainties associated with photometric redshift and galaxy shear estimation