Decay of linkage disequilibrium within genes across HGDP-CEPH human samples: most population isolates do not show increased LD

Background: It is well known that the pattern of linkage disequilibrium varies between human populations, with remarkable geographical stratification. Indirect association studies routinely exploit linkage disequilibrium around genes, particularly in isolated populations where it is assumed to be higher. Here, we explore both the amount and the decay of linkage disequilibrium with physical distance along 211 gene regions, most of them related to complex diseases, across 39 HGDP-CEPH population samples, focusing particularly on the populations defined as isolates. Within each gene region and population we use r2 between all possible single nucleotide polymorphism (SNP) pairs as a measure of linkage disequilibrium and focus on the proportion of SNP pairs with r2 greater than 0.8.Results: Although the average r2 was found to be significantly different both between and within continental regions, a much higher proportion of r2 variance could be attributed to differences between continental regions (2.8% vs. 0.5%, respectively). Similarly, while the proportion of SNP pairs with r2 > 0.8 was significantly different across continents for all distance classes, it was generally much more homogenous within continents, except in the case of Africa and the Americas. The only isolated populations with consistently higher LD in all distance classes with respect to their continent are the Kalash (Central South Asia) and the Surui (America). Moreover, isolated populations showed only slightly higher proportions of SNP pairs with r2 > 0.8 per gene region than non-isolated populations in the same continent. Thus, the number of SNPs in isolated populations that need to be genotyped may be only slightly less than in non-isolates. Conclusion: The "isolated population" label by itself does not guarantee a greater genotyping efficiency in association studies, and properties other than increased linkage disequilibrium may make these populations interesting in genetic epidemiology.

Patterns of evolutionary constraints on genes in humans

Background: Different regions in a genome evolve at different rates depending on structural and functional constraints. Some genomic regions are highly conserved during metazoan evolution, while other regions may evolve rapidly, either in all species or in a lineage-specific manner. A strong or even moderate change in constraints in functional regions, for example in coding regions, can have significant evolutionary consequences. Results: Here we discuss a novel framework, 'BaseDiver', to classify groups of genes in humans based on the patterns of evolutionary constraints on polymorphic positions in their coding regions. Comparing the nucleotide-level divergence among mammals with the extent of deviation from the ancestral base in the human lineage, we identify patterns of evolutionary pressure on nonsynonymous base-positions in groups of genes belonging to the same functional category. Focussing on groups of genes in functional categories, we find that transcription factors contain a significant excess of nonsynonymous base-positions that are conserved in other mammals but changed in human, while immunity related genes harbour mutations at base-positions that evolve rapidly in all mammals including humans due to strong preference for advantageous alleles. Genes involved in olfaction also evolve rapidly in all mammals, and in humans this appears to be due to weak negative selection. Conclusion: While recent studies have identified genes under positive selection in humans, our approach identifies evolutionary constraints on Gene Ontology groups identifying changes in humans relative to some of the other mammals.

Comparative analysis of cancer genes in the human and chimpanzee genomes

Background: Cancer is a major medical problem in modern societies. However, the incidence of this disease in non-human primates is very low. To study whether genetic differences between human and chimpanzee could contribute to their distinct cancer susceptibility, we have examined in the chimpanzee genome the orthologous genes of a set of 333 human cancer genes. Results: This analysis has revealed that all examined human cancer genes are present in chimpanzee, contain intact open reading frames and show a high degree of conservation between both species. However, detailed analysis of this set of genes has shown some differences in genes of special relevance for human cancer. Thus, the chimpanzee gene encoding p53 contains a Pro residue at codon 72, while this codon is polymorphic in humans and can code for Arg or Pro, generating isoforms with different ability to induce apoptosis or interact with p73. Moreover, sequencing of the BRCA1 gene has shown an 8 Kb deletion in the chimpanzee sequence that prematurely truncates the co-regulated NBR2 gene. Conclusion: These data suggest that small differences in cancer genes, as those found in tumor suppressor genes, might influence the differences in cancer susceptibility between human and chimpanzee. Nevertheless, further analysis will be required to determine the exact contribution of the genetic changes identified in this study to the different cancer incidence in non-human primates.

Structural and functional properties of genes involved in human cancer

Background: One of the main goals of cancer genetics is to identify the causative elements at the molecular level leading to cancer.Results: We have conducted an analysis of a set of genes known to be involved in cancer in order to unveil their unique features that can assist towards the identification of new candidate cancer genes. Conclusion: We have detected key patterns in this group of genes in terms of the molecular function or the biological process in which they are involved as well as sequence properties. Based on these features we have developed an accurate Bayesian classification model with which human genes have been scored for their likelihood of involvement in cancer.

Predicting cancer involvement of genes from heterogeneous data

Background: Systematic approaches for identifying proteins involved in different types of cancer are needed. Experimental techniques such as microarrays are being used to characterize cancer, but validating their results can be a laborious task. Computational approaches are used to prioritize between genes putatively involved in cancer, usually based on further analyzing experimental data. Results: We implemented a systematic method using the PIANA software that predicts cancer involvement of genes by integrating heterogeneous datasets. Specifically, we produced lists of genes likely to be involved in cancer by relying on: (i) protein-protein interactions; (ii) differential expression data; and (iii) structural and functional properties of cancer genes. The integrative approach that combines multiple sources of data obtained positive predictive values ranging from 23% (on a list of 811 genes) to 73% (on a list of 22 genes), outperforming the use of any of the data sources alone. We analyze a list of 20 cancer gene predictions, finding that most of them have been recently linked to cancer in literature. Conclusion: Our approach to identifying and prioritizing candidate cancer genes can be used to produce lists of genes likely to be involved in cancer. Our results suggest that differential expression studies yielding high numbers of candidate cancer genes can be filtered using protein interaction networks.

OSIRISv1.2: a named entity recognition system for sequence variants of genes in biomedical literature

Background: Single Nucleotide Polymorphisms, among other type of sequence variants, constitute key elements in genetic epidemiology and pharmacogenomics. While sequence data about genetic variation is found at databases such as dbSNP, clues about the functional and phenotypic consequences of the variations are generally found in biomedical literature. The identification of the relevant documents and the extraction of the information from them are hampered by the large size of literature databases and the lack of widely accepted standard notation for biomedical entities. Thus, automatic systems for the identification of citations of allelic variants of genes in biomedical texts are required. Results: Our group has previously reported the development of OSIRIS, a system aimed at the retrieval of literature about allelic variants of genes http://ibi.imim.es/osirisform.html. Here we describe the development of a new version of OSIRIS (OSIRISv1.2, http://ibi.imim.es/OSIRISv1.2.html webcite) which incorporates a new entity recognition module and is built on top of a local mirror of the MEDLINE collection and HgenetInfoDB: a database that collects data on human gene sequence variations. The new entity recognition module is based on a pattern-based search algorithm for the identification of variation terms in the texts and their mapping to dbSNP identifiers. The performance of OSIRISv1.2 was evaluated on a manually annotated corpus, resulting in 99% precision, 82% recall, and an F-score of 0.89. As an example, the application of the system for collecting literature citations for the allelic variants of genes related to the diseases intracranial aneurysm and breast cancer is presented. Conclusion: OSIRISv1.2 can be used to link literature references to dbSNP database entries with high accuracy, and therefore is suitable for collecting current knowledge on gene sequence variations and supporting the functional annotation of variation databases. The application of OSIRISv1.2 in combination with controlled vocabularies like MeSH provides a way to identify associations of biomedical interest, such as those that relate SNPs with diseases.

Tumor necrosis factor and transforming growth factor β regulate clock genes by controlling the expression of the cold inducible RNA-binding protein (CIRBP).

The circadian clock drives the rhythmic expression of a broad array of genes that orchestrate metabolism, sleep wake behavior, and the immune response. Clock genes are transcriptional regulators engaged in the generation of circadian rhythms. The cold inducible RNA-binding protein (CIRBP) guarantees high amplitude expression of clock. The cytokines TNF and TGFβ impair the expression of clock genes, namely the period genes and the proline- and acidic amino acid-rich basic leucine zipper (PAR-bZip) clock-controlled genes. Here, we show that TNF and TGFβ impair the expression of Cirbp in fibroblasts and neuronal cells. IL-1β, IL-6, IFNα, and IFNγ do not exert such effects. Depletion of Cirbp is found to increase the susceptibility of cells to the TNF-mediated inhibition of high amplitude expression of clock genes and modulates the TNF-induced cytokine response. Our findings reveal a new mechanism of cytokine-regulated expression of clock genes.

Genetik der essentiellen arteriellen Hypertonie [Essential hypertension - which genes can be held responsible for it?].

Hypertension is a common heritable cardiovascular risk factor. Some rare monogenic forms of hypertension have been described, but the majority of patients suffer from essential hypertension, for whom the underlying genetic mechanisms are not clear. Essential hypertension is a complex trait, involving multiple genes and environmental factors. Recently, progress in the identification of common genetic variants associated with essential hypertension has been made due to large-scale international collaborative projects. In this article we review the new research methods used as well as selected recent findings in this field.

A novel family of dehydrin-like proteins is involved in stress response in the human fungal pathogen Aspergillus fumigatus.

 During a search for genes controlling conidial dormancy in Aspergillus fumigatus, two dehydrin-like genes, DprA and DprB, were identified. The deduced proteins had repeated stretches of 23 amino acids that contained a conserved dehydrin-like protein (DPR) motif. Disrupted DprAΔ mutants were hypersensitive to oxidative stress and to phagocytic killing, whereas DprBΔ mutants were impaired in osmotic and pH stress responses. However, no effect was observed on their pathogenicity in our experimental models of invasive aspergillosis. Molecular dissection of the signaling pathways acting upstream showed that expression of DprA was dependent on the stress-activated kinase SakA and the cyclic AMP-protein kinase A (cAMP-PKA) pathways, which activate the bZIP transcription factor AtfA, while expression of DprB was dependent on the SakA mitogen-activated protein kinase (MAPK) pathway, and the zinc finger transcription factor PacC. Fluorescent protein fusions showed that both proteins were associated with peroxisomes and the cytosol. Accordingly, DprA and DprB were important for peroxisome function. Our findings reveal a novel family of stress-protective proteins in A. fumigatus and, potentially, in filamentous ascomycetes.

Repeated exposure to Lutzomyia intermedia sand fly saliva induces local expression of interferon-inducible genes both at the site of injection in mice and in human blood.

During a blood meal, Lutzomyia intermedia sand flies transmit Leishmania braziliensis, a parasite causing tegumentary leishmaniasis. In experimental leishmaniasis, pre-exposure to saliva of most blood-feeding sand flies results in parasite establishment in absence of any skin damages in mice challenged with dermotropic Leishmania species together with saliva. In contrast, pre-immunization with Lu. intermedia salivary gland sonicate (SGS) results in enhanced skin inflammatory exacerbation upon co-inoculation of Lu. intermedia SGS and L. braziliensis. These data highlight potential unique features of both L. braziliensis and Lu. intermedia. In this study, we investigated the genes modulated by Lu. intermedia SGS immunization to understand their potential impact on the subsequent cutaneous immune response following inoculation of both SGS and L. braziliensis. The cellular recruitment and global gene expression profile was analyzed in mice repeatedly inoculated or not with Lu. intermedia. Microarray gene analysis revealed the upregulation of a distinct set of IFN-inducible genes, an immune signature not seen to the same extent in control animals. Of note this INF-inducible gene set was not induced in SGS pre-immunized mice subsequently co-inoculated with SGS and L. braziliensis. These data suggest the parasite prevented the upregulation of this Lu. intermedia saliva-related immune signature. The presence of these IFN-inducible genes was further analyzed in peripheral blood mononuclear cells (PBMCs) sampled from uninfected human individuals living in a L. braziliensis-endemic region of Brazil thus regularly exposed to Lu. intermedia bites. PBMCs were cultured in presence or absence of Lu. intermedia SGS. Using qRT-PCR we established that the IFN-inducible genes induced in the skin of SGS pre-immunized mice, were also upregulated by SGS in PBMCs from human individuals regularly exposed to Lu. intermedia bites, but not in PBMCs of control subjects. These data demonstrate that repeated exposure to Lu. intermedia SGS induces the expression of potentially host-protective IFN-inducible genes.

Domain shuffling in a sensor protein contributed to the evolution of insect pathogenicity in plant-beneficial Pseudomonas protegens.

Pseudomonas protegens is a biocontrol rhizobacterium with a plant-beneficial and an insect pathogenic lifestyle, but it is not understood how the organism switches between the two states. Here, we focus on understanding the function and possible evolution of a molecular sensor that enables P. protegens to detect the insect environment and produce a potent insecticidal toxin specifically during insect infection but not on roots. By using quantitative single cell microscopy and mutant analysis, we provide evidence that the sensor histidine kinase FitF is a key regulator of insecticidal toxin production. Our experimental data and bioinformatic analyses indicate that FitF shares a sensing domain with DctB, a histidine kinase regulating carbon uptake in Proteobacteria. This suggested that FitF has acquired its specificity through domain shuffling from a common ancestor. We constructed a chimeric DctB-FitF protein and showed that it is indeed functional in regulating toxin expression in P. protegens. The shuffling event and subsequent adaptive modifications of the recruited sensor domain were critical for the microorganism to express its potent insect toxin in the observed host-specific manner. Inhibition of the FitF sensor during root colonization could explain the mechanism by which P. protegens differentiates between the plant and insect host. Our study establishes FitF of P. protegens as a prime model for molecular evolution of sensor proteins and bacterial pathogenicity.

Interactions of Adeno-associated virus Rep78 protein with the host cell

Abstract : Adeno-associated virus (AAV) is a small DNA virus belonging to the familiy of Parvoviridae. Its genome contains two genes : the rep gene encoding four non structural proteins (Rep78, 68, 52 and 40) implicated in transcription, replication and site-specific integration of the viral DNA and the cap gene encoding three capsid proteins. AAV does not cause any disease, but is studied in view of its potential use to treat several diseases. An interesting property of AAV is its antiproliferative effect. Two elements of AAV can inhibit cell growth. Firstly, the single stranded viral DNA is recognized in cells as damaged DNA leading to either a G2 block or cell death depending on p53 status. Secondly, the two larger Rep proteins (Rep78 and 68) also arrest the cell cycle when they are expressed at high levels. Rep78 in particular induces a complete cell cycle arrest in all the phases, including S phase. Such a strong S phase arrest is rarely seen in other conditions. It was thus interesting to determine how Rep78 could induce it. We found that this strong block is the consequence of Rep78's effects on at least two pathways. Rep78 induces a DNA damage response by producing nicks in the cellular chromatin. Furthermore, Rep78 can bind to the cellular phosphatase Cdc25A and prevent its binding to its substrates CDK2 and CDK1, thus inhibiting its activity. A mutational analysis of Rep78 protein determined that its endonuclease activity is responsible for the DNA damage response and its zinc finger domain for Cdc25A inhibition. The combined expression of two mutants each defective for one of these activities, or these two activities obtained independently of Rep78, could restore the complete cell cycle block, indicating that these two effects of Rep78 are likely to explain completely the cell cycle block it induces. Secondly, the lack of pathogenicity of AAV, its broad range of infection and its ability to integrate site-specifically in human chromosome 19 make it an interesting potential vector for gene therapy. However site-specific integration is only possible in the presence of Rep78/68 whose gene is removed in recombinant AAV vectors. In this part of the study, we tried to introduce Rep protein separately from recombinant AAV vectors to promote their site-specific integration. For that purpose, a fusion protein, TAT-Rep, comprising Rep78/68 joined to the human immunodeficiency virus Tat protein was produced. It had the ability to enter cells and remain active there for a short period. Its activity was sufficient to mediate transcription from the p5 promoter, second-strand synthesis of a recombinant AAV and probably site-specific integration. Résumé : Le virus associé à l'adénovirus (AAV) est un petit virus à ADN qui fait partie de la famille des Parvoviridae. Son génome contient deux gènes : le gène rep code pour quatre protéines (Rep78, 68, 52 et 40) qui participent à la transcription, la réplication et l'intégration du virus et le gène cap code pour les trois protéines de capside. AAV ne produit pas de maladie, mais pourrait au contraire être utilisé pour en soigner. Sa bénignité, sa capacité à infecter différents types de cellules et son intégration spécifique en font un vecteur potentiel pour la thérapie génique. Pour qu'il puisse s'intégrer spécifiquement, il a besoin de la protéine Rep78 ou 68, mais ce gène doit être enlevé des vecteurs pour la thérapie génique. Le but de la première partie de cette étude était d'introduire Rep78 ou 68 dans des cellules en même temps qu'un AAV recombinant, mais indépendamment afin de permettre une intégration spécifique. La stratégie utilisée était de produire une protéine de fusion (TAT-Rep) qui peut entrer dans des cellules si elle est présente dans leur milieu. Cette protéine entrait bien dans les cellules et y était active favorisant ainsi l'intégration spécifique. Une deuxième propriété d'AAV, son effet anti-prolifératif, est intéressante dans le cadre de certaines maladies comme le cancer. Deux éléments d'AAV en sont responsables. D'abord, son ADN simple brin active une réponse cellulaire à l'ADN endommagé et arrête les cellules en G2 ou provoque leur mort. De plus, la protéine Rep78 d'AAV peut fortement bloquer le cycle cellulaire à toutes les phases, même en phase S, ce qui est rare. C'est pourquoi nous avons essayé de comprendre cet effet. Nous avons remarqué que Rep78 doit agir sur deux fronts pour obtenir ce fort bloc. D'un côté, Rep78 introduit des coupures simple brin sur l'ADN de la cellule ce qui active une réponse cellulaire à l'ADN endommagé qui passe par ATM. D'un autre côté, Rep78 lie une phosphatase cellulaire, Cdc25A, et l'empêche ainsi de lier ses substrats CDK2 et CDK1 et donc d'être active. Finalement, à l'aide de mutants de Rep78, nous avons déterminé que l'activité endonuclease de Rep78 était nécessaire pour induire une réponse cellulaire via ATM et que le domaine C-terminal appelé «zinc finger » était responsable de la liaison avec Cdc25A. En co-exprimant deux mutants, qui n'ont chacun qu'un des effets de Rep78, ou en obtenant les deux effets de Rep78 indépendamment d'elle, nous avons obtenu un bloc complet du cycle cellulaire similaire à celui obtenu avec Rep78. Il est donc probable que ces deux effets de Rep78 sont suffisants pour expliquer comment elle arrive à arrêter le cycle cellulaire si efficacement.

Association of candidate genes with phenotypic traits relevant to anorexia nervosa.

This analysis is a follow-up to an earlier investigation of 182 genes selected as likely candidate genetic variations conferring susceptibility to anorexia nervosa (AN). As those initial case-control results revealed no statistically significant differences in single nucleotide polymorphisms, herein, we investigate alternative phenotypes associated with AN. In 1762 females, using regression analyses, we examined the following: (i) lowest illness-related attained body mass index; (ii) age at menarche; (iii) drive for thinness; (iv) body dissatisfaction; (v) trait anxiety; (vi) concern over mistakes; and (vii) the anticipatory worry and pessimism versus uninhibited optimism subscale of the harm avoidance scale. After controlling for multiple comparisons, no statistically significant results emerged. Although results must be viewed in the context of limitations of statistical power, the approach illustrates a means of potentially identifying genetic variants conferring susceptibility to AN because less complex phenotypes associated with AN are more proximal to the genotype and may be influenced by fewer genes. Copyright © 2011 John Wiley & Sons, Ltd and Eating Disorders Association.

Identifying protein-coding genes in genomic sequences.

The vast majority of the biology of a newly sequenced genome is inferred from the set of encoded proteins. Predicting this set is therefore invariably the first step after the completion of the genome DNA sequence. Here we review the main computational pipelines used to generate the human reference protein-coding gene sets.

Ferripyochelin uptake genes are involved in pyochelin-mediated signalling in Pseudomonas aeruginosa.

In response to iron starvation, Pseudomonas aeruginosa produces the siderophore pyochelin. When secreted to the extracellular environment, pyochelin chelates iron and transports it to the bacterial cytoplasm via its specific outer-membrane receptor FptA and the inner-membrane permease FptX. Exogenously added pyochelin also acts as a signal which induces the expression of the pyochelin biosynthesis and uptake genes by activating PchR, a cytoplasmic regulatory protein of the AraC/XylS family. The importance of ferripyochelin uptake genes in this regulation was evaluated. The fptA and fptX genes were shown to be part of the fptABCX ferripyochelin transport operon, which is conserved in Burkholderia sp. and Rhodospirillum rubrum. The fptB and fptC genes were found to be dispensable for utilization of pyochelin as an iron source, for signalling and for pyochelin production. By contrast, mutations in fptA and fptX not only interfered with pyochelin utilization, but also affected signalling and diminished siderophore production. It is concluded from this that pyochelin-mediated signalling operates to a large extent via the ferripyochelin transport system.