Splicing and the evolution of proteins in mammals.

It is often supposed that a protein's rate of evolution and its amino acid content are determined by the function and anatomy of the protein. Here we examine an alternative possibility, namely that the requirement to specify in the unprocessed RNA, in the vicinity of intron-exon boundaries, information necessary for removal of introns (e.g., exonic splice enhancers) affects both amino acid usage and rates of protein evolution. We find that the majority of amino acids show skewed usage near intron-exon boundaries, and that differences in the trends for the 2-fold and 4-fold blocks of both arginine and leucine show this to be owing to effects mediated at the nucleotide level. More specifically, there is a robust relationship between the extent to which an amino acid is preferred/avoided near boundaries and its enrichment/paucity in splice enhancers. As might then be expected, the rate of evolution is lowest near intron-exon boundaries, at least in part owing to splice enhancers, such that domains flanking intron-exon junctions evolve on average at under half the rate of exon centres from the same gene. In contrast, the rate of evolution of intronless retrogenes is highest near the domains where intron-exon junctions previously resided. The proportion of sequence near intron-exon boundaries is one of the stronger predictors of a protein's rate of evolution in mammals yet described. We conclude that after intron insertion selection favours modification of amino acid content near intron-exon junctions, so as to enable efficient intron removal, these changes then being subject to strong purifying selection even if nonoptimal for protein function. Thus there exists a strong force operating on protein evolution in mammals that is not explained directly in terms of the biology of the protein.

Mutations in FGF17, IL17RD, DUSP6, SPRY4, and FLRT3 Are Identified in Individuals with Congenital Hypogonadotropic Hypogonadism.

Congenital hypogonadotropic hypogonadism (CHH) and its anosmia-associated form (Kallmann syndrome [KS]) are genetically heterogeneous. Among the >15 genes implicated in these conditions, mutations in FGF8 and FGFR1 account for ∼12% of cases; notably, KAL1 and HS6ST1 are also involved in FGFR1 signaling and can be mutated in CHH. We therefore hypothesized that mutations in genes encoding a broader range of modulators of the FGFR1 pathway might contribute to the genetics of CHH as causal or modifier mutations. Thus, we aimed to (1) investigate whether CHH individuals harbor mutations in members of the so-called "FGF8 synexpression" group and (2) validate the ability of a bioinformatics algorithm on the basis of protein-protein interactome data (interactome-based affiliation scoring [IBAS]) to identify high-quality candidate genes. On the basis of sequence homology, expression, and structural and functional data, seven genes were selected and sequenced in 386 unrelated CHH individuals and 155 controls. Except for FGF18 and SPRY2, all other genes were found to be mutated in CHH individuals: FGF17 (n = 3 individuals), IL17RD (n = 8), DUSP6 (n = 5), SPRY4 (n = 14), and FLRT3 (n = 3). Independently, IBAS predicted FGF17 and IL17RD as the two top candidates in the entire proteome on the basis of a statistical test of their protein-protein interaction patterns to proteins known to be altered in CHH. Most of the FGF17 and IL17RD mutations altered protein function in vitro. IL17RD mutations were found only in KS individuals and were strongly linked to hearing loss (6/8 individuals). Mutations in genes encoding components of the FGF pathway are associated with complex modes of CHH inheritance and act primarily as contributors to an oligogenic genetic architecture underlying CHH.

Neofunctionalization in vertebrates: the example of retinoic acid receptors.

Understanding the role of gene duplications in establishing vertebrate innovations is one of the main challenges of Evo-Devo (evolution of development) studies. Data on evolutionary changes in gene expression (i.e., evolution of transcription factor-cis-regulatory elements relationships) tell only part of the story; protein function, best studied by biochemical and functional assays, can also change. In this study, we have investigated how gene duplication has affected both the expression and the ligand-binding specificity of retinoic acid receptors (RARs), which play a major role in chordate embryonic development. Mammals have three paralogous RAR genes--RAR alpha, beta, and gamma--which resulted from genome duplications at the origin of vertebrates. By using pharmacological ligands selective for specific paralogues, we have studied the ligand-binding capacities of RARs from diverse chordates species. We have found that RAR beta-like binding selectivity is a synapomorphy of all chordate RARs, including a reconstructed synthetic RAR representing the receptor present in the ancestor of chordates. Moreover, comparison of expression patterns of the cephalochordate amphioxus and the vertebrates suggests that, of all the RARs, RAR beta expression has remained most similar to that of the ancestral RAR. On the basis of these results together, we suggest that while RAR beta kept the ancestral RAR role, RAR alpha and RAR gamma diverged both in ligand-binding capacity and in expression patterns. We thus suggest that neofunctionalization occurred at both the expression and the functional levels to shape RAR roles during development in vertebrates.

RPS4Y gene family evolution in primates

Background: The RPS4 gene codifies for ribosomal protein S4, a very well-conserved protein present in all kingdoms. In primates, RPS4 is codified by two functional genes located on both sex chromosomes: the RPS4X and RPS4Y genes. In humans, RPS4Y is duplicated and the Y chromosome therefore carries a third functional paralog: RPS4Y2, which presents a testis-specific expression pattern. Results: DNA sequence analysis of the intronic and cDNA regions of RPS4Y genes from species covering the entire primate phylogeny showed that the duplication event leading to the second Y-linked copy occurred after the divergence of New World monkeys, about 35 million years ago. Maximum likelihood analyses of the synonymous and non-synonymous substitutions revealed that positive selection was acting on RPS4Y2 gene in the human lineage, which represents the first evidence of positive selection on a ribosomal protein gene. Putative positive amino acid replacements affected the three domains of the protein: one of these changes is located in the KOW protein domain and affects the unique invariable position of this motif, and might thus have a dramatic effect on the protein function.Conclusion: Here, we shed new light on the evolutionary history of RPS4Y gene family, especially on that of RPS4Y2. The results point that the RPS4Y1 gene might be maintained to compensate gene dosage between sexes, while RPS4Y2 might have acquired a new function, at least in the lineage leading to humans.

Genetic evidence for a role of adiponutrin in the metabolism of apolipoprotein B-containing lipoproteins.

Adiponutrin (PNPLA3) is a predominantly liver-expressed transmembrane protein with phospholipase activity that is regulated by fasting and feeding. Recent genome-wide association studies identified PNPLA3 to be associated with hepatic fat content and liver function, thus pointing to a possible involvement in the hepatic lipoprotein metabolism. The aim of this study was to examine the association between two common variants in the adiponutrin gene and parameters of lipoprotein metabolism in 23,274 participants from eight independent West-Eurasian study populations including six population-based studies [Bruneck (n = 800), KORA S3/F3 (n = 1644), KORA S4/F4 (n = 1814), CoLaus (n = 5435), SHIP (n = 4012), Rotterdam (n = 5967)], the SAPHIR Study as a healthy working population (n = 1738) and the Utah Obesity Case-Control Study including a group of 1037 severely obese individuals (average BMI 46 kg/m2) and 827 controls from the same geographical region of Utah. We observed a strong additive association of a common non-synonymous variant within adiponutrin (rs738409) with age-, gender-, and alanine-aminotransferase-adjusted lipoprotein concentrations: each copy of the minor allele decreased levels of total cholesterol on average by 2.43 mg/dl (P = 8.87 x 10(-7)), non-HDL cholesterol levels by 2.35 mg/dl (P = 2.27 x 10(-6)) and LDL cholesterol levels by 1.48 mg/dl (P = 7.99 x 10(-4)). These associations remained significant after correction for multiple testing. We did not observe clear evidence for associations with HDL cholesterol or triglyceride concentrations. In conclusion, our study suggests that adiponutrin is involved in the metabolism of apoB-containing lipoproteins.

Defining and searching for structural motifs using DeepView/Swiss-PdbViewer.

BACKGROUND: Today, recognition and classification of sequence motifs and protein folds is a mature field, thanks to the availability of numerous comprehensive and easy to use software packages and web-based services. Recognition of structural motifs, by comparison, is less well developed and much less frequently used, possibly due to a lack of easily accessible and easy to use software. RESULTS: In this paper, we describe an extension of DeepView/Swiss-PdbViewer through which structural motifs may be defined and searched for in large protein structure databases, and we show that common structural motifs involved in stabilizing protein folds are present in evolutionarily and structurally unrelated proteins, also in deeply buried locations which are not obviously related to protein function. CONCLUSIONS: The possibility to define custom motifs and search for their occurrence in other proteins permits the identification of recurrent arrangements of residues that could have structural implications. The possibility to do so without having to maintain a complex software/hardware installation on site brings this technology to experts and non-experts alike.

Identification of GPx8 as a novel cellular substrate of the hepatitis C virus NS3-4A protease

Background: The hepatitis C virus (HCV) NS3-4A protease is not only an essential component of the viral replication complex and a prime target for a ntiviral intervention but also a key player i n the persistence and pathogenesis of HCV. It cleaves and thereby inactivates two crucial adaptor proteins in viral RNA sensing and innate immunity (MAVS and TRIF) as well as a phosphatase involved in growth factor signaling (TCPTP). T he aim of this study was to identify novel cellular substrates o f the N S3-4A protease and to investigate their role in the replication and pathogenesis of HCV. Methods: Cell lines inducibly expressing t he NS3-4A protease were analyzed in basal as well as interferon-α-stimulated states by stable isotopic l abeling using amino acids in cell culture (SILAC) coupled with protein separation and mass spectrometry. Candidates fulfilling stringent criteria for potential substrates or products of the NS3-4A protease were further i nvestigated in different experimental systems as well a s in liver biopsies from patients with chronic hepatitis C. Results: SILAC coupled with protein separation and mass spectrometry yielded > 5000 proteins of which 18 candidates were selected for further analyses. These allowed us to identify GPx8, a membrane-associated peroxidase involved in disulfide bond formation in the endoplasmic reticulum, as a n ovel cellular substrate of the H CV NS3-4A protease. Cleavage occurs at cysteine in position 11, removing the cytosolic tip of GPx8, and was observed in different experimental systems as well as in liver biopsies from patients with chronic hepatitis C. Further functional studies, involving overexpression and RNA silencing, revealed that GPx8 is a p roviral factor involved in viral particle production but not in HCV entry or HCV RNA replication. Conclusions: GPx8 is a proviral host factor cleaved by the HCV NS3-4A protease. Studies investigating the consequences of GPx8 cleavage for protein function are underway. The identification of novel cellular substrates o f the HCV N S3-4A protease should yield new insights i nto the HCV life cycle and the pathogenesis of hepatitis C and may reveal novel targets for antiviral intervention.

Ubiquitination and cysteine nitrosylation during aging and Alzheimer's disease.

Protein oxidation and ubiquitination of brain proteins are part of mechanisms that modulate protein function or that inactivate proteins and target misfolded proteins to degradation. In this study, we focused on brain aging and on mechanism involved in neurodegeneration such as events occurring in Alzheimer's disease (AD). The goal was to identify differences in nitrosylated proteins - at cysteine residues, and in the composition of ubiquinated proteins between aging and Alzheimer's samples by using a proteomic approach. A polyclonal anti-S-nitrosyl-cysteine, a mono- and a polyclonal anti-ubiquitin antibody were used for the detection of modified or ubiquitinated proteins in middle-aged and aged human entorhinal autopsy brains tissues of 14 subjects without neurological signs and 8 Alzheimer's patients. Proteins were separated by one- and two-dimensional gel electrophoresis and analyzed by Coomassie blue and immuno-blot staining. We identified that the glial fibrillary acidic and tau proteins are more ubiquitinated in brain tissues of Alzheimer's patients. Furthermore, glial fibrillary proteins were also found in nitrosylated state and further characterized by 2D Western blots and identified. Since reactive astrocytes localized prominently around senile plaques one can speculate that elements of plaques such as beta-amyloid proteins may activate surrounding glial elements and proteins.

Mutations in the human UBR1 gene and the associated phenotypic spectrum.

Johanson-Blizzard syndrome (JBS) is a rare, autosomal recessive disorder characterized by exocrine pancreatic insufficiency, typical facial features, dental anomalies, hypothyroidism, sensorineural hearing loss, scalp defects, urogenital and anorectal anomalies, short stature, and cognitive impairment of variable degree. This syndrome is caused by a defect of the E3 ubiquitin ligase UBR1, which is part of the proteolytic N-end rule pathway. Herein, we review previously reported (n = 29) and a total of 31 novel UBR1 mutations in relation to the associated phenotype in patients from 50 unrelated families. Mutation types include nonsense, frameshift, splice site, missense, and small in-frame deletions consistent with the hypothesis that loss of UBR1 protein function is the molecular basis of JBS. There is an association of missense mutations and small in-frame deletions with milder physical abnormalities and a normal intellectual capacity, thus suggesting that at least some of these may represent hypomorphic UBR1 alleles. The review of clinical data of a large number of molecularly confirmed JBS cases allows us to define minimal clinical criteria for the diagnosis of JBS. For all previously reported and novel UBR1 mutations together with their clinical data, a mutation database has been established at LOVD.

Oxidation proteomics: the role of thiol modifications

Identification of thiol modifications has gained significant importance. It is increasingly recognized that cysteines play an important role in protein function under both physiological and patho-physiological conditions. Here we reviewed different approaches that are used to identify oxidized proteins and discuss different fluorescent labeling techniques, differential two-dimensional gel electrophoresis and matrix-assisted laser desorption ionization - time of flight identification, in short MALDI-TOF. We illuminate processes that depend on protein oxidation of cysteines and we look into consequences of thiol oxidation during aging and in a variety of diseases, with a special reference to Alzheimer's disease. There is an urgent need for methods that detect specifically oxidized proteins and are able to distinguish different oxidation types.

ENFIN--A European network for integrative systems biology.

Integration of biological data of various types and the development of adapted bioinformatics tools represent critical objectives to enable research at the systems level. The European Network of Excellence ENFIN is engaged in developing an adapted infrastructure to connect databases, and platforms to enable both the generation of new bioinformatics tools and the experimental validation of computational predictions. With the aim of bridging the gap existing between standard wet laboratories and bioinformatics, the ENFIN Network runs integrative research projects to bring the latest computational techniques to bear directly on questions dedicated to systems biology in the wet laboratory environment. The Network maintains internally close collaboration between experimental and computational research, enabling a permanent cycling of experimental validation and improvement of computational prediction methods. The computational work includes the development of a database infrastructure (EnCORE), bioinformatics analysis methods and a novel platform for protein function analysis FuncNet.

Computational toolbox towards evolutionary domain mapping of membrane proteins

Membrane proteins account for about 20% to 30% of all proteins encoded in a typical genome. They play central roles in multiple cellular processes mediating the interaction of the cell with its surrounding. Over 60% of all drug targets contain a membrane domain. The experimental difficulties of obtaining a crystal structural severely limits our ability or understanding of membrane protein function. Computational evolutionary studies of proteins are crucial for the prediction of 3D structures. In this project, we construct a tool able to quantify the evolutionary positive selective pressure on each residue of membrane proteins through maximum likelihood phylogeny reconstruction. The conservation plot combined with a structural homology model is also a potent tool to predict those residues that have essentials roles in the structure and function of a membrane protein and can be very useful in the design of validation experiments.

The novel p.E89K mutation in the SRY gene inhibits DNA binding and causes the 46,XY disorder of sex development

Male sex determination in humans is controlled by the SRY gene, which encodes a transcriptional regulator containing a conserved high mobility group box domain (HMG-box) required for DNA binding. Mutations in the SRY HMG-box affect protein function, causing sex reversal phenotypes. In the present study, we describe a 19-year-old female presenting 46,XY karyotype with hypogonadism and primary amenorrhea that led to the diagnosis of 46,XY complete gonadal dysgenesis. The novel p.E89K missense mutation in the SRY HMG-box was identified as a de novo mutation. Electrophoretic mobility shift assays showed that p.E89K almost completely abolished SRY DNA-binding activity, suggesting that it is the cause of SRY function impairment. In addition, we report the occurrence of the p.G95R mutation in a 46,XY female with complete gonadal dysgenesis. According to the three-dimensional structure of the human SRY HMG-box, the substitution of the conserved glutamic acid residue by the basic lysine at position 89 introduces an extra positive charge adjacent to and between the positively charged residues R86 and K92, important for stabilizing the HMG-box helix 2 with DNA. Thus, we propose that an electrostatic repulsion caused by the proximity of these positive charges could destabilize the tip of helix 2, abrogating DNA interaction.

Étude des interactions protéiques impliquant NPM-MLF1 dans la leucémie myéloïde aiguë.

Différentes translocations génomiques sont fréquemment associées à l'apparition de leucémies myéloïdes aiguës (LMA). Ces translocations génomiques résultent de l’assemblage de deux gènes conduisant à la production d'une protéine de fusion. C'est le cas de la translocation t (3; 5) (q25.1; q34) impliquant le suppresseur tumoral NPM et l'oncogène MLF1 donnant naissance à la protéine de fusion NPM-MLF1. Généralement, les gènes impliqués dans ces translocations contrôlent la croissance cellulaire, la différenciation ou la survie cellulaire. Cependant, pour NPM-MLF1 les causes du gain ou de la perte de fonction associée à la translocation demeurent inconnues car nous ne savons pas comment cette translocation peut favoriser ou participer à l'avènement de la LMA. Le but de ce travail est d’analyser le rôle de NPM-MLF1 dans le cancer et d’examiner comment son activité contribue à la leucémie en faisant des études d’interactions protéine/protéine. En effet, l’étude de la fonction d’une protéine implique souvent de connaître ses partenaires d’interactions. Pour ce faire, la technique de double hybride dans la souche de levure AH109 a été utilisée. Tout d’abord, les ADN complémentaires (ADNc) de MLF1, NPM1 et de NPM-MLF1, MLF1-Like (une partie de MLF1 de l’acide aminé 94 à 157) normaux et mutés du domaine MTG8-Like constitué des acides aminés (a.a.) 151 à 164 de MLF1 (excepté NPM) ont été clonés dans un vecteur d'expression de levure pGBKT7. Les ADNc de GFI-1, mSin3A, PLZF, HDAC1 et HDAC3 ont été clonés dans le plasmide pGADT7 de façon à créer des protéines de fusion synthétiques avec le domaine de liaison à l'ADN et de trans-activation de la protéine GAL4. Le plasmide pGBKT7 possède un gène TRP1 et pGADT7 un gène LEU2 qui permettent la sélection des clones insérés dans la levure. Aussi, le pGBKT7 a un épitope c-myc et pGADT7 un épitope HA qui permet de voir l’expression des protéines par buvardage de type Western. Après la transformation des levures les interactions protéine/protéine ont été observées en vérifiant l’expression des gènes rapporteurs HIS3, LacZ, MEL1, ADE2 de la levure en utilisant des milieux de sélection YPD/-Leu/-Trp, YPD/-Leu/-Trp/-His, YPD/-Leu/-Trp/-His/-Ade, YPD/-Leu/-Trp/+ X-Gal, YPD/-Leu/-Trp/ + X-α-Gal. Ensuite, les interactions trouvées par double-hybride ont été vérifiées dans les cellules érythroleucémiques K562 par immuno-précipitation (IP) de protéines suivies de buvardages Westerns avec les anticorps appropriés. NPM-MLF1, MLF1, MTG8, MLF1-Like surexprimés dans les cellules K562 ont été clonés dans le plasmide pOZ-FH-N. pOZ-FH-N possède un récepteur IL-2 qui permet de sélectionner les cellules qui l’expriment ainsi qu’un tag Flag-HA qui permet de voir l’expression des protéines par buvardage-Western. Les résultats du double-hybride suggèrent une interaction faible de NPM-MLF1 avec HDAC1, HDAC3 et mSin3A ainsi qu’une interaction qui semble plus évidente entre NPM-MLF1 et PLZF, GFI-1. NPM interagit avec GFI-1 et mSin3A. Aussi, MLF1 et MLF1-Like interagissent avec HDAC1, HDAC3, GFI-1, PLZF mais pas avec mSin3A. Les IP suggèrent que NPM-MLF1 interagit avec HDAC1, HDAC3, mSin3A et PLZF. MLF1 et MLF1-Like interagissent avec HDAC1, HDAC3 et mSin3A. L’interaction de NPM-MLF1 avec GFI-1, MLF1 et MLF1-Like avec PLZF et GFI-1 n’a pas encore été vérifiée par IP. Ainsi, nos observations permettent de suggérer que NPM-MF1, MLF1 et NPM pourraient jouer un rôle dans la transcription et la régulation de l’expression de certains gènes importants dans l’hématopoïèse et une variété de processus cellulaires parce qu’ils interagissent avec différents corépresseurs. En déterminant les partenaires protéiques de MLF1, NPM et NPM-MLF1, leurs fonctions et comment NPM-MLF1 influence et modifie le fonctionnement cellulaire normal; il sera possible de renverser le processus de LMA favorisé par la t (3; 5) NPM-MLF1 par la technologie d’interférence à l’ARN.

Modulation de l’expression du gène CFTR par le produit du gène FIC1 responsable de la cholestase familiale intra-hépatique progressive de type 1 : Identification des mécanismes moléculaires impliqués

Réalisé en cotutelle avec l'Université de Cergy-Pontoise