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.

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).

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.

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.

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.

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.

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.

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.

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).

Effects of rhamnolipids from Pseudomonas aeruginosa DS10-129 on luminescent bacteria: toxicity and modulation of cadmium bioavailability.

In this study, the mixture of mono- and di-rhamnolipids produced by Pseudomonas aeruginosa DS10-129 was characterized for its toxicity and modulatory effects on Cd availability to different bacteria. Gram-negative naturally bioluminescent Vibrio fischeri and recombinant bioluminescent Pseudomonas fluorescens, P. aeruginosa, Escherichia coli, and Gram-positive Bacillus subtilis were used as model organisms. Rhamnolipids reduced the bioluminescence of these bacteria in less than a second of exposure even in relatively low concentrations (30-min EC(50) 45-167 mg l(-1)). Toxicity of Cd to Gram-negative bacteria (30-min EC(50) values 0.16 mg l(-1) for E. coli, 0.96 mg l(-1) for P. fluorescens, and 4.4 mg l(-1) for V. fischeri) was remarkably (up to 10-fold) reduced in the presence of 50 mg l(-1) rhamnolipids. Interestingly, the toxicity of Cd to Gram-positive B. subtilis (30-min EC(50) value 0.49 mg l(-1)) was not affected by rhamnolipids. Rhamnolipids had an effect on desorption of Cd from soil: 40 mg l(-1) rhamnolipids increased the water-extracted fraction of Cd twice compared with untreated control. However, this additionally desorbed fraction of Cd remained bound with rhamnolipids and was not available to bacteria. Hence, in carefully chosen concentrations (still effectively complexing heavy metals but not yet toxic to soil bacteria), rhamnolipids could be applied in remediation of polluted areas.

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.

The C-532T polymorphism of the angiotensinogen gene is associated with pulse pressure: a possible explanation for heterogeneity in genetic association studies of AGT and hypertension.

BACKGROUND: Many previous studies have investigated whether there is an association between genotypes at the angiotensinogen (AGT) gene and hypertensive status, but few have incorporated quantitative data. Although meta-analyses support a possible effect of AGT variants on blood pressure (BP), substantial unexplained between-study heterogeneity has been observed. We hypothesized that a primary effect of AGT variants on arterial stiffness (and thus pulse pressure) might explain such heterogeneity, and tested for such an effect in a family study. METHODS: We studied 1425 individuals from 248 families ascertained through a proband with essential hypertension. BP was measured using 24 h ambulatory monitoring, and polymorphisms of the AGT gene that had been previously associated with hypertension and/or plasma angiotensinogen levels were typed. Pulse pressure was used as a measurement of arterial stiffness. RESULTS: We observed a highly significant association between genotypes at the AGT C-532T polymorphism and pulse pressure (p = 0.00006). Each T allele was associated with a 5% lower pulse pressure (that is, an additive effect). This resulted from opposing genotypic effects to (slightly) lower systolic BP and (slightly) elevate DBP. CONCLUSIONS: These results suggest that genetic variation at the angiotensinogen locus may primarily affect arterial stiffness, and therefore pulse pressure. The heterogeneity between previous genetic studies of AGT and hypertension status could in part be explained by this finding, since case selection criteria based on systolic BP, diastolic BP, or both would result in different levels of selection for the -532T allele.

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.

INSIG-2 promoter polymorphism and obesity related phenotypes: association study in 1428 members of 248 families.

BACKGROUND: Obesity is a major public health problem. Body mass index (BMI) is a highly heritable phenotype but robust associations of genetic polymorphisms to BMI or other obesity-related phenotypes have been difficult to establish. Recently a large genetic association study showed evidence for association of the single nucleotide polymorphism (SNP) rs7566605, which lies 10 Kb 5' to the first exon of the insulin-induced gene 2 (INSIG-2), with obesity in several cohorts. We tested this polymorphism for association with body mass related phenotypes in a large family study whose mean BMI was consistent with moderate overweight. METHODS: We studied 1428 members of 248 British Caucasian families who had been ascertained through a proband with hypertension. We measured BMI, waist and hip circumference, and plasma levels of leptin. We genotyped the rs7566605 SNP using a restriction fragment length polymorphism assay, and carried out a family-based association test for quantitative traits related to obesity using the statistical programs MERLIN and QTDT. RESULTS: We observed no significant association between genotype at rs7566605 and covariate-adjusted (for age, sex, alcohol consumption, smoking and exercise habit) log-transformed BMI, waist measurement, hip measurement, waist-to-hip ratio, or plasma levels of leptin. CONCLUSION: There was no association between genotype at rs7566605 and obesity-related phenotypes in this British Caucasian population. These families were in general moderately overweight, few members being severely obese. Our result indicates that this polymorphism has little if any effect on BMI within the normal to moderately overweight range. The effects of this polymorphism on body mass may be restricted to those already predisposed to at least moderate obesity as a result of environmental factors and other predisposing genotypes.

The degradation of n-hexadecane in soil by thermophilic geobacilli.

Soil microcosms have been used to demonstrate the ability of indigenous soil thermophiles to degrade effectively a representative alkane (n-hexadecane). A fragment of the alkane mono-oxygenase gene (alkB) was amplified from thermophilic Geobacillus thermoleovorans strain T70 by PCR using degenerate primers. The amplicon demonstrated 96% sequence similarity with the alkB gene from Rhodococcus erythropolis. Critical controls ensured that the positive PCR signal detected was not a result of mesophilic soil organisms. A reverse transcription PCR (RT-PCR) assay was developed to determine if expression of the gene was inducible in the presence of an alkane or constitutively expressed in soil. In the presence of n-hexadecane, expression of the alkane mono-oxygenase gene was induced in pure cultures and soil samples and was dependent on temperature. No positive RT-PCR signal was detected at mesophilic growth temperatures either in pure cultures or in soil microcosms, whereas at 55 degrees C positive RT-PCR signals were obtained for both pure cultures of T70 and soil samples. Many different amplicons of the alkB gene fragment were obtained from the soil used in the microcosms. Thirty cloned fragments yielded 27 different sequences showing 85-96% sequence similarity with the alkB sequence of T70. To establish that the amplified alkB gene sequences from soil were derived from thermophilic geobacilli, additional strains were isolated on a selective medium containing n-hexadecane as sole carbon source. The 16S rRNA gene sequences were determined to identify the 50 isolates obtained (G. thermoleovorans, 27; G. caldoxylosilyticus, 17; G. pallidus, 2; G. toebiii, 1; Geobacillus sp., 3) representing 18 different strains and alkB gene sequences determined and deposited with the European Bioinformatics Institute.

Rhamnolipid biosurfactant production by strains of Pseudomonas aeruginosa using low-cost raw materials.

This study was aimed at the development of economical methods for higher yields of biosurfactant by suggesting the use of low-cost raw materials. Two oil-degrading strains, Pseudomonas aeruginosa GS9-119 and DS10-129, were used to optimize a substrate for maximum rhamnolipid production. Among the two strains, the latter produced maxima of 4.31, 2.98, and 1.77 g/L rhamnolipid biosurfactant using soybean oil, safflower oil, and glycerol, respectively. The yield of biosurfactant steadily increased even after the bacterial cultures reached the stationary phase of growth. Characterization of rhamnolipids using mass spectrometry revealed the presence of dirhamnolipids (Rha-Rha-C(10)-C(10)). Emulsification activity of the rhamnolipid biosurfactant produced by P. aeruginosa DS10-129 was greater than 70% using all the hydrocarbons tested, including xylene, benzene, hexane, crude oil, kerosene, gasoline, and diesel. P. aeruginosa GS9-119 emulsified only hexane and kerosene to that level.

Enhanced bioremediation of n-alkane in petroleum sludge using bacterial consortium amended with rhamnolipid and micronutrients.

The purpose of the present study was to investigate possible methods to enhance the rate of biodegradation of oil sludge from crude oil tank bottom, thus reducing the time usually required for bioremediation. Enhancement of biodegradation was achieved through bioaugmentation and biostimulation. About 10% and 20% sludge contaminated sterile and non-sterile soil samples were treated with bacterial consortium (BC), rhamnolipid biosurfactant (RL) and nitrogen, phosphorus and potassium (NPK) solution. Maximum n-alkane degradation occurred in the 10% sludge contaminated soil samples. The effects of treatment carried out with the non-sterile soil samples were more pronounced than in the sterile soils. Maximum degradation was achieved after the 56th day of treatment. n-Alkanes in the range of nC8-nC11 were degraded completely followed by nC12-nC21, nC22-nC31 and nC32-nC40 with percentage degradations of 100%, 83-98%, 80-85% and 57-73% respectively. Statistical analysis using analysis of variance and Duncan's multiple range test revealed that the level of amendments, incubation time and combination of amendments significantly influenced bacterial growth, protein concentration and surface tension at a 1% probability level. All tested additives BC, NPK and RL had significant positive effects on the bioremediation of n-alkane in petroleum sludge.

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.