Genetic and functional implications of an exonic TRIM55 variant in heart failure.

BACKGROUND: To tackle the missing heritability of sporadic heart failure, we screened for novel heart failure-associated genetic variants in the Finnish population and functionally characterized a novel variant in vitro and in vivo. METHODS AND RESULTS: Heart failure-associated variants were screened in genotyping array data of the FINRISK study, consisting of 994 cases and 20,118 controls. Based on logistic regression analysis, a potentially damaging variant in TRIM55 (rs138811034), encoding an E140K variant, was selected for validations. In HL-1 cardiomyocytes, we used CRISPR/Cas9 technology to introduce the variant in the endogenous locus, and additionally TRIM55 wildtype or E140K was overexpressed from plasmid. Functional responses were profiled using whole-genome RNA sequencing, RT-PCR and Western analyses, cell viability and cell cycle assays and cell surface area measurements. In zebrafish embryos, cardiac contractility was measured using videomicroscopy after CRISPR-mediated knockout of trim55a or plasmid overexpression of TRIM55 WT or E140K. Genes related to muscle contraction and cardiac stress were highly regulated in Trim55 E140K/- cardiomyocytes. When compared to the WT/WT cells, the variant cells demonstrated reduced viability, significant hypertrophic response to isoproterenol, p21 protein overexpression and impaired cell cycle progression. In zebrafish embryos, the deletion of trim55a or overexpression of TRIM55 E140K reduced cardiac contractility as compared to embryos with wildtype genotype or overexpression of WT TRIM55, respectively. CONCLUSIONS: A previously uncharacterized TRIM55 E140K variant demonstrated a number of functional implications for cardiomyocyte functions in vitro and in vivo. These findings suggest a novel role for TRIM55 polymorphism in predisposing to heart failure.

Immune cell contexture in the bone marrow tumor microenvironment impacts therapy response in CML.

Increasing evidence suggests that the immune system affects prognosis of chronic myeloid leukemia (CML), but the detailed immunological composition of the leukemia bone marrow (BM) microenvironment is unknown. We aimed to characterize the immune landscape of the CML BM and predict the current treatment goal of tyrosine kinase inhibitor (TKI) therapy, molecular remission 4.0 (MR4.0). Using multiplex immunohistochemistry (mIHC) and automated image analysis, we studied BM tissues of CML patients (n = 56) and controls (n = 14) with a total of 30 immunophenotype markers essential in cancer immunology. CML patients' CD4+ and CD8+ T-cells expressed higher levels of putative exhaustion markers PD1, TIM3, and CTLA4 when compared to control. PD1 expression was higher in BM compared to paired peripheral blood (PB) samples, and decreased during TKI therapy. By combining clinical parameters and immune profiles, low CD4+ T-cell proportion, high proportion of PD1+TIM3-CD8+ T cells, and high PB neutrophil count were most predictive of lower MR4.0 likelihood. Low CD4+ T-cell proportion and high PB neutrophil counts predicted MR4.0 also in a validation cohort (n = 52) analyzed with flow cytometry. In summary, the CML BM is characterized by immune suppression and immune biomarkers predicted MR4.0, thus warranting further testing of immunomodulatory drugs in CML treatment.

Enrichment of low-frequency functional variants revealed by whole-genome sequencing of multiple isolated European populations.

The genetic features of isolated populations can boost power in complex-trait association studies, and an in-depth understanding of how their genetic variation has been shaped by their demographic history can help leverage these advantageous characteristics. Here, we perform a comprehensive investigation using 3,059 newly generated low-depth whole-genome sequences from eight European isolates and two matched general populations, together with published data from the 1000 Genomes Project and UK10K. Sequencing data give deeper and richer insights into population demography and genetic characteristics than genotype-chip data, distinguishing related populations more effectively and allowing their functional variants to be studied more fully. We demonstrate relaxation of purifying selection in the isolates, leading to enrichment of rare and low-frequency functional variants, using novel statistics, DVxy and SVxy. We also develop an isolation-index (Isx) that predicts the overall level of such key genetic characteristics and can thus help guide population choice in future complex-trait association studies.

An Expanded Genome-Wide Association Study of Type 2 Diabetes in Europeans.

To characterize type 2 diabetes (T2D)-associated variation across the allele frequency spectrum, we conducted a meta-analysis of genome-wide association data from 26,676 T2D case and 132,532 control subjects of European ancestry after imputation using the 1000 Genomes multiethnic reference panel. Promising association signals were followed up in additional data sets (of 14,545 or 7,397 T2D case and 38,994 or 71,604 control subjects). We identified 13 novel T2D-associated loci (P < 5 × 10-8), including variants near the GLP2R, GIP, and HLA-DQA1 genes. Our analysis brought the total number of independent T2D associations to 128 distinct signals at 113 loci. Despite substantially increased sample size and more complete coverage of low-frequency variation, all novel associations were driven by common single nucleotide variants. Credible sets of potentially causal variants were generally larger than those based on imputation with earlier reference panels, consistent with resolution of causal signals to common risk haplotypes. Stratification of T2D-associated loci based on T2D-related quantitative trait associations revealed tissue-specific enrichment of regulatory annotations in pancreatic islet enhancers for loci influencing insulin secretion and in adipocytes, monocytes, and hepatocytes for insulin action-associated loci. These findings highlight the predominant role played by common variants of modest effect and the diversity of biological mechanisms influencing T2D pathophysiology.

Whole-genome view of the consequences of a population bottleneck using 2926 genome sequences from Finland and United Kingdom.

Isolated populations with enrichment of variants due to recent population bottlenecks provide a powerful resource for identifying disease-associated genetic variants and genes. As a model of an isolate population, we sequenced the genomes of 1463 Finnish individuals as part of the Sequencing Initiative Suomi (SISu) Project. We compared the genomic profiles of the 1463 Finns to a sample of 1463 British individuals that were sequenced in parallel as part of the UK10K Project. Whereas there were no major differences in the allele frequency of common variants, a significant depletion of variants in the rare frequency spectrum was observed in Finns when comparing the two populations. On the other hand, we observed >2.1 million variants that were twice as frequent among Finns compared with Britons and 800 000 variants that were more than 10 times more frequent in Finns. Furthermore, in Finns we observed a relative proportional enrichment of variants in the minor allele frequency range between 2 and 5% (P<2.2 × 10-16). When stratified by their functional annotations, loss-of-function variants showed the highest proportional enrichment in Finns (P=0.0291). In the non-coding part of the genome, variants in conserved regions (P=0.002) and promoters (P=0.01) were also significantly enriched in the Finnish samples. These functional categories represent the highest a priori power for downstream association studies of rare variants using population isolates.

Crowdsourced assessment of common genetic contribution to predicting anti-TNF treatment response in rheumatoid arthritis.

Rheumatoid arthritis (RA) affects millions world-wide. While anti-TNF treatment is widely used to reduce disease progression, treatment fails in ∼one-third of patients. No biomarker currently exists that identifies non-responders before treatment. A rigorous community-based assessment of the utility of SNP data for predicting anti-TNF treatment efficacy in RA patients was performed in the context of a DREAM Challenge (http://www.synapse.org/RA_Challenge). An open challenge framework enabled the comparative evaluation of predictions developed by 73 research groups using the most comprehensive available data and covering a wide range of state-of-the-art modelling methodologies. Despite a significant genetic heritability estimate of treatment non-response trait (h(2)=0.18, P value=0.02), no significant genetic contribution to prediction accuracy is observed. Results formally confirm the expectations of the rheumatology community that SNP information does not significantly improve predictive performance relative to standard clinical traits, thereby justifying a refocusing of future efforts on collection of other data.

Rare loss-of-function variants in SETD1A are associated with schizophrenia and developmental disorders.

By analyzing the whole-exome sequences of 4,264 schizophrenia cases, 9,343 controls and 1,077 trios, we identified a genome-wide significant association between rare loss-of-function (LoF) variants in SETD1A and risk for schizophrenia (P = 3.3 × 10(-9)). We found only two heterozygous LoF variants in 45,376 exomes from individuals without a neuropsychiatric diagnosis, indicating that SETD1A is substantially depleted of LoF variants in the general population. Seven of the ten individuals with schizophrenia carrying SETD1A LoF variants also had learning difficulties. We further identified four SETD1A LoF carriers among 4,281 children with severe developmental disorders and two more carriers in an independent sample of 5,720 Finnish exomes, both with notable neuropsychiatric phenotypes. Together, our observations indicate that LoF variants in SETD1A cause a range of neurodevelopmental disorders, including schizophrenia. Combining these data with previous common variant evidence, we suggest that epigenetic dysregulation, specifically in the histone H3K4 methylation pathway, is an important mechanism in the pathogenesis of schizophrenia.

Genetic fine mapping and genomic annotation defines causal mechanisms at type 2 diabetes susceptibility loci.

We performed fine mapping of 39 established type 2 diabetes (T2D) loci in 27,206 cases and 57,574 controls of European ancestry. We identified 49 distinct association signals at these loci, including five mapping in or near KCNQ1. 'Credible sets' of the variants most likely to drive each distinct signal mapped predominantly to noncoding sequence, implying that association with T2D is mediated through gene regulation. Credible set variants were enriched for overlap with FOXA2 chromatin immunoprecipitation binding sites in human islet and liver cells, including at MTNR1B, where fine mapping implicated rs10830963 as driving T2D association. We confirmed that the T2D risk allele for this SNP increases FOXA2-bound enhancer activity in islet- and liver-derived cells. We observed allele-specific differences in NEUROD1 binding in islet-derived cells, consistent with evidence that the T2D risk allele increases islet MTNR1B expression. Our study demonstrates how integration of genetic and genomic information can define molecular mechanisms through which variants underlying association signals exert their effects on disease.

Targeted resequencing of the pericentromere of chromosome 2 linked to constitutional delay of growth and puberty.

Constitutional delay of growth and puberty (CDGP) is the most common cause of pubertal delay. CDGP is defined as the proportion of the normal population who experience pubertal onset at least 2 SD later than the population mean, representing 2.3% of all adolescents. While adolescents with CDGP spontaneously enter puberty, they are at risk for short stature, decreased bone mineral density, and psychosocial problems. Genetic factors contribute heavily to the timing of puberty, but the vast majority of CDGP cases remain biologically unexplained, and there is no definitive test to distinguish CDGP from pathological absence of puberty during adolescence. Recently, we published a study identifying significant linkage between a locus at the pericentromeric region of chromosome 2 (chr 2) and CDGP in Finnish families. To investigate this region for causal variation, we sequenced chr 2 between the genomic coordinates of 79-124 Mb (genome build GRCh37) in the proband and affected parent of the 13 families contributing most to this linkage signal. One gene, DNAH6, harbored 6 protein-altering low-frequency variants (< 6% in the Finnish population) in 10 of the CDGP probands. We sequenced an additional 135 unrelated Finnish CDGP subjects and utilized the unique Sequencing Initiative Suomi (SISu) population reference exome set to show that while 5 of these variants were present in the CDGP set, they were also present in the Finnish population at similar frequencies. Additional variants in the targeted region could not be prioritized for follow-up, possibly due to gaps in sequencing coverage or lack of functional knowledge of non-genic genomic regions. Thus, despite having a well-characterized sample collection from a genetically homogeneous population with a large population-based reference sequence dataset, we were unable to pinpoint variation in the linked region predisposing delayed puberty. This study highlights the difficulties of detecting genetic variants under linkage regions for complex traits and suggests that advancements in annotation of gene function and regulatory regions of the genome will be critical for solving the genetic background of complex phenotypes like CDGP.