Biallelic TET2 mutations confer sensitivity to 5'-azacitidine in acute myeloid leukemia.

Precision medicine can significantly improve outcomes for patients with cancer, but implementation requires comprehensive characterization of tumor cells to identify therapeutically exploitable vulnerabilities. Here, we describe somatic biallelic TET2 mutations in an elderly patient with acute myeloid leukemia (AML) that was chemoresistant to anthracycline and cytarabine but acutely sensitive to 5'-azacitidine (5'-Aza) hypomethylating monotherapy, resulting in long-term morphological remission. Given the role of TET2 as a regulator of genomic methylation, we hypothesized that mutant TET2 allele dosage affects response to 5'-Aza. Using an isogenic cell model system and an orthotopic mouse xenograft, we demonstrate that biallelic TET2 mutations confer sensitivity to 5'-Aza compared with cells with monoallelic mutations. Our data argue in favor of using hypomethylating agents for chemoresistant disease or as first-line therapy in patients with biallelic TET2-mutated AML and demonstrate the importance of considering mutant allele dosage in the implementation of precision medicine for patients with cancer.

Genome-wide association study identifies susceptibility loci for acute myeloid leukemia.

Acute myeloid leukemia (AML) is a hematological malignancy with an undefined heritable risk. Here we perform a meta-analysis of three genome-wide association studies, with replication in a fourth study, incorporating a total of 4018 AML cases and 10488 controls. We identify a genome-wide significant risk locus for AML at 11q13.2 (rs4930561; P = 2.15 × 10-8; KMT5B). We also identify a genome-wide significant risk locus for the cytogenetically normal AML sub-group (N = 1287) at 6p21.32 (rs3916765; P = 1.51 × 10-10; HLA). Our results inform on AML etiology and identify putative functional genes operating in histone methylation (KMT5B) and immune function (HLA).

Genome-wide association study identifies risk loci for progressive chronic lymphocytic leukemia.

Prognostication in patients with chronic lymphocytic leukemia (CLL) is challenging due to heterogeneity in clinical course. We hypothesize that constitutional genetic variation affects disease progression and could aid prognostication. Pooling data from seven studies incorporating 842 cases identifies two genomic locations associated with time from diagnosis to treatment, including 10q26.13 (rs736456, hazard ratio (HR) = 1.78, 95% confidence interval (CI) = 1.47-2.15; P = 2.71 × 10-9) and 6p (rs3778076, HR = 1.99, 95% CI = 1.55-2.55; P = 5.08 × 10-8), which are particularly powerful prognostic markers in patients with early stage CLL otherwise characterized by low-risk features. Expression quantitative trait loci analysis identifies putative functional genes implicated in modulating B-cell receptor or innate immune responses, key pathways in CLL pathogenesis. In this work we identify rs736456 and rs3778076 as prognostic in CLL, demonstrating that disease progression is determined by constitutional genetic variation as well as known somatic drivers.

In silico engineering of Pseudomonas metabolism reveals new biomarkers for increased biosurfactant production.

BACKGROUND: Rhamnolipids, biosurfactants with a wide range of biomedical applications, are amphiphilic molecules produced on the surfaces of or excreted extracellularly by bacteria including Pseudomonas aeruginosa. However, Pseudomonas putida is a non-pathogenic model organism with greater metabolic versatility and potential for industrial applications. METHODS: We investigate in silico the metabolic capabilities of P. putida for rhamnolipids biosynthesis using statistical, metabolic and synthetic engineering approaches after introducing key genes (RhlA and RhlB) from P. aeruginosa into a genome-scale model of P. putida. This pipeline combines machine learning methods with multi-omic modelling, and drives the engineered P. putida model toward an optimal production and export of rhamnolipids out of the membrane. RESULTS: We identify a substantial increase in synthesis of rhamnolipids by the engineered model compared to the control model. We apply statistical and machine learning techniques on the metabolic reaction rates to identify distinct features on the structure of the variables and individual components driving the variation of growth and rhamnolipids production. We finally provide a computational framework for integrating multi-omics data and identifying latent pathways and genes for the production of rhamnolipids in P. putida. CONCLUSIONS: We anticipate that our results will provide a versatile methodology for integrating multi-omics data for topological and functional analysis of P. putida toward maximization of biosurfactant production.

Genome-wide association analysis of chronic lymphocytic leukaemia, Hodgkin lymphoma and multiple myeloma identifies pleiotropic risk loci.

B-cell malignancies (BCM) originate from the same cell of origin, but at different maturation stages and have distinct clinical phenotypes. Although genetic risk variants for individual BCMs have been identified, an agnostic, genome-wide search for shared genetic susceptibility has not been performed. We explored genome-wide association studies of chronic lymphocytic leukaemia (CLL, N = 1,842), Hodgkin lymphoma (HL, N = 1,465) and multiple myeloma (MM, N = 3,790). We identified a novel pleiotropic risk locus at 3q22.2 (NCK1, rs11715604, P = 1.60 × 10-9) with opposing effects between CLL (P = 1.97 × 10-8) and HL (P = 3.31 × 10-3). Eight established non-HLA risk loci showed pleiotropic associations. Within the HLA region, Ser37 + Phe37 in HLA-DRB1 (P = 1.84 × 10-12) was associated with increased CLL and HL risk (P = 4.68 × 10-12), and reduced MM risk (P = 1.12 × 10-2), and Gly70 in HLA-DQB1 (P = 3.15 × 10-10) showed opposing effects between CLL (P = 3.52 × 10-3) and HL (P = 3.41 × 10-9). By integrating eQTL, Hi-C and ChIP-seq data, we show that the pleiotropic risk loci are enriched for B-cell regulatory elements, as well as an over-representation of binding of key B-cell transcription factors. These data identify shared biological pathways influencing the development of CLL, HL and MM. The identification of these risk loci furthers our understanding of the aetiological basis of BCMs.

Using population data for assessing next-generation sequencing performance.

MOTIVATION: During the past 4 years, whole-exome sequencing has become a standard tool for finding rare variants causing Mendelian disorders. In that time, there has also been a proliferation of both sequencing platforms and approaches to analyse their output. This requires approaches to assess the performance of different methods. Traditionally, criteria such as comparison with microarray data or a number of known polymorphic sites have been used. Here we expand such approaches, developing a maximum likelihood framework and using it to estimate the sensitivity and specificity of whole-exome sequencing data. RESULTS: Using whole-exome sequencing data for a panel of 19 individuals, we show that estimated sensitivity and specificity are similar to those calculated using microarray data as a reference. We explore the effect of frequency misspecification arising from using an inappropriately selected population and find that, although the estimates are affected, the rankings across procedures remain the same. AVAILABILITY AND IMPLEMENTATION: An implementation using Perl and R can be found at busso.ncl.ac.uk (Username: igm101; Password: Z1z1nts).

Functionally significant, rare transcription factor variants in tetralogy of Fallot.

OBJECTIVE: Rare variants in certain transcription factors involved in cardiac development cause Mendelian forms of congenital heart disease. The purpose of this study was to systematically assess the frequency of rare transcription factor variants in sporadic patients with the cardiac outflow tract malformation tetralogy of Fallot (TOF). METHODS AND RESULTS: We sequenced the coding, 5'UTR, and 3'UTR regions of twelve transcription factor genes implicated in cardiac outflow tract development (NKX2.5, GATA4, ISL1, TBX20, MEF2C, BOP/SMYD1, HAND2, FOXC1, FOXC2, FOXH, FOXA2 and TBX1) in 93 non-syndromic, non-Mendelian TOF cases. We also analysed Illumina Human 660W-Quad SNP Array data for copy number variants in these genes; none were detected. Four of the rare variants detected have previously been shown to affect transactivation in in vitro reporter assays: FOXC1 p.P297S, FOXC2 p.Q444R, FOXH1 p.S113T and TBX1 p.P43_G61del PPPPRYDPCAAAAPGAPGP. Two further rare variants, HAND2 p.A25_A26insAA and FOXC1 p.G378_G380delGGG, A488_491delAAAA, affected transactivation in in vitro reporter assays. Each of these six functionally significant variants was present in a single patient in the heterozygous state; each of the four for which parental samples were available were maternally inherited. Thus in the 93 TOF cases we identified six functionally significant mutations in the secondary heart field transcriptional network. SIGNIFICANCE: This study indicates that rare genetic variants in the secondary heart field transcriptional network with functional effects on protein function occur in 3-13% of patients with TOF. This is the first report of a functionally significant HAND2 mutation in a patient with congenital heart disease.

Mutations in FAM111B cause hereditary fibrosing poikiloderma with tendon contracture, myopathy, and pulmonary fibrosis.

Congenital poikiloderma is characterized by a combination of mottled pigmentation, telangiectasia, and epidermal atrophy in the first few months of life. We have previously described a South African European-descent family affected by a rare autosomal-dominant form of hereditary fibrosing poikiloderma accompanied by tendon contracture, myopathy, and pulmonary fibrosis. Here, we report the identification of causative mutations in FAM111B by whole-exome sequencing. In total, three FAM111B missense mutations were identified in five kindreds of different ethnic backgrounds. The mutation segregated with the disease in one large pedigree, and mutations were de novo in two other pedigrees. All three mutations were absent from public databases and were not observed on Sanger sequencing of 388 ethnically matched control subjects. The three single-nucleotide mutations code for amino acid changes that are clustered within a putative trypsin-like cysteine/serine peptidase domain of FAM111B. These findings provide evidence of the involvement of FAM111B in congenital poikiloderma and multisystem fibrosis.

Association between C677T polymorphism of methylene tetrahydrofolate reductase and congenital heart disease: meta-analysis of 7697 cases and 13,125 controls.

BACKGROUND: Association between the C677T polymorphism of the methylene tetrahydrofolate reductase (MTHFR) gene and congenital heart disease (CHD) is contentious. METHODS AND RESULTS: We compared genotypes between CHD cases and controls and between mothers of CHD cases and controls. We placed our results in context by conducting meta-analyses of previously published studies. Among 5814 cases with primary genotype data and 10 056 controls, there was no evidence of association between MTHFR C677T genotype and CHD risk (odds ratio [OR], 0.96 [95% confidence interval, 0.87-1.07]). A random-effects meta-analysis of all studies (involving 7697 cases and 13 125 controls) suggested the presence of association (OR, 1.25 [95% confidence interval, 1.03-1.51]; P=0.022) but with substantial heterogeneity among contributing studies (I(2)=64.4%) and evidence of publication bias. Meta-analysis of large studies only (defined by a variance of the log OR <0.05), which together contributed 83% of all cases, yielded no evidence of association (OR, 0.97 [95% confidence interval, 0.91-1.03]) without significant heterogeneity (I(2)=0). Moreover, meta-analysis of 1781 mothers of CHD cases (829 of whom were genotyped in this study) and 19 861 controls revealed no evidence of association between maternal C677T genotype and risk of CHD in offspring (OR, 1.13 [95% confidence interval, 0.87-1.47]). There was no significant association between MTHFR genotype and CHD risk in large studies from regions with different levels of dietary folate. CONCLUSIONS: The MTHFR C677T polymorphism, which directly influences plasma folate levels, is not associated with CHD risk. Publication biases appear to substantially contaminate the literature with regard to this genetic association.

Genome-wide association study of multiple congenital heart disease phenotypes identifies a susceptibility locus for atrial septal defect at chromosome 4p16.

We carried out a genome-wide association study (GWAS) of congenital heart disease (CHD). Our discovery cohort comprised 1,995 CHD cases and 5,159 controls and included affected individuals from each of the 3 major clinical CHD categories (with septal, obstructive and cyanotic defects). When all CHD phenotypes were considered together, no region achieved genome-wide significant association. However, a region on chromosome 4p16, adjacent to the MSX1 and STX18 genes, was associated (P = 9.5 × 10⁻7) with the risk of ostium secundum atrial septal defect (ASD) in the discovery cohort (N = 340 cases), and this association was replicated in a further 417 ASD cases and 2,520 controls (replication P = 5.0 × 10⁻5; odds ratio (OR) in replication cohort = 1.40, 95% confidence interval (CI) = 1.19-1.65; combined P = 2.6 × 10⁻¹°). Genotype accounted for ~9% of the population-attributable risk of ASD.

Common variation neighbouring micro-RNA 22 is associated with increased left ventricular mass.

AIMS: Previous genome-wide linkage analysis has suggested that chromosomal region 17p13.3 may harbour genes influencing left ventricular mass (LVM) in man. To date, the genetic factors accounting for LVM variability remain largely unknown but a non-coding RNA gene within this region, micro-RNA 22 (miR-22), has been implicated in cardiac hypertrophy and heart failure in animal models. We thus investigated the relationship between common genetic polymorphisms surrounding miR-22 and left ventricular mass in a family-based association study. METHODS AND RESULTS: We studied a cohort of 255 families comprising 1,425 individuals ascertained via a hypertensive proband. Ten single nucleotide polymorphisms which together tagged common genetic variation surrounding the miR-22 gene were genotyped. There was evidence of association between the rs7223247 polymorphism, which lies within the 3'UTR of a gene of unknown function, TLCD2, immediately downstream from miR-22, and left ventricular mass determined by Sokolow-Lyon voltage (Bonferroni corrected p-value = 0.038). The T allele at rs7223247 was associated with an 0.272 standard deviation higher Sokolow-Lyon voltage. Genotype was responsible for ~1% of the population variability in LVM. CONCLUSIONS: Genotype at the rs7223247 polymorphism affects left ventricular mass determined by Sokolow-Lyon voltage. The neighbouring genes miR-22 and TLCD2 are strong candidates to account for this observation.

Genome-wide association study identifies loci on 12q24 and 13q32 associated with tetralogy of Fallot.

We conducted a genome-wide association study to search for risk alleles associated with Tetralogy of Fallot (TOF), using a northern European discovery set of 835 cases and 5159 controls. A region on chromosome 12q24 was associated (P = 1.4 × 10(-7)) and replicated convincingly (P = 3.9 × 10(-5)) in 798 cases and 2931 controls [per allele odds ratio (OR) = 1.27 in replication cohort, P = 7.7 × 10(-11) in combined populations]. Single nucleotide polymorphisms in the glypican 5 gene on chromosome 13q32 were also associated (P = 1.7 × 10(-7)) and replicated convincingly (P = 1.2 × 10(-5)) in 789 cases and 2927 controls (per allele OR = 1.31 in replication cohort, P = 3.03 × 10(-11) in combined populations). Four additional regions on chromosomes 10, 15 and 16 showed suggestive association accompanied by nominal replication. This study, the first genome-wide association study of a congenital heart malformation phenotype, provides evidence that common genetic variation influences the risk of TOF.

Contribution of global rare copy-number variants to the risk of sporadic congenital heart disease.

Previous studies have shown that copy-number variants (CNVs) contribute to the risk of complex developmental phenotypes. However, the contribution of global CNV burden to the risk of sporadic congenital heart disease (CHD) remains incompletely defined. We generated genome-wide CNV data by using Illumina 660W-Quad SNP arrays in 2,256 individuals with CHD, 283 trio CHD-affected families, and 1,538 controls. We found association of rare genic deletions with CHD risk (odds ratio [OR] = 1.8, p = 0.0008). Rare deletions in study participants with CHD had higher gene content (p = 0.001) with higher haploinsufficiency scores (p = 0.03) than they did in controls, and they were enriched with Wnt-signaling genes (p = 1 × 10(-5)). Recurrent 15q11.2 deletions were associated with CHD risk (OR = 8.2, p = 0.02). Rare de novo CNVs were observed in ~5% of CHD trios; 10 out of 11 occurred on the paternally transmitted chromosome (p = 0.01). Some of the rare de novo CNVs spanned genes known to be involved in heart development (e.g., HAND2 and GJA5). Rare genic deletions contribute ~4% of the population-attributable risk of sporadic CHD. Second to previously described CNVs at 1q21.1, deletions at 15q11.2 and those implicating Wnt signaling are the most significant contributors to the risk of sporadic CHD. Rare de novo CNVs identified in CHD trios exhibit paternal origin bias.

A common variant in the PTPN11 gene contributes to the risk of tetralogy of Fallot.

BACKGROUND: Tetralogy of Fallot (TOF) is the commonest cyanotic form of congenital heart disease. In 80% of cases, TOF behaves as a complex genetic condition exhibiting significant heritability. As yet, no common genetic variants influencing TOF risk have been robustly identified. METHODS AND RESULTS: Two hundred and seven haplotype-tagging single nucleotide polymorphisms in 22 candidate genes were genotyped in a test cohort comprising 362 nonsyndromic British white patients with TOF together with 717 unaffected parents of patients and 183 unrelated healthy controls. Single nucleotide polymorphisms with suggestive evidence of association in the test cohort (P<0.01) were taken forward for genotyping in an independent replication cohort comprising 392 cases of TOF, 218 unaffected parents of patients, and 1319 controls. Significant association was observed for 1 single nucleotide polymorphism, rs11066320 in the PTPN11 gene, in both the test and the replication cohort. Genotype at rs11066320 was associated with a per-allele odds ratio of 1.34 (95% confidence interval [CI], 1.19 to 1.52; P=2.9 × 10(-6)) in the total cohort of TOF cases and controls; this remained highly significant after Bonferroni correction for 207 analyses (corrected P=0.00061). Genotype at rs11066320 was responsible for a population-attributable risk of TOF of approximately 10%. CONCLUSIONS: Common variation in the linkage disequilibrium block including the PTPN11 gene contributes to the risk of nonsyndromic TOF. Rare mutations in PTPN11 are known to cause the autosomal dominant condition Noonan syndrome, which includes congenital heart disease, by upregulating Ras/mitogen-activated protein kinase (MAPK) signaling. Our results suggest a role for milder perturbations in PTPN11 function in sporadic, nonsyndromic congenital heart disease.

A new phenotype of brain iron accumulation with dystonia, optic atrophy, and peripheral neuropathy.

BACKGROUND: Neurodegeneration with brain iron accumulation is clinically and genetically heterogeneous because of mutations in at least 7 nuclear genes. METHODS: We performed homozygosity mapping and whole-exome sequencing in 2 brothers with brain iron accumulation from a consanguineous family. RESULTS: We identified a homozygous missense mutation in both brothers in the very recently identified chromosome 19 open-reading frame 12 gene. The disease presented before age 10 with slowly progressive tremor, dystonia, and spasticity. Additional features were optic atrophy, peripheral neuropathy, and learning difficulties. A raised serum creatine kinase indicated neuromuscular involvement, and compensatory mitochondrial proliferation implicated mitochondrial dysfunction as a pathological mechanism. CONCLUSIONS: Further studies are needed to explore the function of the chromosome 19 open-reading frame 12 gene, and extended genetic analysis on larger patient cohorts will provide more information about the presentation and frequency of this disease.

Phenotype-specific effect of chromosome 1q21.1 rearrangements and GJA5 duplications in 2436 congenital heart disease patients and 6760 controls.

Recurrent rearrangements of chromosome 1q21.1 that occur via non-allelic homologous recombination have been associated with variable phenotypes exhibiting incomplete penetrance, including congenital heart disease (CHD). However, the gene or genes within the ~1 Mb critical region responsible for each of the associated phenotypes remains unknown. We examined the 1q21.1 locus in 948 patients with tetralogy of Fallot (TOF), 1488 patients with other forms of CHD and 6760 ethnically matched controls using single nucleotide polymorphism genotyping arrays (Illumina 660W and Affymetrix 6.0) and multiplex ligation-dependent probe amplification. We found that duplication of 1q21.1 was more common in cases of TOF than in controls [odds ratio (OR) 30.9, 95% confidence interval (CI) 8.9-107.6); P = 2.2 × 10(-7)], but deletion was not. In contrast, deletion of 1q21.1 was more common in cases of non-TOF CHD than in controls [OR 5.5 (95% CI 1.4-22.0); P = 0.04] while duplication was not. We also detected rare (n = 3) 100-200 kb duplications within the critical region of 1q21.1 in cases of TOF. These small duplications encompassed a single gene in common, GJA5, and were enriched in cases of TOF in comparison to controls [OR = 10.7 (95% CI 1.8-64.3), P = 0.01]. These findings show that duplication and deletion at chromosome 1q21.1 exhibit a degree of phenotypic specificity in CHD, and implicate GJA5 as the gene responsible for the CHD phenotypes observed with copy number imbalances at this locus.

Genotype at the P554L variant of the hexose-6 phosphate dehydrogenase gene is associated with carotid intima-medial thickness.

OBJECTIVE: The combined thickness of the intima and media of the carotid artery (carotid intima-medial thickness, CIMT) is associated with cardiovascular disease and stroke. Previous studies indicate that carotid intima-medial thickness is a significantly heritable phenotype, but the responsible genes are largely unknown. Hexose-6 phosphate dehydrogenase (H6PDH) is a microsomal enzyme whose activity regulates corticosteroid metabolism in the liver and adipose tissue; variability in measures of corticosteroid metabolism within the normal range have been associated with risk factors for cardiovascular disease. We performed a genetic association study in 854 members of 224 families to assess the relationship between polymorphisms in the gene coding for hexose-6 phosphate dehydrogenase (H6PD) and carotid intima-medial thickness. METHODS: Families were ascertained via a hypertensive proband. CIMT was measured using B-mode ultrasound. Single nucleotide polymorphisms (SNPs) tagging common variation in the H6PD gene were genotyped. Association was assessed following adjustment for significant covariates including "classical" cardiovascular risk factors. Functional studies to determine the effect of particular SNPs on H6PDH were performed. RESULTS: There was evidence of association between the single nucleotide polymorphism rs17368528 in exon five of the H6PD gene, which encodes an amino-acid change from proline to leucine in the H6PDH protein, and mean carotid intima-medial thickness (p = 0.00065). Genotype was associated with a 5% (or 0.04 mm) higher mean carotid intima-medial thickness measurement per allele, and determined 2% of the population variability in the phenotype. CONCLUSIONS: Our results suggest a novel role for the H6PD gene in atherosclerosis susceptibility.

Exome sequencing identifies GATA-2 mutation as the cause of dendritic cell, monocyte, B and NK lymphoid deficiency.

The human syndrome of dendritic cell, monocyte, B and natural killer lymphoid deficiency presents as a sporadic or autosomal dominant trait causing susceptibility to mycobacterial and other infections, predisposition to myelodysplasia and leukemia, and, in some cases, pulmonary alveolar proteinosis. Seeking a genetic cause, we sequenced the exomes of 4 unrelated persons, 3 with sporadic disease, looking for novel, heterozygous, and probably deleterious variants. A number of genes harbored novel variants in person, but only one gene, GATA2, was mutated in all 4 persons. Each person harbored a different mutation, but all were predicted to be highly deleterious and to cause loss or mutation of the C-terminal zinc finger domain. Because GATA2 is the only common mutated gene in 4 unrelated persons, it is highly probable to be the cause of dendritic cell, monocyte, B, and natural killer lymphoid deficiency. This disorder therefore constitutes a new genetic form of heritable immunodeficiency and leukemic transformation.

Common variation at the 11-ß hydroxysteroid dehydrogenase type 1 gene is associated with left ventricular mass.

BACKGROUND: Polymorphisms in 11-ß hydroxysteroid dehydrogenase type 1 (11ß-HSD1, encoded by HSD11B1) have been reported to be associated with obesity-related cardiovascular risk factors, such as type II diabetes and hypertension. Left ventricular hypertrophy (LVH) is an independent risk factor for cardiovascular death associated with these factors but has significant additional heritability, the cause of which is undetermined. The 11ß-HSD1 is believed to maintain tonic inhibition of the mineralocorticoid receptor in cardiomyocytes, and mineralocorticoid receptor activation is involved in the pathophysiology of LVH. We assessed the association between polymorphisms in the HSD11B1 gene and left ventricular mass (LVM) in 248 families ascertained through a proband with hypertension. METHODS AND RESULTS: LVM was measured by electrocardiography and echocardiography in 868 and 829 participants, respectively. Single-nucleotide polymorphisms (SNPs) tagging common variation in the HSD11B1 gene were genotyped by mass spectrometry. The rs846910 SNP, which lies in the flanking region 5' to exon 1B of HSD11B1, was associated with LVM both by electrocardiography (≈5% lower LVM per copy of the rare allele, P=0.02) and by echocardiography (≈10% lower LVM per copy of the rare allele, P=0.003). Genotype explained 1% to 2% of the population variability in LVM, or approximately 5% of the heritable fraction. There were no significant associations between any HSD11B1 SNP and blood pressure or body mass index that could have confounded the association with LVM. CONCLUSIONS: Genotype at HSD11B1 has a small, but significant effect on LVM, apparently independently of any effect on obesity-related traits. These findings suggest a novel action of 11ß-HSD1 in the human cardiomyocyte, which may be of therapeutic importance.