Interactions of the eNOS and ACE genes and cigarette smoking in chronic obstructive pulmonary disease.

Background: Chronic obstructive pulmonary disease (COPD) is a complex disorder with unexplained heritability. Interactions of genetic and environmental factors are thought to be crucial in COPD. So, we aim to examine interactions of the endothelial nitric oxide synthase (eNOS) and angiotensin converting enzyme (ACE) genes and cigarette smoking in COPD. Methods: The eNOS G 894T and ACE ID variants were analyzed in 122 COPD patients and 200 controls from Serbia. The effect of the variants on COPD was assessed by logistic regression. Interactions between eNOS, ACE and cigarette smoking in COPD were evaluated using a case-control model. Interaction between the genes was analyzed in silico. Results: No effect of the eNOS G 894T and ACE ID variants on COPD was found in our study. Gene-gene interaction between the eN OS T T and A CE D was identified (p=0.033) in COPD. The interaction is realized within the complex network of biochemical pathways. Gene-environment interactions between the eNOS T and cigarette smoking (p=0.013), and the ACE II and cigarette smoking (p=0.009) were detected in COPD in our study. Conclusions: This is the first research to reveal interactions of the eNOS and ACE genes and cigarette smoking in COPD progressing our understanding of COPD heritability and contributing to the development of appropriate treatments.

Transcription factor site dependencies in human, mouse and rat genomes.

BACKGROUND: It is known that transcription factors frequently act together to regulate gene expression in eukaryotes. In this paper we describe a computational analysis of transcription factor site dependencies in human, mouse and rat genomes. RESULTS: Our approach for quantifying tendencies of transcription factor binding sites to co-occur is based on a binding site scoring function which incorporates dependencies between positions, the use of information about the structural class of each transcription factor (major/minor groove binder), and also considered the possible implications of varying GC content of the sequences. Significant tendencies (dependencies) have been detected by non-parametric statistical methodology (permutation tests). Evaluation of obtained results has been performed in several ways: reports from literature (many of the significant dependencies between transcription factors have previously been confirmed experimentally); dependencies between transcription factors are not biased due to similarities in their DNA-binding sites; the number of dependent transcription factors that belong to the same functional and structural class is significantly higher than would be expected by chance; supporting evidence from GO clustering of targeting genes. Based on dependencies between two transcription factor binding sites (second-order dependencies), it is possible to construct higher-order dependencies (networks). Moreover results about transcription factor binding sites dependencies can be used for prediction of groups of dependent transcription factors on a given promoter sequence. Our results, as well as a scanning tool for predicting groups of dependent transcription factors binding sites are available on the Internet. CONCLUSION: We show that the computational analysis of transcription factor site dependencies is a valuable complement to experimental approaches for discovering transcription regulatory interactions and networks. Scanning promoter sequences with dependent groups of transcription factor binding sites improve the quality of transcription factor predictions.

Matrix metalloproteinases gene variants in idiopathic disseminated bronchiectasis.

BACKGROUND: The excess of matrix metalloproteinases (MMPs) might be associated with the airways destruction or dilatation in bronchiectasis. The functional promoter polymorphisms of MMP1 and MMP9 genes, involved in the extracellular matrix remodeling, might increase the expression of MMPs leading to the development of bronchiectasis. METHODS: Detection of MMP1 G-1607GG and MMP9 C-1562T gene variants was performed on 37 patients with idiopathic disseminated bronchiectasis and 102 control subjects. We also described a novel method for simple and rapid detection of MMP1 G-1607GG polymorphism. RESULTS: The frequency of -1607GG allele was significantly higher in the group of patients than in control subjects (P = 0.014). The heterozygote genotype showed association with bronchiectasis (odds ratio, 5.3; 95% confidence intervals, 1.4-20.0). The association was even stronger in homozygotes for -1607GG allele (odds ration, 8.7; 95% confidence intervals, 1.9-41.0). The allelic and genotype frequencies of MMP9 C-1562T variant did not show significant differences between the groups. CONCLUSIONS: This is the first report concerning a role of MMP1 G-1607GG and MMP9 C-1562T variants in pathogenesis of idiopathic disseminated bronchiectasis. The results of our study revealed the association of -1607GG allele and the lack of association of MMP9 C-1562T variant with the disease.

Quality estimation of multiple sequence alignments by Bayesian hypothesis testing.

UNLABELLED: In this work we present a web-based tool for estimating multiple alignment quality using Bayesian hypothesis testing. The proposed method is very simple, easily implemented and not time consuming with a linear complexity. We evaluated method against a series of different alignments (a set of random and biologically derived alignments) and compared the results with tools based on classical statistical methods (such as sFFT and csFFT). Taking correlation coefficient as an objective criterion of the true quality, we found that Bayesian hypothesis testing performed better on average than the classical methods we tested. This approach may be used independently or as a component of any tool in computational biology which is based on the statistical estimation of alignment quality. AVAILABILITY: http://www.fmi.ch/groups/functional.genomics/tool.htm. SUPPLEMENTARY INFORMATION: Supplementary data are available from http://www.fmi.ch/groups/functional.genomics/tool-Supp.htm.

Position dependencies in transcription factor binding sites.

MOTIVATION: Most of the available tools for transcription factor binding site prediction are based on methods which assume no sequence dependence between the binding site base positions. Our primary objective was to investigate the statistical basis for either a claim of dependence or independence, to determine whether such a claim is generally true, and to use the resulting data to develop improved scoring functions for binding-site prediction. RESULTS: Using three statistical tests, we analyzed the number of binding sites showing dependent positions. We analyzed transcription factor-DNA crystal structures for evidence of position dependence. Our final conclusions were that some factors show evidence of dependencies whereas others do not. We observed that the conformational energy (Z-score) of the transcription factor-DNA complexes was lower (better) for sequences that showed dependency than for those that did not (P < 0.02). We suggest that where evidence exists for dependencies, these should be modeled to improve binding-site predictions. However, when no significant dependency is found, this correction should be omitted. This may be done by converting any existing scoring function which assumes independence into a form which includes a dependency correction. We present an example of such an algorithm and its implementation as a web tool. AVAILABILITY: http://promoterplot.fmi.ch/cgi-bin/dep.html

The CFTR M470V gene variant as a potential modifier of COPD severity: study of Serbian population.

Chronic obstructive pulmonary disease (COPD) is a complex disease influenced by genetic and environmental factors. Cystic fibrosis transmembrane conductance regulator (CFTR) protein is an important component of the lung tissue homeostasis, involved in the regulation of the rate of mucociliary clearance. As it is known that certain CFTR variants have consequences on the function of CFTR protein, the aim of this study was to examine the possible role of F508del, M470V, Tn locus, and R75Q variants in COPD development and modulation. Total number of 86 COPD patients and 102 control subjects were included in the study. Possible association between COPD susceptibility, severity, and onset of the disease and allele or genotype of four analyzed CFTR variants was examined. No associations were detected between COPD development, onset of the disease and tested CFTR alleles and genotypes. However, VV470 genotype was associated with mild/moderate COPD stages in comparison to severe/very severe ones (OR = 0.29, 95%CI = 0.11-0.80, p = 0.016). Our study showed that patients with VV470 genotype had a 3.4-fold decreased risk for the appearance of severe/very severe COPD symptoms, and the obtained results indicate that this genotype may have a protective role. These results also suggest the importance of studying CFTR gene as a modifier of this disease.

n-gram-based classification and unsupervised hierarchical clustering of genome sequences.

In this paper we address the problem of automated classification of isolates, i.e., the problem of determining the family of genomes to which a given genome belongs. Additionally, we address the problem of automated unsupervised hierarchical clustering of isolates according only to their statistical substring properties. For both of these problems we present novel algorithms based on nucleotide n-grams, with no required preprocessing steps such as sequence alignment. Results obtained experimentally are very positive and suggest that the proposed techniques can be successfully used in a variety of related problems. The reported experiments demonstrate better performance than some of the state-of-the-art methods. We report on a new distance measure between n-gram profiles, which shows superior performance compared to many other measures, including commonly used Euclidean distance.

Gene-environment interaction between the MMP9 C-1562T promoter variant and cigarette smoke in the pathogenesis of chronic obstructive pulmonary disease.

The aetiology of chronic obstructive pulmonary disease (COPD) is complex. While cigarette smoking is a well-established cause of COPD, a myriad of assessed genetic factors has given conflicting data. Since gene-environment interactions are thought to be implicated in aetiopathogenesis of COPD, we aimed to examine the matrix metalloproteinase (MMP) 9 C-1562T (rs3918242) functional variant and cigarette smoke in the pathogenesis of this disease. The distribution of the MMP9 C-1562T variant was analyzed in COPD patients and controls with normal pulmonary function from Serbia. Interaction between the C-1562T genetic variant and cigarette smoking was assessed using a case-control model. The response of the C-1562T promoter variant to cigarette smoke condensate (CSC) exposure was examined using a dual luciferase reporter assay. The frequency of T allele carriers was higher in the COPD group than in smoker controls (38.4% vs. 20%; OR = 2.7, P = 0.027). Interaction between the T allele and cigarette smoking was identified in COPD occurrence (OR = 4.38, P = 0.005) and severity (P = 0.001). A functional analysis of the C-1562T variant demonstrated a dose-dependent and allele-specific response (P < 0.01) to CSC. Significantly higher MMP9 promoter activity following CSC exposure was found for the promoter harboring the T allele compared to the promoter harboring the C allele (P < 0.05). Our study is the first to reveal an interaction between the MMP9-1562T allele and cigarette smoke in COPD, emphasising gene-environment interactions as a possible cause of lung damage in the pathogenesis of COPD. Environ. Mol. Mutagen. 57:447-454, 2016. © 2016 Wiley Periodicals, Inc.

SARS-CoV genome polymorphism: a bioinformatics study.

A dataset of 103 SARS-CoV isolates (101 human patients and 2 palm civets) was investigated on different aspects of genome polymorphism and isolate classification. The number and the distribution of single nucleotide variations (SNVs) and insertions and deletions, with respect to a "profile", were determined and discussed ("profile" being a sequence containing the most represented letter per position). Distribution of substitution categories per codon positions, as well as synonymous and non-synonymous substitutions in coding regions of annotated isolates, was determined, along with amino acid (a.a.) property changes. Similar analysis was performed for the spike (S) protein in all the isolates (55 of them being predicted for the first time). The ratio Ka/Ks confirmed that the S gene was subjected to the Darwinian selection during virus transmission from animals to humans. Isolates from the dataset were classified according to genome polymorphism and genotypes. Genome polymorphism yields to two groups, one with a small number of SNVs and another with a large number of SNVs, with up to four subgroups with respect to insertions and deletions. We identified three basic nine-locus genotypes: TTTT/TTCGG, CGCC/TTCAT, and TGCC/TTCGT, with four subgenotypes. Both classifications proposed are in accordance with the new insights into possible epidemiological spread, both in space and time.

Association of Functional Variants of Phase I and II Genes with Chronic Obstructive Pulmonary Disease in a Serbian Population.

BACKGROUND: Chronic obstructive pulmonary disease (COPD) is a complex disorder characterized by increased oxidative stress. Functional genetic variants of phase I and II genes are implicated in oxidants-antioxidants imbalance and may be involved in COPD development. In this study, we aimed to investigate the role of cytochrome P450 (CYP), glutathione S-transferase (GST) and microsomal epoxide hydrolase (mEH) functional variants in the pathogenesis of COPD in a Serbian population. METHODS: The genotypes of 122 COPD patients and 100 controls with normal lung function were determined for CYP1A1 *1A/*2A, CYP2E1 *1A/*5B, GSTM1 null, GSTT1 null GSTP1 Ile105Val, mEH Tyr113His and mEH His139Arg gene variants. RESULTS: Results obtained showed that GSTM1 null variant was significantly more represented in COPD patients than in controls (61.5% vs. 47.0%; OR=1.80; p=0.042). Also, a significant difference was observed for combinations of GSTM1 null and GSTP1 105Val/(Val) (38.5% vs. 24.0%; OR=1.98; p=0.029), as well as for CYP1A1 *1A/*2A, GSTM1 null and mEH 113His/(His) genotypes (7.4% vs. 1.0%; OR=7.88; p=0.025). CONCLUSIONS: These are the first data concerning the analysis of the variants of phase I and II genes in the pathogenesis of COPD in a Serbian population. Results obtained in this study open up the possibility for thorough analyses of the role of genetic factors in COPD on larger cohorts. Also, they implicate the importance of previously described genetic associations with COPD in our population, as well as reveal a new one, not reported so far.

Computational structural analysis: multiple proteins bound to DNA.

BACKGROUND: With increasing numbers of crystal structures of proteinratioDNA and proteinratioproteinratioDNA complexes publically available, it is now possible to extract sufficient structural, physical-chemical and thermodynamic parameters to make general observations and predictions about their interactions. In particular, the properties of macromolecular assemblies of multiple proteins bound to DNA have not previously been investigated in detail. METHODOLOGY/PRINCIPAL FINDINGS: We have performed computational structural analyses on macromolecular assemblies of multiple proteins bound to DNA using a variety of different computational tools: PISA; PROMOTIF; X3DNA; ReadOut; DDNA and DCOMPLEX. Additionally, we have developed and employed an algorithm for approximate collision detection and overlapping volume estimation of two macromolecules. An implementation of this algorithm is available at http://promoterplot.fmi.ch/Collision1/. The results obtained are compared with structural, physical-chemical and thermodynamic parameters from proteinratioprotein and single proteinratioDNA complexes. Many of interface properties of multiple proteinratioDNA complexes were found to be very similar to those observed in binary proteinratioDNA and proteinratioprotein complexes. However, the conformational change of the DNA upon protein binding is significantly higher when multiple proteins bind to it than is observed when single proteins bind. The water mediated contacts are less important (found in less quantity) between the interfaces of components in ternary (proteinratioproteinratioDNA) complexes than in those of binary complexes (proteinratioprotein and proteinratioDNA).The thermodynamic stability of ternary complexes is also higher than in the binary interactions. Greater specificity and affinity of multiple proteins binding to DNA in comparison with binary protein-DNA interactions were observed. However, protein-protein binding affinities are stronger in complexes without the presence of DNA. CONCLUSIONS/SIGNIFICANCE: Our results indicate that the interface properties: interface area; number of interface residues/atoms and hydrogen bonds; and the distribution of interface residues, hydrogen bonds, van der Walls contacts and secondary structure motifs are independent of whether or not a protein is in a binary or ternary complex with DNA. However, changes in the shape of the DNA reduce the off-rate of the proteins which greatly enhances the stability and specificity of ternary complexes compared to binary ones.