Malin Larsson

Network-Based Multiple Sclerosis Pathway Analysis with GWAS Data from 15,000 Cases and 30,000 Controls

Multiple sclerosis (MS) is an inflammatory CNS disease with a substantial genetic component, originally mapped to only the human leukocyte antigen (HLA) region. In the last 5 years, a total of seven genome-wide association studies and one meta-analysis successfully identified 57 non-HLA susceptibility loci. Here, we merged nominal statistical evidence of association and physical evidence of interaction to conduct a protein-interaction-network-based pathway analysis (PINBPA) on two large genetic MS studies comprising a total of 15,317 cases and 29,529 controls. The distribution of nominally significant loci at the gene level matched the patterns of extended linkage disequilibrium in regions of interest. We found that products of genome-wide significantly associated genes are more likely to interact physically and belong to the same or related pathways. We next searched for subnetworks (modules) of genes (and their encoded proteins) enriched with nominally associated loci within each study and identified those modules in common between the two studies. We demonstrate that these modules are more likely to contain genes with bona fide susceptibility variants and, in addition, identify several high-confidence candidates (including BCL10, CD48, REL, TRAF3, and TEC). PINBPA is a powerful approach to gaining further insights into the biology of associated genes and to prioritizing candidates for subsequent genetic studies of complex traits.

Sex-Specific Effects of Adiponectin on Carotid Intima-Media Thickness and Incident Cardiovascular Disease

Background Plasma adiponectin levels have previously been inversely associated with carotid intima-media thickness (IMT), a marker of subclinical atherosclerosis. In this study, we used a sex-stratified Mendelian randomization approach to investigate whether adiponectin has a causal protective influence on IMT. Methods and Results Baseline plasma adiponectin concentration was tested for association with baseline IMT, IMT progression over 30 months, and occurrence of cardiovascular events within 3 years in 3430 participants (women, n =1777; men, n =1653) with high cardiovascular risk but no prevalent disease. Plasma adiponectin levels were inversely associated with baseline mean bifurcation IMT after adjustment for established risk factors (β =−0.018, P<0.001) in men but not in women (β =−0.006, P =0.185; P for interaction =0.061). Adiponectin levels were inversely associated with progression of mean common carotid IMT in men (β =−0.0022, P =0.047), whereas no association was seen in women (0.0007, P =0.475; P for interaction =0.018). Moreover, we observed that adiponectin levels were inversely associated with coronary events in women (hazard ratio 0.57, 95% CI 0.37 to 0.87) but not in men (hazard ratio 0.82, 95% CI 0.54 to 1.25). A gene score of adiponectin-raising alleles in 6 loci, reported recently in a large multi-ethnic meta-analysis, was inversely associated with baseline mean bifurcation IMT in men (β =−0.0008, P =0.004) but not in women (β =−0.0003, P =0.522; P for interaction =0.007). Conclusions This report provides some evidence for adiponectin protecting against atherosclerosis, with effects being confined to men; however, compared with established cardiovascular risk factors, the effect of plasma adiponectin was modest. Further investigation involving mechanistic studies is warranted.

Ebf1 heterozygosity results in increased DNA damage in pro-B cells and their synergistic transformation by Pax5 haploinsufficiency.

Early B-cell factor 1 (Ebf1) is a transcription factor with documented dose-dependent functions in normal and malignant B-lymphocyte development. To understand more about the roles of Ebf1 in malignant transformation, we investigated the impact of reduced functional Ebf1 dosage on mouse B-cell progenitors. Gene expression analysis suggested that Ebf1 was involved in the regulation of genes important for DNA repair and cell survival. Investigation of the DNA damage in steady state, as well as after induction of DNA damage by UV light, confirmed that pro-B cells lacking 1 functional allele of Ebf1 display signs of increased DNA damage. This correlated to reduced expression of DNA repair genes including Rad51, and chromatin immunoprecipitation data suggested that Rad51 is a direct target for Ebf1. Although reduced dosage of Ebf1 did not significantly increase tumor formation in mice, a dramatic increase in the frequency of pro-B cell leukemia was observed in mice with combined heterozygous mutations in the Ebf1 and Pax5 genes, revealing a synergistic effect of combined dose reduction of these proteins. Our data suggest that Ebf1 controls DNA repair in a dose-dependent manner providing a possible explanation to the frequent involvement of EBF1 gene loss in human leukemia.

Early B-cell Factor 1 Regulates the Expansion of B-cell Progenitors in a Dose-dependent Manner

Background: Transcription factor doses play important roles in normal and malignant B-lymphocyte development. Results: We show dose-dependent regulation of B-cell specification and expansion of committed progenitors. Conclusion: Transcription factor dose impacts several aspects of B-cell development. Significance: Knowing the effects of reduced transcription factor dose aids our understanding of the molecular events underlying leukemia and B-cell development. Transcription factor doses are of importance for normal and malignant B-lymphocyte development; however, the understanding of underlying mechanisms and functional consequences of reduced transcription factor levels is limited. We have analyzed progenitor and B-lineage compartments in mice carrying heterozygote mutations in the E2a, Ebf1, or Pax5 gene. Although lymphoid progenitors from Ebf1 or Pax5 heterozygote mice were specified and lineage-restricted in a manner comparable with Wt progenitors, this process was severely impaired in E2a heterozygote mutant mice. This defect was not significantly enhanced upon combined deletion of E2a with Ebf1 or Pax5. Analysis of the pre-B-cell compartment in Ebf1 heterozygote mice revealed a reduction in cell numbers. These cells expressed Pax5 and other B-lineage-associated genes, and global gene expression analysis suggested that the reduction of the pre-B-cell compartment was a result of impaired pre-B-cell expansion. This idea was supported by a reduction in IL2Rα-expressing late pre-B-cells as well as by cell cycle analysis and by the finding that the complexity of the VDJ rearrangement patterns was comparable in Wt and Ebf1+/− pre-B-cells, although the number of progenitors was reduced. Heterozygote deletion of Ebf1 resulted in impaired response to IL7 in vitro and reduced expression levels of pre-BCR on the cell surface, providing possible explanations for the observed stage-specific reduction in cellular expansion. Thus, transcription factor doses are critical for specification as well as expansion of B-lymphoid progenitors, providing increased insight into the molecular regulation of B-cell development.

EPH/ephrin profile and EPHB2 expression predicts patient survival in breast cancer

The EPH and ephrins function as both receptor and ligands and the output on their complex signaling is currently investigated in cancer. Previous work shows that some EPH family members have clinical value in breast cancer, suggesting that this family could be a source of novel clinical targets. Here we quantified the mRNA expression levels of EPH receptors and their ligands, ephrins, in 65 node positive breast cancer samples by RT-PCR with TaqMan® Micro Fluidics Cards Microarray. Upon hierarchical clustering of the mRNA expression levels, we identified a subgroup of patients with high expression, and poor clinical outcome. EPHA2, EPHA4, EFNB1, EFNB2, EPHB2 and EPHB6 were significantly correlated with the cluster groups and particularly EPHB2 was an independent prognostic factor in multivariate analysis and in four public databases. The EPHB2 protein expression was also analyzed by immunohistochemistry in paraffin embedded material (cohort 2). EPHB2 was detected in the membrane and cytoplasmic cell compartments and there was an inverse correlation between membranous and cytoplasmic EPHB2. Membranous EPHB2 predicted longer breast cancer survival in both univariate and multivariate analysis while cytoplasmic EPHB2 indicated shorter breast cancer survival in univariate analysis. Concluding: the EPH/EFN cluster analysis revealed that high EPH/EFN mRNA expression is an independent prognostic factor for poor survival. Especially EPHB2 predicted poor breast cancer survival in several materials and EPHB2 protein expression has also prognostic value depending on cell localization.

Frequent low-level mutations of Protein Kinase D2 in angiolipoma

Tumours displaying differentiation towards normal fat constitute the most common subgroup of soft tissue neoplasms. A series of such tumours was investigated by whole‐exome sequencing followed by targeted ultra‐deep sequencing. Eighty per cent of angiolipomas, but not any other tumour type, displayed mutations in the protein kinase D2 (PRKD2) gene, typically in the part encoding the catalytic domain. The absence of other aberrations at the chromosome or RNA level suggests that PRKD2 mutations are critical for angiolipoma development. Consistently, the mutated PRKD2 alleles were present at low (3–15%) frequencies, indicating that only a subset of the tumour cells is affected. Indeed, by sequencing mature fat cells and other cells separately, the former typically showed the highest mutation frequencies. Thus, we hypothesize that altered PRKD2 signalling in the adipocytic cells drives tumourigenesis and, in agreement with its pivotal role in angiogenesis, induces the vessel formation that is characteristic for angiolipoma. Copyright

Somatic mutagenesis in satellite cells associates with human skeletal muscle aging

Human aging is associated with a decline in skeletal muscle (SkM) function and a reduction in the number and activity of satellite cells (SCs), the resident stem cells. To study the connection between SC aging and muscle impairment, we analyze the whole genome of single SC clones of the leg muscle vastus lateralis from healthy individuals of different ages (21–78 years). We find an accumulation rate of 13 somatic mutations per genome per year, consistent with proliferation of SCs in the healthy adult muscle. SkM-expressed genes are protected from mutations, but aging results in an increase in mutations in exons and promoters, targeting genes involved in SC activity and muscle function. In agreement with SC mutations affecting the whole tissue, we detect a missense mutation in a SC propagating to the muscle. Our results suggest somatic mutagenesis in SCs as a driving force in the age-related decline of SkM function.Aging skeletal muscle shows declining numbers and activity of satellite cells. Here, Franco et al. show that in satellite cells of the human leg muscle vastus lateralis, somatic mutations accumulate with age and that these mutations become enriched in exons and promoters of genes involved in muscle function.

Sarek: A portable workflow for whole-genome sequencing analysis of germline and somatic variants

Summary Whole-genome sequencing (WGS) is a cornerstone of precision medicine, but portable and reproducible open-source workflows for WGS analyses of germline and somatic variants are lacking. We present Sarek, a modular, comprehensive, and easy-to-install workflow, combining a range of software for the identification and annotation of single-nucleotide variants (SNVs), insertion and deletion variants (indels), structural variants, tumor sample heterogeneity, and karyotyping from germline or paired tumor/normal samples. Sarek is implemented in a bioinformatics workflow language (Nextflow) with Docker and Singularity compatible containers, ensuring easy deployment and full reproducibility at any Linux based compute cluster or cloud computing environment. Sarek supports the human reference genomes GRCh37 and GRCh38, and can readily be used both as a core production workflow at sequencing facilities and as a powerful stand-alone tool for individual research groups. Availability Source code and instructions for local installation are available at GitHub (https://github.com/SciLifeLab/Sarek) under the MIT open-source license, and we invite the research community to contribute additional functionality as a collaborative open-source development project.