The penetrance of rare variants in cardiomyopathy-associated genes: A cross-sectional approach to estimating penetrance for secondary findings.

Understanding the penetrance of pathogenic variants identified as secondary findings (SFs) is of paramount importance with the growing availability of genetic testing. We estimated penetrance through large-scale analyses of individuals referred for diagnostic sequencing for hypertrophic cardiomyopathy (HCM; 10,400 affected individuals, 1,332 variants) and dilated cardiomyopathy (DCM; 2,564 affected individuals, 663 variants), using a cross-sectional approach comparing allele frequencies against reference populations (293,226 participants from UK Biobank and gnomAD). We generated updated prevalence estimates for HCM (1:543) and DCM (1:220). In aggregate, the penetrance by late adulthood of rare, pathogenic variants (23% for HCM, 35% for DCM) and likely pathogenic variants (7% for HCM, 10% for DCM) was substantial for dominant cardiomyopathy (CM). Penetrance was significantly higher for variant subgroups annotated as loss of function or ultra-rare and for males compared to females for variants in HCM-associated genes. We estimated variant-specific penetrance for 316 recurrent variants most likely to be identified as SFs (found in 51% of HCM- and 17% of DCM-affected individuals). 49 variants were observed at least ten times (14% of affected individuals) in HCM-associated genes. Median penetrance was 14.6% (±14.4% SD). We explore estimates of penetrance by age, sex, and ancestry and simulate the impact of including future cohorts. This dataset reports penetrance of individual variants at scale and will inform the management of individuals undergoing genetic screening for SFs. While most variants had low penetrance and the costs and harms of screening are unclear, some individuals with highly penetrant variants may benefit from SFs.

Multi-response Mendelian randomization: Identification of shared and distinct exposures for multimorbidity and multiple related disease outcomes.

The existing framework of Mendelian randomization (MR) infers the causal effect of one or multiple exposures on one single outcome. It is not designed to jointly model multiple outcomes, as would be necessary to detect causes of more than one outcome and would be relevant to model multimorbidity or other related disease outcomes. Here, we introduce multi-response Mendelian randomization (MR2), an MR method specifically designed for multiple outcomes to identify exposures that cause more than one outcome or, conversely, exposures that exert their effect on distinct responses. MR2 uses a sparse Bayesian Gaussian copula regression framework to detect causal effects while estimating the residual correlation between summary-level outcomes, i.e., the correlation that cannot be explained by the exposures, and vice versa. We show both theoretically and in a comprehensive simulation study how unmeasured shared pleiotropy induces residual correlation between outcomes irrespective of sample overlap. We also reveal how non-genetic factors that affect more than one outcome contribute to their correlation. We demonstrate that by accounting for residual correlation, MR2 has higher power to detect shared exposures causing more than one outcome. It also provides more accurate causal effect estimates than existing methods that ignore the dependence between related responses. Finally, we illustrate how MR2 detects shared and distinct causal exposures for five cardiovascular diseases in two applications considering cardiometabolic and lipidomic exposures and uncovers residual correlation between summary-level outcomes reflecting known relationships between cardiovascular diseases.

Molecular characterisation of 36 multilocus imprinting disturbance (MLID) patients: a comprehensive approach.

BACKGROUND: Imprinting disorders (ImpDis) comprise diseases which are caused by aberrant regulation of monoallelically and parent-of-origin-dependent expressed genes. A characteristic molecular change in ImpDis patients is aberrant methylation signatures at disease-specific loci, without an obvious DNA change at the specific differentially methylated region (DMR). However, there is a growing number of reports on multilocus imprinting disturbances (MLIDs), i.e. aberrant methylation at different DMRs in the same patient. These MLIDs account for a significant number of patients with specific ImpDis, and several reports indicate a central role of pathogenic maternal effect variants in their aetiology by affecting the maturation of the oocyte and the early embryo. Though several studies on the prevalence and the molecular causes of MLID have been conducted, homogeneous datasets comprising both genomic and methylation data are still lacking. RESULTS: Based on a cohort of 36 MLID patients, we here present both methylation data obtained from next-generation sequencing (NGS, ImprintSeq) approaches and whole-exome sequencing (WES). The compilation of methylation data did not reveal a disease-specific MLID episignature, and a predisposition for the phenotypic modification was not obvious as well. In fact, this lack of epigenotype-phenotype correlation might be related to the mosaic distribution of imprinting defects and their functional relevance in specific tissues. CONCLUSIONS: Due to the higher sensitivity of NGS-based approaches, we suggest that ImprintSeq might be offered at reference centres in case of ImpDis patients with unusual phenotypes but MLID negative by conventional tests. By WES, additional MLID causes than the already known maternal effect variants could not be identified, neither in the patients nor in the maternal exomes. In cases with negative WES results, it is currently unclear to what extent either environmental factors or undetected genetic variants contribute to MLID.

Accurate Measurement of DNA Methylation: Challenges and Bias Correction.

DNA methylation is a key epigenetic modification involved in gene regulation whose contribution to disease susceptibility is still not fully understood. As the cost of genome sequencing technologies continues to drop, it will soon become commonplace to perform genome-wide quantification of DNA methylation at a single base-pair resolution. However, the demand for its accurate quantification might vary across studies. When the scope of the analysis is to detect regions of the genome with different methylation patterns between two or more conditions, e.g., case vs control; treatments vs placebo, accuracy is not crucial. This is the case in epigenome-wide association studies used as genome-wide screening of methylation changes to detect new candidate genes and regions associated with a specific disease or condition. If the aim of the analysis is to use DNA methylation measurements as a biomarker for diseases diagnosis and treatment (Laird, Nat Rev Cancer 3:253-266, 2003; Bock, Epigenomics 1:99-110, 2009), it is instead recommended to produce accurate methylation measurements. Furthermore, if the objective is the detection of DNA methylation in subclonal tumor cell populations or in circulating tumor DNA or in any case of mosaicism, the importance of accuracy becomes critical. The aim of this chapter is to describe the factors that could affect the precise measurement of methylation levels and a recent Bayesian statistical method called MethylCal and its extension that have been proposed to minimize this problem.

ImprintSeq, a novel tool to interrogate DNA methylation at human imprinted regions and diagnose multilocus imprinting disturbance.

PURPOSE: Disruptions of genomic imprinting are associated with congenital imprinting disorders (CIDs) and other disease states, including cancer. CIDs are most often associated with altered methylation at imprinted differentially methylated regions (iDMRs). In some cases, multiple iDMRs are affected causing multilocus imprinting disturbances (MLIDs). The availability of accurate, quantitative, and scalable high-throughput methods to interrogate multiple iDMRs simultaneously would enhance clinical diagnostics and research. METHODS: We report the development of a custom targeted methylation sequencing panel that covered most relevant 63 iDMRs for CIDs and the detection of MLIDs. We tested it in 70 healthy controls and 147 individuals with CIDs. We distinguished loss and gain of methylation per differentially methylated region and classified high and moderate methylation alterations. RESULTS: Across a range of CIDs with a variety of molecular mechanisms, ImprintSeq performed at 98.4% sensitivity, 99.9% specificity, and 99.9% accuracy (when compared with previous diagnostic testing). ImprintSeq was highly sensitive for detecting MLIDs and enabled diagnostic criteria for MLID to be proposed. In a child with extreme MLID profile a probable genetic cause was identified. CONCLUSION: ImprintSeq provides a novel assay for clinical diagnostic and research studies of CIDs, MLIDs, and the role of disordered imprinting in human disease states.

Leveraging Genetic Data to Elucidate the Relationship Between COVID-19 and Ischemic Stroke.

Background The relationship between COVID-19 and ischemic stroke is poorly understood due to potential unmeasured confounding and reverse causation. We aimed to leverage genetic data to triangulate reported associations. Methods and Results Analyses primarily focused on critical COVID-19, defined as hospitalization with COVID-19 requiring respiratory support or resulting in death. Cross-trait linkage disequilibrium score regression was used to estimate genetic correlations of critical COVID-19 with ischemic stroke, other related cardiovascular outcomes, and risk factors common to both COVID-19 and cardiovascular disease (body mass index, smoking and chronic inflammation, estimated using C-reactive protein). Mendelian randomization analysis was performed to investigate whether liability to critical COVID-19 was associated with increased risk of any cardiovascular outcome for which genetic correlation was identified. There was evidence of genetic correlation between critical COVID-19 and ischemic stroke (rg=0.29, false discovery rate [FDR]=0.012), body mass index (rg=0.21, FDR=0.00002), and C-reactive protein (rg=0.20, FDR=0.00035), but no other trait investigated. In Mendelian randomization, liability to critical COVID-19 was associated with increased risk of ischemic stroke (odds ratio [OR] per logOR increase in genetically predicted critical COVID-19 liability 1.03, 95% CI 1.00-1.06, P-value=0.03). Similar estimates were obtained for ischemic stroke subtypes. Consistent estimates were also obtained when performing statistical sensitivity analyses more robust to the inclusion of pleiotropic variants, including multivariable Mendelian randomization analyses adjusting for potential genetic confounding through body mass index, smoking, and chronic inflammation. There was no evidence to suggest that genetic liability to ischemic stroke increased the risk of critical COVID-19. Conclusions These data support that liability to critical COVID-19 is associated with an increased risk of ischemic stroke. The host response predisposing to severe COVID-19 is likely to increase the risk of ischemic stroke, independent of other potentially mitigating risk factors.

A trans locus causes a ribosomopathy in hypertrophic hearts that affects mRNA translation in a protein length-dependent fashion.

BACKGROUND: Little is known about the impact of trans-acting genetic variation on the rates with which proteins are synthesized by ribosomes. Here, we investigate the influence of such distant genetic loci on the efficiency of mRNA translation and define their contribution to the development of complex disease phenotypes within a panel of rat recombinant inbred lines. RESULTS: We identify several tissue-specific master regulatory hotspots that each control the translation rates of multiple proteins. One of these loci is restricted to hypertrophic hearts, where it drives a translatome-wide and protein length-dependent change in translational efficiency, altering the stoichiometric translation rates of sarcomere proteins. Mechanistic dissection of this locus across multiple congenic lines points to a translation machinery defect, characterized by marked differences in polysome profiles and misregulation of the small nucleolar RNA SNORA48. Strikingly, from yeast to humans, we observe reproducible protein length-dependent shifts in translational efficiency as a conserved hallmark of translation machinery mutants, including those that cause ribosomopathies. Depending on the factor mutated, a pre-existing negative correlation between protein length and translation rates could either be enhanced or reduced, which we propose to result from mRNA-specific imbalances in canonical translation initiation and reinitiation rates. CONCLUSIONS: We show that distant genetic control of mRNA translation is abundant in mammalian tissues, exemplified by a single genomic locus that triggers a translation-driven molecular mechanism. Our work illustrates the complexity through which genetic variation can drive phenotypic variability between individuals and thereby contribute to complex disease.

Prioritizing the Role of Major Lipoproteins and Subfractions as Risk Factors for Peripheral Artery Disease.

BACKGROUND: Lipoprotein-related traits have been consistently identified as risk factors for atherosclerotic cardiovascular disease, largely on the basis of studies of coronary artery disease (CAD). The relative contributions of specific lipoproteins to the risk of peripheral artery disease (PAD) have not been well defined. We leveraged large-scale genetic association data to investigate the effects of circulating lipoprotein-related traits on PAD risk. METHODS: Genome-wide association study summary statistics for circulating lipoprotein-related traits were used in the mendelian randomization bayesian model averaging framework to prioritize the most likely causal major lipoprotein and subfraction risk factors for PAD and CAD. Mendelian randomization was used to estimate the effect of apolipoprotein B (ApoB) lowering on PAD risk using gene regions proxying lipid-lowering drug targets. Genes relevant to prioritized lipoprotein subfractions were identified with transcriptome-wide association studies. RESULTS: ApoB was identified as the most likely causal lipoprotein-related risk factor for both PAD (marginal inclusion probability, 0.86; P=0.003) and CAD (marginal inclusion probability, 0.92; P=0.005). Genetic proxies for ApoB-lowering medications were associated with reduced risk of both PAD (odds ratio,0.87 per 1-SD decrease in ApoB [95% CI, 0.84-0.91]; P=9×10-10) and CAD (odds ratio,0.66 [95% CI, 0.63-0.69]; P=4×10-73), with a stronger predicted effect of ApoB lowering on CAD (ratio of effects, 3.09 [95% CI, 2.29-4.60]; P<1×10-6). Extra-small very-low-density lipoprotein particle concentration was identified as the most likely subfraction associated with PAD risk (marginal inclusion probability, 0.91; P=2.3×10-4), whereas large low-density lipoprotein particle concentration was the most likely subfraction associated with CAD risk (marginal inclusion probability, 0.95; P=0.011). Genes associated with extra-small very-low-density lipoprotein particle and large low-density lipoprotein particle concentration included canonical ApoB pathway components, although gene-specific effects were variable. Lipoprotein(a) was associated with increased risk of PAD independently of ApoB (odds ratio, 1.04 [95% CI, 1.03-1.04]; P=1.0×10-33). CONCLUSIONS: ApoB was prioritized as the major lipoprotein fraction causally responsible for both PAD and CAD risk. However, ApoB-lowering drug targets and ApoB-containing lipoprotein subfractions had diverse associations with atherosclerotic cardiovascular disease, and distinct subfraction-associated genes suggest possible differences in the role of lipoproteins in the pathogenesis of PAD and CAD.

EPISPOT: An epigenome-driven approach for detecting and interpreting hotspots in molecular QTL studies.

We present EPISPOT, a fully joint framework which exploits large panels of epigenetic annotations as variant-level information to enhance molecular quantitative trait locus (QTL) mapping. Thanks to a purpose-built Bayesian inferential algorithm, EPISPOT accommodates functional information for both cis and trans actions, including QTL hotspot effects. It effectively couples simultaneous QTL analysis of thousands of genetic variants and molecular traits with hypothesis-free selection of biologically interpretable annotations which directly contribute to the QTL effects. This unified, epigenome-aided learning boosts statistical power and sheds light on the regulatory basis of the uncovered hits; EPISPOT therefore marks an essential step toward improving the challenging detection and functional interpretation of trans-acting genetic variants and hotspots. We illustrate the advantages of EPISPOT in simulations emulating real-data conditions and in a monocyte expression QTL study, which confirms known hotspots and finds other signals, as well as plausible mechanisms of action. In particular, by highlighting the role of monocyte DNase-I sensitivity sites from >150 epigenetic annotations, we clarify the mediation effects and cell-type specificity of major hotspots close to the lysozyme gene. Our approach forgoes the daunting and underpowered task of one-annotation-at-a-time enrichment analyses for prioritizing cis and trans QTL hits and is tailored to any transcriptomic, proteomic, or metabolomic QTL problem. By enabling principled epigenome-driven QTL mapping transcriptome-wide, EPISPOT helps progress toward a better functional understanding of genetic regulation.

Leveraging genetic data to elucidate the relationship between Covid-19 and ischemic stroke.

BACKGROUND: The relationship between coronavirus disease 2019 (Covid-19) and ischemic stroke is poorly defined. We aimed to leverage genetic data to investigate reported associations. METHODS: Genetic association estimates for liability to Covid-19 and cardiovascular traits were obtained from large-scale consortia. Analyses primarily focused on critical Covid-19, defined as hospitalization with Covid-19 requiring respiratory support or resulting in death. Cross-trait linkage disequilibrium score regression was used to estimate genetic correlations of critical Covid-19 with ischemic stroke, other related cardiovascular outcomes, and risk factors common to both Covid-19 and cardiovascular disease (body mass index, smoking and chronic inflammation, estimated using C-reactive protein). Mendelian randomization analysis was performed to investigate whether liability to critical Covid-19 was associated with increased risk of any of the cardiovascular outcomes for which genetic correlation was identified. RESULTS: There was evidence of genetic correlation between critical Covid-19 and ischemic stroke (r g =0.29, FDR p -value=4.65×10 -3 ), body mass index (r g =0.21, FDR- p -value = 6.26×10 -6 ) and C-reactive protein (r g =0.20, FDR- p -value=1.35×10 -4 ), but none of the other considered traits. In Mendelian randomization analysis, liability to critical Covid-19 was associated with increased risk of ischemic stroke (odds ratio [OR] per logOR increase in genetically predicted critical Covid-19 liability 1.03, 95% confidence interval 1.00-1.06, p -value=0.03). Similar estimates were obtained when considering ischemic stroke subtypes. Consistent estimates were also obtained when performing statistical sensitivity analyses more robust to the inclusion of pleiotropic variants, including multivariable Mendelian randomization analyses adjusting for potential genetic confounding through body mass index, smoking and chronic inflammation. There was no evidence to suggest that genetic liability to ischemic stroke increased the risk of critical Covid-19. CONCLUSIONS: These data support that liability to critical Covid-19 is associated with an increased risk of ischemic stroke. The host response predisposing to severe Covid-19 is likely to increase the risk of ischemic stroke, independent of other potentially mitigating risk factors.

A computationally efficient Bayesian seemingly unrelated regressions model for high-dimensional quantitative trait loci discovery.

Our work is motivated by the search for metabolite quantitative trait loci (QTL) in a cohort of more than 5000 people. There are 158 metabolites measured by NMR spectroscopy in the 31-year follow-up of the Northern Finland Birth Cohort 1966 (NFBC66). These metabolites, as with many multivariate phenotypes produced by high-throughput biomarker technology, exhibit strong correlation structures. Existing approaches for combining such data with genetic variants for multivariate QTL analysis generally ignore phenotypic correlations or make restrictive assumptions about the associations between phenotypes and genetic loci. We present a computationally efficient Bayesian seemingly unrelated regressions model for high-dimensional data, with cell-sparse variable selection and sparse graphical structure for covariance selection. Cell sparsity allows different phenotype responses to be associated with different genetic predictors and the graphical structure is used to represent the conditional dependencies between phenotype variables. To achieve feasible computation of the large model space, we exploit a factorisation of the covariance matrix. Applying the model to the NFBC66 data with 9000 directly genotyped single nucleotide polymorphisms, we are able to simultaneously estimate genotype-phenotype associations and the residual dependence structure among the metabolites. The R package BayesSUR with full documentation is available at https://cran.r-project.org/web/packages/BayesSUR/.

Disease-specific variant pathogenicity prediction significantly improves variant interpretation in inherited cardiac conditions.

PURPOSE: Accurate discrimination of benign and pathogenic rare variation remains a priority for clinical genome interpretation. State-of-the-art machine learning variant prioritization tools are imprecise and ignore important parameters defining gene-disease relationships, e.g., distinct consequences of gain-of-function versus loss-of-function variants. We hypothesized that incorporating disease-specific information would improve tool performance. METHODS: We developed a disease-specific variant classifier, CardioBoost, that estimates the probability of pathogenicity for rare missense variants in inherited cardiomyopathies and arrhythmias. We assessed CardioBoost's ability to discriminate known pathogenic from benign variants, prioritize disease-associated variants, and stratify patient outcomes. RESULTS: CardioBoost has high global discrimination accuracy (precision recall area under the curve [AUC] 0.91 for cardiomyopathies; 0.96 for arrhythmias), outperforming existing tools (4-24% improvement). CardioBoost obtains excellent accuracy (cardiomyopathies 90.2%; arrhythmias 91.9%) for variants classified with >90% confidence, and increases the proportion of variants classified with high confidence more than twofold compared with existing tools. Variants classified as disease-causing are associated with both disease status and clinical severity, including a 21% increased risk (95% confidence interval [CI] 11-29%) of severe adverse outcomes by age 60 in patients with hypertrophic cardiomyopathy. CONCLUSIONS: A disease-specific variant classifier outperforms state-of-the-art genome-wide tools for rare missense variants in inherited cardiac conditions ( https://www.cardiodb.org/cardioboost/ ), highlighting broad opportunities for improved pathogenicity prediction through disease specificity.

High-throughput multivariable Mendelian randomization analysis prioritizes apolipoprotein B as key lipid risk factor for coronary artery disease.

BACKGROUND: Genetic variants can be used to prioritize risk factors as potential therapeutic targets via Mendelian randomization (MR). An agnostic statistical framework using Bayesian model averaging (MR-BMA) can disentangle the causal role of correlated risk factors with shared genetic predictors. Here, our objective is to identify lipoprotein measures as mediators between lipid-associated genetic variants and coronary artery disease (CAD) for the purpose of detecting therapeutic targets for CAD. METHODS: As risk factors we consider 30 lipoprotein measures and metabolites derived from a high-throughput metabolomics study including 24 925 participants. We fit multivariable MR models of genetic associations with CAD estimated in 453 595 participants (including 113 937 cases) regressed on genetic associations with the risk factors. MR-BMA assigns to each combination of risk factors a model score quantifying how well the genetic associations with CAD are explained. Risk factors are ranked by their marginal score and selected using false-discovery rate (FDR) criteria. We perform supplementary and sensitivity analyses varying the dataset for genetic associations with CAD. RESULTS: In the main analysis, the top combination of risk factors ranked by the model score contains apolipoprotein B (ApoB) only. ApoB is also the highest ranked risk factor with respect to the marginal score (FDR <0.005). Additionally, ApoB is selected in all sensitivity analyses. No other measure of cholesterol or triglyceride is consistently selected otherwise. CONCLUSIONS: Our agnostic genetic investigation prioritizes ApoB across all datasets considered, suggesting that ApoB, representing the total number of hepatic-derived lipoprotein particles, is the primary lipid determinant of CAD.

Leptin-Mediated Changes in the Human Metabolome.

CONTEXT: While severe obesity due to congenital leptin deficiency is rare, studies in patients before and after treatment with leptin can provide unique insights into the role that leptin plays in metabolic and endocrine function. OBJECTIVE: The aim of this study was to characterize changes in peripheral metabolism in people with congenital leptin deficiency undergoing leptin replacement therapy, and to investigate the extent to which these changes are explained by reduced caloric intake. DESIGN: Ultrahigh performance liquid chromatography-tandem mass spectroscopy (UPLC-MS/MS) was used to measure 661 metabolites in 6 severely obese people with congenital leptin deficiency before, and within 1 month after, treatment with recombinant leptin. Data were analyzed using unsupervised and hypothesis-driven computational approaches and compared with data from a study of acute caloric restriction in healthy volunteers. RESULTS: Leptin replacement was associated with class-wide increased levels of fatty acids and acylcarnitines and decreased phospholipids, consistent with enhanced lipolysis and fatty acid oxidation. Primary and secondary bile acids increased after leptin treatment. Comparable changes were observed after acute caloric restriction. Branched-chain amino acids and steroid metabolites decreased after leptin, but not after acute caloric restriction. Individuals with severe obesity due to leptin deficiency and other genetic obesity syndromes shared a metabolomic signature associated with increased BMI. CONCLUSION: Leptin replacement was associated with changes in lipolysis and substrate utilization that were consistent with negative energy balance. However, leptin's effects on branched-chain amino acids and steroid metabolites were independent of reduced caloric intake and require further exploration.

Sphingolipid, fatty acid and phospholipid metabolites are associated with disease severity and mTOR inhibition in lymphangioleiomyomatosis.

BACKGROUND: Lymphangioleiomyomatosis (LAM) is a rare multisystem disease almost exclusively affecting women which causes loss of lung function, lymphatic abnormalities and angiomyolipomas. LAM occurs sporadically and in people with tuberous sclerosis complex (TSC). Loss of TSC gene function leads to dysregulated mechanistic target of rapamycin (mTOR) signalling. As mTOR is a regulator of lipid and nucleotide synthesis, we hypothesised that the serum metabolome would be altered in LAM and related to disease severity and activity. METHODS: Ultrahigh performance liquid chromatography-tandem mass spectroscopy was used to examine the serum metabolome of 79 closely phenotyped women with LAM, including 29 receiving treatment with an mTOR inhibitor and 43 healthy control women. RESULTS: Sphingolipid, fatty acid and phospholipid metabolites were associated with FEV1 in women with LAM (eg, behenoyl sphingomyelin adjusted (adj.) p=8.10 × 10-3). Those with higher disease-burden scores had abnormalities in fatty acid, phospholipid and lysolipids. Rate of loss of FEV1 was associated with differences in acyl-carnitine, acyl-glycines, acyl-glutamine, fatty acids, endocanbinoids and sphingolipids (eg, myristoleoylcarnitine adj. p=0.07). In TSC-LAM, rapamycin affected modules of interrelated metabolites which comprised linoleic acid, the tricarboxylic acid cycle, aminoacyl-tRNA biosynthesis, cysteine, methionine, arginine and proline metabolism. Metabolomic pathway analysis within modules reiterated the importance of glycerophospholipid metabolites (adj. p=0.047). CONCLUSIONS: Women with LAM have altered lipid metabolism. The associations between these metabolites, multiple markers of disease activity and their potential biological roles in cell survival and signalling, suggest that lipid species may be both disease-relevant biomarkers and potential therapeutic targets for LAM.

Metabolic effects of bezafibrate in mitochondrial disease.

Mitochondrial disorders affect 1/5,000 and have no cure. Inducing mitochondrial biogenesis with bezafibrate improves mitochondrial function in animal models, but there are no comparable human studies. We performed an open-label observational experimental medicine study of six patients with mitochondrial myopathy caused by the m.3243A>G MTTL1 mutation. Our primary aim was to determine the effects of bezafibrate on mitochondrial metabolism, whilst providing preliminary evidence of safety and efficacy using biomarkers. The participants received 600-1,200 mg bezafibrate daily for 12 weeks. There were no clinically significant adverse events, and liver function was not affected. We detected a reduction in the number of complex IV-immunodeficient muscle fibres and improved cardiac function. However, this was accompanied by an increase in serum biomarkers of mitochondrial disease, including fibroblast growth factor 21 (FGF-21), growth and differentiation factor 15 (GDF-15), plus dysregulation of fatty acid and amino acid metabolism. Thus, although potentially beneficial in short term, inducing mitochondrial biogenesis with bezafibrate altered the metabolomic signature of mitochondrial disease, raising concerns about long-term sequelae.

A Global-Local Approach for Detecting Hotspots in Multiple-Response Regression.

We tackle modelling and inference for variable selection in regression problems with many predictors and many responses. We focus on detecting hotspots, that is, predictors associated with several responses. Such a task is critical in statistical genetics, as hotspot genetic variants shape the architecture of the genome by controlling the expression of many genes and may initiate decisive functional mechanisms underlying disease endpoints. Existing hierarchical regression approaches designed to model hotspots suffer from two limitations: their discrimination of hotspots is sensitive to the choice of top-level scale parameters for the propensity of predictors to be hotspots, and they do not scale to large predictor and response vectors, for example, of dimensions 103-105 in genetic applications. We address these shortcomings by introducing a flexible hierarchical regression framework that is tailored to the detection of hotspots and scalable to the above dimensions. Our proposal implements a fully Bayesian model for hotspots based on the horseshoe shrinkage prior. Its global-local formulation shrinks noise globally and, hence, accommodates the highly sparse nature of genetic analyses while being robust to individual signals, thus leaving the effects of hotspots unshrunk. Inference is carried out using a fast variational algorithm coupled with a novel simulated annealing procedure that allows efficient exploration of multimodal distributions.

Author Correction: WWP2 regulates pathological cardiac fibrosis by modulating SMAD2 signaling.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

WWP2 regulates pathological cardiac fibrosis by modulating SMAD2 signaling.

Cardiac fibrosis is a final common pathology in inherited and acquired heart diseases that causes cardiac electrical and pump failure. Here, we use systems genetics to identify a pro-fibrotic gene network in the diseased heart and show that this network is regulated by the E3 ubiquitin ligase WWP2, specifically by the WWP2-N terminal isoform. Importantly, the WWP2-regulated pro-fibrotic gene network is conserved across different cardiac diseases characterized by fibrosis: human and murine dilated cardiomyopathy and repaired tetralogy of Fallot. Transgenic mice lacking the N-terminal region of the WWP2 protein show improved cardiac function and reduced myocardial fibrosis in response to pressure overload or myocardial infarction. In primary cardiac fibroblasts, WWP2 positively regulates the expression of pro-fibrotic markers and extracellular matrix genes. TGFß1 stimulation promotes nuclear translocation of the WWP2 isoforms containing the N-terminal region and their interaction with SMAD2. WWP2 mediates the TGFß1-induced nucleocytoplasmic shuttling and transcriptional activity of SMAD2.

MethylCal: Bayesian calibration of methylation levels.

Bisulfite amplicon sequencing has become the primary choice for single-base methylation quantification of multiple targets in parallel. The main limitation of this technology is a preferential amplification of an allele and strand in the PCR due to methylation state. This effect, known as 'PCR bias', causes inaccurate estimation of the methylation levels and calibration methods based on standard controls have been proposed to correct for it. Here, we present a Bayesian calibration tool, MethylCal, which can analyse jointly all CpGs within a CpG island (CGI) or a Differentially Methylated Region (DMR), avoiding 'one-at-a-time' CpG calibration. This enables more precise modeling of the methylation levels observed in the standard controls. It also provides accurate predictions of the methylation levels not considered in the controlled experiment, a feature that is paramount in the derivation of the corrected methylation degree. We tested the proposed method on eight independent assays (two CpG islands and six imprinting DMRs) and demonstrated its benefits, including the ability to detect outliers. We also evaluated MethylCal's calibration in two practical cases, a clinical diagnostic test on 18 patients potentially affected by Beckwith-Wiedemann syndrome, and 17 individuals with celiac disease. The calibration of the methylation levels obtained by MethylCal allows a clearer identification of patients undergoing loss or gain of methylation in borderline cases and could influence further clinical or treatment decisions.