Conditional mutual inclusive information enables accurate quantification of associations in gene regulatory networks

Mutual information (MI), a quantity describing the nonlinear dependence between two random variables, has been widely used to construct gene regulatory networks (GRNs). Despite its good performance, MI cannot separate the direct regulations from indirect ones among genes. Although the conditional mutual information (CMI) is able to identify the direct regulations, it generally underestimates the regulation strength, i.e. it may result in false negatives when inferring gene regulations. In this work, to overcome the problems, we propose a novel concept, namely conditional mutual inclusive information (CMI2), to describe the regulations between genes. Furthermore, with CMI2, we develop a new approach, namely CMI2NI (CMI2-based network inference), for reverse-engineering GRNs. In CMI2NI, CMI2 is used to quantify the mutual information between two genes given a third one through calculating the Kullback–Leibler divergence between the postulated distributions of including and excluding the edge between the two genes. The benchmark results on the GRNs from DREAM challenge as well as the SOS DNA repair network in Escherichia coli demonstrate the superior performance of CMI2NI. Specifically, even for gene expression data with small sample size, CMI2NI can not only infer the correct topology of the regulation networks but also accurately quantify the regulation strength between genes. As a case study, CMI2NI was also used to reconstruct cancer-specific GRNs using gene expression data from The Cancer Genome Atlas (TCGA). CMI2NI is freely accessible at http://www.comp-sysbio.org/cmi2ni.     

Concerted and differential actions of two enzymatic domains underlie Rad5 contributions to DNA damage tolerance

Many genome maintenance factors have multiple enzymatic activities. In most cases, how their distinct activities functionally relate with each other is unclear. Here we examined the conserved budding yeast Rad5 protein that has both ubiquitin ligase and DNA helicase activities. The Rad5 ubiquitin ligase activity mediates PCNA poly-ubiquitination and subsequently recombination-based DNA lesion tolerance. Interestingly, the ligase domain is embedded in a larger helicase domain comprising seven consensus motifs. How features of the helicase domain influence ligase function is controversial. To clarify this issue, we use genetic, 2D gel and biochemical analyses and show that a Rad5 helicase motif important for ATP binding is also required for PCNA poly-ubiquitination and recombination-based lesion tolerance. We determine that this requirement is due to a previously unrecognized contribution of the motif to the PCNA and ubiquitination enzyme interaction, and not due to its canonical role in supporting helicase activity. We further show that Rad5′s helicase-mediated contribution to replication stress survival is separable from recombination. These findings delineate how two Rad5 enzymatic domains concertedly influence PCNA modification, and unveil their discrete contributions to stress tolerance.     

Allele-specific copy number profiling by next-generation DNA sequencing

The progression and clonal development of tumors often involve amplifications and deletions of genomic DNA. Estimation of allele-specific copy number, which quantifies the number of copies of each allele at each variant loci rather than the total number of chromosome copies, is an important step in the characterization of tumor genomes and the inference of their clonal history. We describe a new method, falcon, for finding somatic allele-specific copy number changes by next generation sequencing of tumors with matched normals. falcon is based on a change-point model on a bivariate mixed Binomial process, which explicitly models the copy numbers of the two chromosome haplotypes and corrects for local allele-specific coverage biases. By using the Binomial distribution rather than a normal approximation, falcon more effectively pools evidence from sites with low coverage. A modified Bayesian information criterion is used to guide model selection for determining the number of copy number events. Falcon is evaluated on in silico spike-in data and applied to the analysis of a pre-malignant colon tumor sample and late-stage colorectal adenocarcinoma from the same individual. The allele-specific copy number estimates obtained by falcon allows us to draw detailed conclusions regarding the clonal history of the individual''s colon cancer.     

Molecular Mechanism of Nicotine Degradation by a Newly Isolated Strain,Ochrobactrum sp. Strain SJY1

A newly isolated strain, SJY1, identified as Ochrobactrum sp., utilizes nicotine as a sole source of carbon, nitrogen, and energy. Strain SJY1 could efficiently degrade nicotine via a variant of the pyridine and pyrrolidine pathways (the VPP pathway), which highlights bacterial metabolic diversity in relation to nicotine degradation. A 97-kbp DNA fragment containing six nicotine degradation-related genes was obtained by gap closing from the genome sequence of strain SJY1. Three genes, designated vppB, vppD, and vppE, in the VPP pathway were cloned and heterologously expressed, and the related proteins were characterized. The vppB gene encodes a flavin-containing amine oxidase converting 6-hydroxynicotine to 6-hydroxy-N-methylmyosmine. Although VppB specifically catalyzes the dehydrogenation of 6-hydroxynicotine rather than nicotine, it shares higher amino acid sequence identity with nicotine oxidase (38%) from the pyrrolidine pathway than with its isoenzyme (6-hydroxy-l-nicotine oxidase, 24%) from the pyridine pathway. The vppD gene encodes an NADH-dependent flavin-containing monooxygenase, which catalyzes the hydroxylation of 6-hydroxy-3-succinoylpyridine to 2,5-dihydroxypyridine. VppD shows 62% amino acid sequence identity with the hydroxylase (HspB) from Pseudomonas putida strain S16, whereas the specific activity of VppD is ∼10-fold higher than that of HspB. VppE is responsible for the transformation of 2,5-dihydroxypyridine. Sequence alignment and phylogenetic analysis suggested that the VPP pathway, which evolved independently from nicotinic acid degradation, might have a closer relationship with the pyrrolidine pathway. The proteins and functional pathway identified here provide a sound basis for future studies aimed at a better understanding of molecular principles of nicotine degradation.     

Laccaria bicolor aquaporin LbAQP1 is required for Hartig net development in trembling aspen (Populus tremuloides)

The development of ectomycorrhizal associations is crucial for growth of many forest trees. However, the signals that are exchanged between the fungus and the host plant during the colonization process are still poorly understood. In this study, we have identified the relationship between expression patterns of Laccaria bicolor aquaporin LbAQP1 and the development of ectomycorrhizal structures in trembling aspen (Populus tremuloides) seedlings. The peak expression of LbAQP1 was 700‐fold higher in the hyphae within the root than in the free‐living mycelium after 24 h of direct interaction with the roots. Moreover, in LbAQP1 knock‐down strains, a non‐mycorrhizal phenotype was developed without the Hartig net and the expression of the mycorrhizal effector protein MiSSP7 quickly declined after an initial peak on day 5 of interaction of the fungal hyphae with the roots. The increase in the expression of LbAQP1 required a direct contact of the fungus with the root and it modulated the expression of MiSSP7. We have also determined that LbAQP1 facilitated NO, H₂O₂ and CO₂ transport when heterologously expressed in yeast. The report demonstrates that the L. bicolor aquaporin LbAQP1 acts as a molecular signalling channel, which is fundamental for the development of Hartig net in root tips of P. tremuloides.