Inferring the conservative causal core of gene regulatory networks

Background
Inferring gene regulatory networks from large-scale expression data is an important problem that received much attention in recent years. These networks have the potential to gain insights into causal molecular interactions of biological processes. Hence, from a methodological point of view, reliable estimation methods based on observational data are needed to approach this problem practically.

Results
In this paper, we introduce a novel gene regulatory network inference (GRNI) algorithm, called C3NET. We compare C3NET with four well known methods, ARACNE, CLR, MRNET and RN, conducting in-depth numerical ensemble simulations and demonstrate also for biological expression data from E. coli that C3NET performs consistently better than the best known GRNI methods in the literature. In addition, it has also a low computational complexity. Since C3NET is based on estimates of mutual information values in conjunction with a maximization step, our numerical investigations demonstrate that our inference algorithm exploits causal structural information in the data efficiently.

Conclusions
For systems biology to succeed in the long run, it is of crucial importance to establish methods that extract large-scale gene networks from high-throughput data that reflect the underlying causal interactions among genes or gene products. Our method can contribute to this endeavor by demonstrating that an inference algorithm with a neat design permits not only a more intuitive and possibly biological interpretation of its working mechanism but can also result in superior results.

Harnessing naturally randomized transcription to infer regulatory relationships among genes

We develop an approach utilizing randomized genotypes to rigorously infer causal regulatory relationships among genes at the transcriptional level, based on experiments in which genotyping and expression profiling are performed. This approach can be used to build transcriptional regulatory networks and to identify putative regulators of genes. We apply the method to an experiment in yeast, in which genes known to be in the same processes and functions are recovered in the resulting transcriptional regulatory network.

Resequencing of the CCL5 and CCR5 genes and investigation of variants for association with diabetic nephropathy

Chemokine (C-C motif) ligand 5 (CCL5) and chemokine (C-C motif) receptor 5 are implicated in the pathogenesis of diabetic nephropathy (DN). We hypothesize that variants in these genes may be associated with DN. The CCL5 and chemokine receptor type 5 (CCR5) genes were resequenced, variants identified (n=58), allele frequencies determined in 46 individuals (92 chromosomes) and efficient haplotype tag single-nucleotide polymorphisms (htSNPs) selected to effectively evaluate the common variation in these genes. One reportedly functional gene variant and eight htSNPs were genotyped in a case-control association study involving Caucasian individuals with type 1 diabetes (267 cases with DN and 442 non-nephropathic diabetic controls). Genotyping was performed using MassARRAY iPLEX, TaqMan, gel electrophoresis and direct capillary sequencing. After correction for multiple testing, there were no statistically significant associations between variants in the CCL5 and CCR5 genes and DN. Journal of Human Genetics (2010) 55, 248-251; doi:10.1038/jhg.2010.15; published online 5 March 2010

Multiple superoxide dismutase 1/splicing factor serine alanine 15 variants are associated with the development and progression of diabetic nephropathy:the Diabetes Control and Complications Trial/Epidemiology of Diabetes Interventions and Complications Genetics study

Despite familial clustering of nephropathy and retinopathy severity in type 1 diabetes, few gene variants have been consistently associated with these outcomes.

Association of Genetic Variants at 3q22 with Nephropathy in Patients with Type 1 Diabetes Mellitus

Diabetic nephropathy (DN) is the primary cause of morbidity and mortality in patients with type 1 diabetes mellitus (T1DM) and affects about 30% of these patients. We have previously localized a DN locus on chromosome 3q with suggestive linkage in Finnish individuals. Linkage to this region has also been reported earlier by several other groups. To fine map this locus, we conducted a multistage case-control association study in T1DM patients, comprising 1822 cases with nephropathy and 1874 T1DM patients free of nephropathy, from Finland, Iceland, and the British Isles. At the screening stage, we genotyped 3072 tag SNPs, spanning a 28 Mb region, in 234 patients and 215 controls from Finland. SNPs that met the significance threshold of p

Hypoglycaemia, liver necrosis and perinatal death in mice lacking all isoforms of phosphoinositide 3-kinase p85 alpha

Phosphoinositide 3-kinases produce 3'-phosphorylated phosphoinositides that act as second messengers to recruit other signalling proteins to the membrane(1). Pi3ks are activated by many extracellular stimuli and have been implicated in a variety of cellular responses(1). The Pi3k gene family is complex and the physiological roles of different classes and isoforms are not clear. The gene Pik3r1 encodes three proteins (p85 alpha, p55 alpha and p50 alpha) that serve as regulatory subunits of class I-A Pi3ks (ref. 2). Mice lacking only the p85a isoform are viable but display hypoglycaemia and increased insulin sensitivity correlating with upregulation of the p55 alpha and p50 alpha variants(3). Here we report that loss of all protein products of Pik3r1 results in perinatal lethality. We observed, among other abnormalities, extensive hepatocyte necrosis and chylous ascites, We also noted enlarged skeletal muscle fibres, brown fat necrosis and calcification of cardiac tissue. In liver and muscle, loss of the major regulatory isoform caused a great decrease in expression and activity of class I-A Pi3k catalytic subunits: nevertheless, homozygous mice still displayed hypoglycaemia, lower insulin levels and increased glucose tolerance. Our findings reveal that p55 alpha and/or p50 alpha are required for survival, but not for development of hypoglycaemia, in mice lacking p85 alpha.

Regulation of cyclin D1 by the BRCA1-BARD1 complex

Background BRCA1 and cyclin D₁ are both essential for normal breast development and mutation or aberration of their expression is associated with breast cancer [1,2]. Cyclin D₁ is best known as a G₁ cyclin where it regulates the G₁ to S phase transition by acting as a rate-limiting subunit of CDK4/6 kinase activity. More recently, however, Stacey has demonstrated that cyclin D₁ levels in G₂/M determine whether a cell continues to proliferate or exits the cell cycle [3]. The majority of BRCA1 in the cell is bound to BARD1 through their N-terminal RING domains. Heterodimerization is essential for the stability and correct localization of the complex and confers ubiquitin ligase activity to BRCA1. The importance of the ligase activity of BRCA1 to breast cancer development is inferred from the fact that N-terminal diseaseassociated mutations are proposed to reduce ligase activity [4]. Methods Protein–protein interactions were demonstrated using yeast-two-hybrid and coimmunoprecipitation. Protein levels were altered through overexpression, siRNA and antisense technology. The effect of proteasome inhibitors and cycloheximide treatment was also examined. Results We initially identified cyclin D₁ as a binding partner of BARD1 in a yeast-two-hybrid screen and defined the minimal binding region as the N-terminus of BARD1. This interaction was confirmed in vivo by coimmunoprecipitation. The N-terminus of BARD1 also binds BRCA1 and imparts ubiquitin ligase activity to the complex. Covalent modification of proteins with ubiquitin is a common regulatory mechanism in eukaryotic cells. Traditionally polyubiquitin chains linked through lysine 48 target proteins for degradation by the 26 S proteasome. We have demonstrated that cyclin D₁ protein levels are inversely related to BRCA1 and BARD1 levels in several model systems. Furthermore, regulation of cyclin D1 levels occurs through a post-transcriptional mechanism and requires the ligase activity of BRCA1. Interestingly, this phenomenon is cell-cycle regulated, occurring in G₂/M. Conclusion We propose that cyclin D₁ is a potential substrate for BRCA1 ubiquitination and that this targets cyclin D₁ for proteasomal-mediated degradation. Future work will focus on ascertaining the functional consequence of cyclin D₁ regulation by the BRCA1–BARD1 complex; in particular, the impact of BRCA1, mediated through regulation of cyclin D₁, on the proliferation versus differentiation decision.

Common variants at ABCA7, MS4A6A/MS4A4E, EPHA1, CD33 and CD2AP are associated with Alzheimer's disease

We sought to identify new susceptibility loci for Alzheimer's disease through a staged association study (GERAD+) and by testing suggestive loci reported by the Alzheimer's Disease Genetic Consortium (ADGC) in a companion paper. We undertook a combined analysis of four genome-wide association datasets (stage 1) and identified ten newly associated variants with P = 1 × 10(-5). We tested these variants for association in an independent sample (stage 2). Three SNPs at two loci replicated and showed evidence for association in a further sample (stage 3). Meta-analyses of all data provided compelling evidence that ABCA7 (rs3764650, meta P = 4.5 × 10(-17); including ADGC data, meta P = 5.0 × 10(-21)) and the MS4A gene cluster (rs610932, meta P = 1.8 × 10(-14); including ADGC data, meta P = 1.2 × 10(-16)) are new Alzheimer's disease susceptibility loci. We also found independent evidence for association for three loci reported by the ADGC, which, when combined, showed genome-wide significance: CD2AP (GERAD+, P = 8.0 × 10(-4); including ADGC data, meta P = 8.6 × 10(-9)), CD33 (GERAD+, P = 2.2 × 10(-4); including ADGC data, meta P = 1.6 × 10(-9)) and EPHA1 (GERAD+, P = 3.4 × 10(-4); including ADGC data, meta P = 6.0 × 10(-10)).

Genome-wide association of early-onset myocardial infarction with single nucleotide polymorphisms and copy number variants

We conducted a genome-wide association study testing single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) for association with early-onset myocardial infarction in 2,967 cases and 3,075 controls. We carried out replication in an independent sample with an effective sample size of up to 19,492. SNPs at nine loci reached genome-wide significance: three are newly identified (21q22 near MRPS6-SLC5A3-KCNE2, 6p24 in PHACTR1 and 2q33 in WDR12) and six replicated prior observations1-4 (9p21, 1p13 near CELSR2-PSRC1-SORT1, 10q11 near CXCL12, 1q41 in MIA3, 19p13 near LDLR and 1p32 near PCSK9). We tested 554 common copy number polymorphisms (>1% allele frequency) and none met the pre-specified threshold for replication (P < 10-3). We identified 8,065 rare CNVs but did not detect a greater CNV burden in cases compared to controls, in genes compared to the genome as a whole, or at any individual locus. SNPs at nine loci were reproducibly associated with myocardial infarction, but tests of common and rare CNVs failed to identify additional associations with myocardial infarction risk.

Dioxygenase-catalysed cis-dihydroxylation of meta-substituted phenols to yield cyclohexenone cis-diol and derived enantiopure cis-triol metabolites

cis-Dihydroxylation of meta-substituted phenol (m-phenol) substrates, to yield the corresponding cyclohexenone cis-diol metabolites, was catalysed by arene dioxygenases present in mutant and recombinant bacterial strains. The presence of cyclohexenone cis-diol metabolites and several of their cyclohexene and cyclohexane cis-triol derivatives was detected by LC-TOFMS analysis and confirmed by NMR spectroscopy. Structural and stereochemical analyses of chiral ketodiol bioproducts, was carried out using NMR and CD spectroscopy and stereochemical correlation methods. The formation of enantiopure cyclohexenone cis-diol metabolites is discussed in the context of postulated binding interactions of the m-phenol substrates at the active site of toluene dioxygenase (TDO).