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.

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

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.

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.

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.

Common variation in the CD36 (fatty acid translocase) gene is associated with left-ventricular mass.

AIMS: Genetic variation in the fatty acid translocase (CD36) gene has been shown in animal models to affect several risk factors for the development of left-ventricular hypertrophy, but this phenotype has not, thus far, been investigated in humans. We examined the relationship between common genetic polymorphisms in the CD36 gene and left-ventricular mass. METHODS AND RESULTS: We studied a cohort of 255 families comprising 1425 individuals ascertained via a hypertensive proband. Seven single-nucleotide polymorphisms which together tagged common genetic variation in the CD36 gene were genotyped using a SEQUENOM MALDI-TOF instrument. There was evidence of association between the rs1761663 polymorphism in intron 1 of the CD36 gene and left-ventricular mass determined either by echocardiography (P=0.003, N=780) or electrocardiography (P=0.001, N=814). There was also association between rs1761663 genotype and body mass index (P<0.001, N=1354). Genotype was associated with between 2 and 8% differences in these phenotypes per allele. After adjustment for the effect of body mass index, there remained significant associations between genotype and left ventricular mass measured either by echo (P=0.017) or ECG (P=0.007). CONCLUSIONS: Genotype at the rs1761663 polymorphism has independent effects both on body mass index and left-ventricular mass. Genes with such pleiotropic effects may be particularly attractive therapeutic targets for interventions to modify multiple risk factors for cardiovascular events.

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.

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.

Ambulatory blood pressure is associated with polymorphic variation in P2X receptor genes.

The P2X receptor gene family encodes a series of proteins that function as ATP-gated nonselective ion channels. P2X receptor channels are involved in transducing aldosterone-mediated signaling in the distal renal tubule and are potential candidate genes for blood pressure regulation. We performed a family based quantitative genetic association study in 248 families ascertained through a proband with hypertension to investigate the relationship between common genetic variation in the P2X4, P2X6, and P2X7 genes and ambulatory blood pressure. We genotyped 28 single nucleotide polymorphisms, which together captured the common genetic variability in the 3 genes. We corrected our results for multiple comparisons specifying a false discovery rate of 5%. We found significant evidence of association between the single nucleotide polymorphism rs591874 in the first intron of the P2X7 gene and blood pressure. The strongest association was found for nighttime diastolic blood pressure (P=0.0032), although association was present for both systolic and diastolic blood pressures measured by an observer during the day and at night. Genotypes were associated with a 0.2 SD ( approximately 2.5 mm Hg) difference in night diastolic blood pressure per allele and accounted for approximately 1% of the total variability in this measurement. Other suggestive associations were found, but these were nonsignificant after correction for multiple testing. These genetic data suggest that drugs affecting P2X receptor signaling may have promise as clinical antihypertensive agents.

Geobacillus debilis sp. nov., a novel obligately thermophilic bacterium isolated from a cool soil environment, and reassignment of Bacillus pallidus to Geobacillus pallidus comb. nov.

Several aerobic, motile, rod-shaped, thermophilic, spore-forming Geobacillus bacteria predominantly giving a Gram-positive staining reaction were isolated from a cool soil environment in Northern Ireland and taxonomically investigated. Two isolates, F10 and Tf(T), showed low 16S rRNA gene sequence similarity to recognized members of the genus Geobacillus. Phylogenetic tree investigation using neighbour-joining, maximum-likelihood and parsimony methods indicated that strains F10 and Tf(T) represent a single novel species, for which the name Geobacillus debilis sp. nov. is proposed, with type strain Tf(T) (=DSM 16016(T)=NCIMB 13995(T)) and which belongs to a subgroup of the genus Geobacillus comprising Geobacillus toebii and Geobacillus caldoxylosilyticus. However, G. debilis showed closest affinities to Bacillus pallidus, which we propose should become Geobacillus pallidus comb. nov.

The frequency and characteristics of highly thermophilic bacteria in cool soil environments.

Following enrichment at 70 degrees C and 80 degrees C, five highly thermophilic aerobic eubacteria have been isolated from cool soil environments. These organisms show a temperature range for growth of 40-80 degrees C and have optimal and very high growth rates around 70 degrees C with generation times less than 30 min. All isolates are narrow rods, which stain Gram-negative, but have a Gram-positive cell wall structure and only one of five isolates is a spore former. All cultures contain a small proportion of previously unreported extremely long flexuous rods, which can be seen to divide eventually. Biochemical testing of five strains reveals a significant ability to utilize alkanes and some aromatic hydrocarbons. Using polymerase chain reaction (PCR) and restriction fragment length polymorphism (RFLP) of 16S rDNA the five strains were differentiated into three categories, which paralleled the biochemical results. 16S rDNA sequences showed high similarity with thermophilic Bacillus species now reclassified as Geobacillus. These bacteria are present in high numbers in apparently all soils and the question is raised of how these organisms, which are apparently unable to grow at the temperatures experienced in these cool soils, are so prominent.