Family Size Evolution in Drosophila Chemosensory Gene Families: A Comparative Analysis with a Critical Appraisal of Methods

Gene turnover rates and the evolution of gene family sizes are important aspects of genome evolution. Here, we use curated sequence data of the major chemosensory gene families from Drosophila-the gustatory receptor, odorant receptor, ionotropic receptor, and odorant-binding protein families-to conduct a comparative analysis among families, exploring different methods to estimate gene birth and death rates, including an ad hoc simulation study. Remarkably, we found that the state-of-the-art methods may produce very different rate estimates, which may lead to disparate conclusions regarding the evolution of chemosensory gene family sizes in Drosophila. Among biological factors, we found that a peculiarity of D. sechellia's gene turnover rates was a major source of bias in global estimates, whereas gene conversion had negligible effects for the families analyzed herein. Turnover rates vary considerably among families, subfamilies, and ortholog groups although all analyzed families were quite dynamic in terms of gene turnover. Computer simulations showed that the methods that use ortholog group information appear to be the most accurate for the Drosophila chemosensory families. Most importantly, these results reveal the potential of rate heterogeneity among lineages to severely bias some turnover rate estimation methods and the need of further evaluating the performance of these methods in a more diverse sampling of gene families and phylogenetic contexts. Using branch-specific codon substitution models, we find further evidence of positive selection in recently duplicated genes, which attests to a nonneutral aspect of the gene birth-and-death process.

Rarity of DNA sequence alterations in the promoter region of the human androgen receptor gene

The human androgen receptor (AR) gene promoter lies in a GC-rich region containing two principal sites of transcription initiation and a putative Sp1 protein-binding site, without typical "TATA" and "CAAT" boxes. It has been suggested that mutations within the 5'untranslated region (5'UTR) may contribute to the development of prostate cancer by changing the rates of gene transcription and/or translation. In order to investigate this question, the aim of the present study was to search for the presence of mutations or polymorphisms at the AR-5'UTR in 92 prostate cancer patients, where histological diagnosis of adenocarcinoma was established in specimens obtained from transurethral resection or after prostatectomy. The AR-5'UTR was amplified by PCR from genomic DNA samples of the patients and of 100 healthy male blood donors, included as controls. Conformation-sensitive gel electrophoresis was used for DNA sequence alteration screening. Only one band shift was detected in one individual from the blood donor group. Sequencing revealed a new single nucleotide deletion (T) in the most conserved portion of the promoter region at position +36 downstream from the transcription initiation site I. Although the effect of this specific mutation remains unknown, its rarity reveals the high degree of sequence conservation of the human androgen promoter region. Moreover, the absence of detectable variation within the critical 5'UTR in prostate cancer patients indicates a low probability of its involvement in prostate cancer etiology.

New insights into small molecules inhibitors and protein-protein interactions of VirB8 : a critical conserved component of the type IV secretion system.

Les systèmes bactériens de sécrétion de type IV (T4SS) sont constitués d’un ensemble de 8 à 12 protéines conservées. Ces dernières sont utilisées lors de la translocation de protéines, la translocation de complexes ADN-protéines mais aussi pour le transport de ces derniers au travers de la membrane cellulaire. Les T4SS, en tant que facteurs de virulence pour beaucoup de pathogènes comme Brucella suis, sont donc d’excellents modèles cibles pour le développement de médicaments d’antivirulence. Ces médicaments, en privant le pathogène de son facteur essentiel de virulence : le T4SS, constituent une alternative ou encore une amélioration des traitements antibiotiques utilisés actuellement. VirB8, un facteur d’assemblage conservé dans le T4SS, forme des dimères qui sont importants pour la fonction des T4SS dans ces pathogènes. De par ses interactions multiples, VirB8 est un excellent modèle pour l’analyse des facteurs d’assemblage mais aussi en tant que cible de médicaments qui empêcheraient son interaction avec d’autres protéines et qui, in fine, désarmeraient les bactéries en les privant de leur fonctions essentielles de virulence. À ce jour, nous savons qu’il existe un équilibre monomère-dimère et un processus d’homodimerization de VirB8 dont l’importance est vitale pour la fonctionnement biologique des T4SSs. En se basant sur des essais quantitatifs d’interaction, nous avons identifié (i) des sites potentiels d’interaction avec d’autres protéines VirB du T4SS mais aussi (ii) isolé des petites molécules inhibitrices afin de tester la fonction protéique de VirB8. Afin de déterminer les acides aminés importants pour l’hétérodimérization de VirB8 avec VirB10, nous avons effectué des expériences de mutagenèse aléatoire, de phage display et d’arrimage moléculaire in silico. Ces expériences ont démontré l’importance de trois acides aminés localisés sur le feuillet β : R160, S162, T164 et I165. Ces derniers seraient importants pour l’association de VirB8 avec VirB10 étant donné que leur mutagenèse entraine une diminution de la formation du complexe VirB8-VirB10. L’objectif actuel de notre projet de recherche est de pouvoir mieux comprendre mais aussi d’évaluer le rôle de VirB8 dans l’assemblage du T4SS. Grace à un méthode de criblage adaptée à partir de la structure de VirB8, nous avons pu identifié une petite molécule inhibitrice BAR-068, qui aurait un rôle prometteur dans l’inhibition du T4SS. Nous avons utilisé la spectroscopie par fluorescence, l’essai à deux hybrides, le cross-linking et la cristallographie afin de déterminer le mécanisme d'interaction existant entre VirB8 et BAR-068. Ces travaux pourraient permettre de nombreuses avancées, notamment en termes de compréhension des mécanismes d’inhibition du T4SS. Notre objectif ultime est de pouvoir caractériser la séquence d’évènements essentiels à l’assemblage et au fonctionnement du T4SS. De manière globale, notre projet de recherche permettrait de révéler les grands principes d’assemblage des protéines membranaires, les processus de sécrétion de protéines chez les bactéries mais aussi de proposer une nouvelle stratégie lors du développement de drogues antimicrobiennes.

Sequence and phylogenetic analysis of viper venom serine proteases

Snakebites are a major neglected tropical disease responsible for as many as 95000 deaths every year worldwide. Viper venom serine proteases disrupt haemostasis of prey and victims by affecting various stages of the blood coagulation system. A better understanding of their sequence, structure, function and phylogenetic relationships will improve the knowledge on the pathological conditions and aid in the development of novel therapeutics for treating snakebites. A large dataset for all available viper venom serine proteases was developed and analysed to study various features of these enzymes. Despite the large number of venom serine protease sequences available, only a small proportion of these have been functionally characterised. Although, they share some of the common features such as a C-terminal extension, GWG motif and disulphide linkages, they vary widely between each other in features such as isoelectric points, potential N-glycosylation sites and functional characteristics. Some of the serine proteases contain substitutions for one or more of the critical residues in catalytic triad or primary specificity pockets. Phylogenetic analysis clustered all the sequences in three major groups. The sequences with substitutions in catalytic triad or specificity pocket clustered together in separate groups. Our study provides the most complete information on viper venom serine proteases to date and improves the current knowledge on the sequence, structure, function and phylogenetic relationships of these enzymes. This collective analysis of venom serine proteases will help in understanding the complexity of envenomation and potential therapeutic avenues.

Identification of critical residues in the active site of porcine membrane-bound aminopeptidase P

The membrane-bound form of mammalian aminopeptidase P (AP-P; EC 3.4. 11.9) is a mono-zinc-containing enzyme that lacks any of the typical metal binding motifs found in other zinc metalloproteases. To identify residues involved in metal binding and catalysis, sequence and structural information was used to align the sequence of porcine membrane-bound AP-P with other members of the peptidase clan MG, including Escherichia coli AP-P and methionyl aminopeptidases. Residues predicted to be critical for activity were mutated and the resultant proteins were expressed in COS-1 cells. Immunoelectrophoretic blot analysis was used to compare the levels of expression of the mutant proteins, and their ability to hydrolyze bradykinin and Gly-Pro-hydroxyPro was assessed. Asp449, Asp460, His523, Glu554, and Glu568 are predicted to serve as metal ion ligands in the active site, and mutagenesis of these residues resulted in fully glycosylated proteins that were catalytically inactive. Mutation of His429 and His532 also resulted in catalytically inactive proteins, and these residues, by analogy with E. coli AP-P, are likely to play a role in shuttling protons during catalysis. These studies indicate that mammalian membrane-bound AP-P has an active-site configuration similar to that of other members of the peptidase clan MG, which is compatible with either a dual metal ion model or a single metal ion in the active site. The latter model is consistent, however, with the known metal stoichiometry of both the membrane-bound and cytosolic forms of AP-P and with a recently proposed model for methionyl aminopeptidase.

Self-assembly of a Model Peptide incorporating a Hexa-Histidine sequence attached to an Oligo-Alanine sequence, and binding to Gold NTA/Nickel nanoparticles

Amyloid fibrils are formed by a model surfactant-like peptide (Ala)10-(His)6 containing a hexahistidine tag. This peptide undergoes a remarkable two-step self-assembly process with two distinct critical aggregation concentrations (cac’s), probed by fluorescence techniques. A micromolar range cac is ascribed to the formation of prefibrillar structures, whereas a millimolar range cac is associated with the formation of well-defined but more compact fibrils. We examine the labeling of these model tagged amyloid fibrils using Ni-NTA functionalized gold nanoparticles (Nanogold). Successful labeling is demonstrated via electron microscopy imaging. The specificity of tagging does not disrupt the β-sheet structure of the peptide fibrils. Binding of fibrils and Nanogold is found to influence the circular dichroism associated with the gold nanoparticle plasmon absorption band. These results highlight a new approach to the fabrication of functionalized amyloid fibrils and the creation of peptide/nanoparticle hybrid materials.

Critical role for an acidic amino acid region in platelet signaling by the HemITAM (hemi-immunoreceptor tyrosine-based activation motif) containing receptor CLEC-2 (C-type lectin receptor-2)

CLEC-2 is a member of new family of C-type lectin receptors characterized by a cytosolic YXXL downstream of three acidic amino acids in a sequence known as a hemITAM (hemi-immunoreceptor tyrosine-based activation motif). Dimerization of two phosphorylated CLEC-2 molecules leads to recruitment of the tyrosine kinase Syk via its tandem SH2 domains and initiation of a downstream signaling cascade. Using Syk-deficient and Zap-70-deficient cell lines we show that hemITAM signaling is restricted to Syk and that the upstream triacidic amino acid sequence is required for signaling. Using surface plasmon resonance and phosphorylation studies, we demonstrate that the triacidic amino acids are required for phosphorylation of the YXXL. These results further emphasize the distinct nature of the proximal events in signaling by hemITAM relative to ITAM receptors.

Mobius Sequence in a Girl and Arthrogryposis in her Half-Brother: Distinct Phenotypes Caused By Prenatal Injuries

Mobius sequence is a congenital facial and abducens nerve palsy, frequently associated to abnormalities of extremities. Arthrogryposis multiplex congenital is defined as a congenital fixation of multiple joints seldom of neurogenic origin. Both sequences must have a genetic origin, but usually are sporadic cases related to environmental factors such as drugs exposition and maternal trauma. A 5-year-old girl and a 1-year-old boy were born with Mobius sequence and arthrogryposis multiplex congenital, respectively. During pregnancies, the mother had vaginal bleeding at 7 weeks and used crack (free-based cocaine) in the first trimester, respectively. The girl also has equinovarus talipes and autistic behavior. The boy has arthrogryposis with flexion contractures of the feet and knees. A vascular disruption, due to hemorrhage and cocaine exposure, causing a transient ischemic insult to embryos in a critical period of development may be responsible for distinct phenotypes in these cases.

Soliton solutions for quasilinear Schrodinger equations with critical growth

In this paper we establish the existence of standing wave solutions for quasilinear Schrodinger equations involving critical growth. By using a change of variables, the quasilinear equations are reduced to semilinear one. whose associated functionals are well defined in the usual Sobolev space and satisfy the geometric conditions of the mountain pass theorem. Using this fact, we obtain a Cerami sequence converging weakly to a solution v. In the proof that v is nontrivial, the main tool is the concentration-compactness principle due to P.L. Lions together with some classical arguments used by H. Brezis and L. Nirenberg (1983) in [9]. (C) 2009 Elsevier Inc. All rights reserved.

Repair Welding Effects on the Bending Fatigue Strength of AISI 4130 Aeronautical Steel used in a Critical to the Flight-Safety Structure

The aim of this study was to analyze the effect of successive TIG (tungsten inert gas) welding repairs on the reverse bending fatigue strength of AISI 4130 steel, which is widely used in components critical to the flight-safety. In order to simulate the abrupt maneuvers, wind bursts, motor vibration and helixes efforts, which generate cyclic bending loadings at the welded joints of a specific aircraft component called motor cradle, experimental reverse bending fatigue tests were carried out on specimens made from hot-rolled steel plate, 1.10 mm (0.043 in) thick, by mean of a SCHENK PWS equipment, with load ratio R = -1, under constant amplitude, at 30 Hz frequency and room temperature. It was observed that the bending fatigue strength decreases after the TIG (Tungsten Inert Gas) welding process application on AISI 4130 steel, with subsequent decrease due to re-welding sequence as well. Microstructural analyses and microhardness measurements on the base material, heat-affected zone (HAZ) and weld metal, as well as the effects of the weld bead geometry on the obtained results, have complemented this study.

Considerations about the welding repair effects on the structural integrity of an airframe critical to the flight-safety

Structures critical to the flight-safety are commonly submitted to several maintenance repairs at the welded joints in order to prolong the in-service life of aircrafts. The aim of this study is to analyze the effects of Tungsten Inert Gas (TIG) welding repair on the structural integrity of the AISI 4130 aeronautical steel by means of experimental fatigue crack growth tests in base-material, heat-affected zone (HAZ) and weld metal. The tests were performed on hot-rolled steel plate specimens, 0.89 mm thick, with load ratio R = 0.1, constant amplitude, at 10 Hz frequency and room temperature. Increase of the fracture resistance was observed in the weld metal but decreasing in the HAZ after repair. The results were associated to microhardness and microstructural changes with the welding sequence. (C) 2010 Published by Elsevier Ltd.

Quantitative real-time PCR identifies a critical region of deletion on 22q13 related to prognosis in oral cancer

Quantitative real time PCR was performed on genomic DNA from 40 primary oral carcinomas and the normal adjacent tissues. The target genes ECGFB, DIA1, BIK, and PDGFB and the microsatellite markers D22S274 and D22S277, mapped on 22q13, were selected according to our previous loss of heterozygosity findings in head and neck tumors. Quantitative PCR relies on the comparison of the amount of product generated from a target gene and that generated from a disomic reference gene (GAPDH-housekeeping gene). Reactions have been performed with normal control in triplicates, using the 7700 Sequence Detection System (PE Applied Biosystems). Losses in the sequences D22S274 (22q13.31) and in the DIA1 (22q13.2-13.31) gene were detected in 10 out of 40 cases (25%) each. Statistically significant correlations were observed for patients with relative copy number loss of the marker D22S274 and stages T3-T4 of disease (P=0.025), family history of cancer (P = 0.001), and death (P = 0.021). Relative copy number loss involving the DIA1 gene was correlated to family history of cancer (P<0.001), death (P=0.002), and consumption of alcohol (P=0.026). Log-rank test revealed a significant decrease in survival (P=0.0018) for patients with DIA1 gene loss. Relative copy number losses detected in these sequences may be related to disease progression and a worse prognosis in patients with oral cancer.

Sequence, evolution and ligand binding properties of mammalian Duffy antigen/receptor for chemokines

The Duffy antigen/receptor for chemokine, DARC, acts as a widely expressed promiscuous chemokine receptor and as the erythrocyte receptor for Plasmodium vivax. To gain insight into the evolution and structure/function relations of DARC, we analyzed the binding of anti-human Fy monoclonal antibodies (mAbs) and human chemokines to red blood cells (RBCs) from 11 nonhuman primates and two nonprimate mammals, and we elucidated the structures of the DARC genes from gorilla, gibbon, baboon, marmoset, tamarin, night monkey and cattle. CXCL-8 and CCL-5 chemokine binding analysis indicated that the promiscuous binding profile characteristic of DARC is conserved across species. Among three mAbs that detected the Fy6 epitope by flow cytometric analysis of human and chimpanzee RBCs, only one reacted with night monkey and squirrel monkey. Only chimpanzee RBCs bound a significant amount of the anti-Fy3 mAb. Fy3 was also poorly detected on RBCs from gorilla, baboon and rhesus monkey, but not from new world monkeys. Alignment of DARC homologous sequences allowed us to construct a phylogenetic tree in which all branchings were in accordance with current knowledge of primate phylogeny. Although DARC was expected to be under strong internal and external selection pressure, in order to maintain chemokine binding and avoid Plasmodium vivax binding, respectively, our present study did not provide arguments in favor of a selection pressure on the extracellular domains involved in ligand specificity. The amino acid variability of DARC-like polypeptides was found to be well correlated with the hydrophylicity indexes, with the highest divergence on the amino-terminal extracellular domain. Analysis of the deduced amino acid sequences highlighted the conservation of some amino acid residues, which should prove to be critical for the structural and functional properties of DARC.

A clustering method for robust and reliable large scale functional and structural protein sequence annotation

Bioinformatics, in the last few decades, has played a fundamental role to give sense to the huge amount of data produced. Obtained the complete sequence of a genome, the major problem of knowing as much as possible of its coding regions, is crucial. Protein sequence annotation is challenging and, due to the size of the problem, only computational approaches can provide a feasible solution. As it has been recently pointed out by the Critical Assessment of Function Annotations (CAFA), most accurate methods are those based on the transfer-by-homology approach and the most incisive contribution is given by cross-genome comparisons. In the present thesis it is described a non-hierarchical sequence clustering method for protein automatic large-scale annotation, called “The Bologna Annotation Resource Plus” (BAR+). The method is based on an all-against-all alignment of more than 13 millions protein sequences characterized by a very stringent metric. BAR+ can safely transfer functional features (Gene Ontology and Pfam terms) inside clusters by means of a statistical validation, even in the case of multi-domain proteins. Within BAR+ clusters it is also possible to transfer the three dimensional structure (when a template is available). This is possible by the way of cluster-specific HMM profiles that can be used to calculate reliable template-to-target alignments even in the case of distantly related proteins (sequence identity < 30%). Other BAR+ based applications have been developed during my doctorate including the prediction of Magnesium binding sites in human proteins, the ABC transporters superfamily classification and the functional prediction (GO terms) of the CAFA targets. Remarkably, in the CAFA assessment, BAR+ placed among the ten most accurate methods. At present, as a web server for the functional and structural protein sequence annotation, BAR+ is freely available at http://bar.biocomp.unibo.it/bar2.0.

A genome-to-genome analysis of associations between human genetic variation, HIV-1 sequence diversity, and viral control

HIV-1 sequence diversity is affected by selection pressures arising from host genomic factors. Using paired human and viral data from 1071 individuals, we ran >3000 genome-wide scans, testing for associations between host DNA polymorphisms, HIV-1 sequence variation and plasma viral load (VL), while considering human and viral population structure. We observed significant human SNP associations to a total of 48 HIV-1 amino acid variants (p<2.4 × 10−12). All associated SNPs mapped to the HLA class I region. Clinical relevance of host and pathogen variation was assessed using VL results. We identified two critical advantages to the use of viral variation for identifying host factors: (1) association signals are much stronger for HIV-1 sequence variants than VL, reflecting the ‘intermediate phenotype’ nature of viral variation; (2) association testing can be run without any clinical data. The proposed genome-to-genome approach highlights sites of genomic conflict and is a strategy generally applicable to studies of host–pathogen interaction.