The role of gene duplication in the origin and evolution of new biological functions

Abstract : Gene duplication is an essential source of material for the origin of genetic novelty and the evolution of lineage- or species-specific phenotypic traits. The reverse transcription of source gene mRNA followed by the genomic insertion of the resulting cDNA - retroposition - has provided the human genome with a significant number of gene copies during the last ~63 million years (MYA) of primate evolution. We estimated that at least 1 new functional gene (retrogene) per MYA emerged by retroposition in the primate lineage leading to humans. Using a combination of comparative sequencing and evolutionary simulations, we obtained strong evidence of functionality for 7 primate specific retrogenes. Most of these genes are specifically expressed in testis suggesting that retroposition has contributed with genetic raw material necessary for the evolution ofmale-specific functions in primates. We characterized CDC14Bretro (identified in the previous survey) that originated from the retroposition of a cell cycle gene - CDC14B - in the common ancestor of humans and apes. We demonstrate that CDC14Bretro experienced a period of intense positive selection in the African ape ancestor. By virtue of the amino acid substitutions that occurred during this period CDC 14Bretro adapted to a new subcellular compartment in African apes. Further analyses indicate that this subcellular shift reflects the evolution of anew functional role of CDC 14Bretro. Prompted by this result, we used yeast (Saccharomyces cerevisiae) to investigate on a global scale the extent of functional diversification of duplicate genes through the subcellular adaptation of their encoded proteins. We found that duplicate proteins frequently evolved new cellular localization patterns, either by partitioning of ancestral localizations ("sublocalization"), or more frequently by relocalization to previously unoccupied compartments ("neolocalization"). Interestingly, proteins involved in processes with a wider subcellular distribution more frequently evolved new localization patterns suggesting that subcellular localization changes are dependent on progenitor gene functions. Relocated proteins adapted to their new subcellular environments and evolved new functional roles through changes of their physio-chemical properties, expression levels, and interaction partners. Our work suggests an important role of subcellular adaptation for the emergence of new gene functions.

Increased therapeutic window for the R91W mutant form of Rpe65 compared to Rpe65 null backgroud

Purpose:Given the advances of gene therapy studies to cure RPE65-derived Leber Congenital Amaurosis (LCA) (clinical trials phase I) and the heterogeneity of the targeted patients both genetically and phenotypically, it is of prime importance to examine the rescue efficiency of gene transfer in different mutant contexts. Indeed, half of these mutations are missense mutations, leading to potential residual RPE65 activity. Consequently, we wanted to evaluate the effect on retinal activity and cone survival of lentivirus-mediated gene therapy in the R91W knock-in mouse model expressing the mutant Rpe65R91W gene (Samardzija et al. 2008), a mutation found in LCA patients. Notably we investigated whether if the therapeutic window is prolonged in comparison to null mutations. Methods:An HIV-1-derived lentiviral vector (LV) expressing either the GFP or the mouse Rpe65 cDNA under the control of a 0.8 kb fragment of the human Rpe65 promoter (R0.8) was produced by transient transfection of 293T cells. LV-R0.8-RPE65 or GFP was injected into 5-days-old (P5) or 1 month-old R91W mice. Functional rescue was assessed by ERG (1 and 4 months post-injection) and pupillary light response (PLR) recordings and cone survival by histological analysis. Results:Increased light sensitivity was detected by scotopic ERG in animals injected with LV-R0.8-RPE65 at both P5 and 1 month compared to GFP-treated animals or untreated mice. PLR was also improved in some eyes and histological analysis of cone markers showed that the density of cones reached the wild type level in the region of wt RPE65 delivery after treatment at P5. However, the rescue effect of the injection at 1 month was limited and attained 60% of the wild type level, but still more cones were observed in the treated area than in 1 month-old untreated Rpe65R91W mice. Conclusions:We were able to show that lentivirus-mediated Rpe65 gene transfer not only increases retinal activity of the Rpe65R91W mouse and survival of cones after treatment at P5 but also after treatment at 1 month. However even if the treatment at 1 month is more limited (60% of the wild type level) than treatment at P5, the amount of cone markers is increased compared to the proportion found at 1 month of age in untreated animals. This results contrast with the lack of cone rescue by treatment at 1 month of age in Rpe65-/- (Bemelmans et al, 2006). Thus patient suffering from R91W mutation might benefit from a prolonged therapeutic window.

Compositions in life science data

The aim of this talk is to convince the reader that there are a lot of interesting statisticalproblems in presentday life science data analysis which seem ultimately connected withcompositional statistics.Key words: SAGE, cDNA microarrays, (1D-)NMR, virus quasispecies

Bases moléculaires de la dysplasie pseudorhumatoïde progressive

La dysplasie pseudorhumatoïde progressive (PPRD) est une arthropathie non- inflammatoire causée par des mutations récessives du gène WISP3. Elle se manifeste pendant la petite enfance avec une raideur progressive des articulations et des douleurs. Les patients sont référées habituellement aux rhumatologues pédiatres et aux orthopédistes, et seulement dans un deuxième temps, vers les généticiens et/ou experts dans les dysplasies osseuses. Pour cette raison le diagnostic clinique, qui repose sur les signes radiologiques typiques et l'expérience clinique, est souvent retardé et les patients peuvent recevoir des traitements anti-inflammatoires et immunosuppresseurs innécessaires. Nous reportons ici une large série de patients atteints de PPRD avec confirmation du diagnostic au niveau moléculaire, et nous soulignons les caractéristiques cliniques et radiologiques de la maladie. Il existe une fenêtre d'âge dans laquelle les signes radiologiques sont spécifiquement reconnaissables. Des anomalies spondyloépiphyseales très légères peuvent apparaître avant l'âge de 9 ans et une perte de cartilage non-spécifique qui ressemble à une ostéoarthrite avancée est présente chez les jeunes adultes. Les articulations interphalangiennes sont les premières à être atteintes, suivis par les genoux et les hanches. L'atteint de la colonne arrive chez les enfants plus grands et dans l'adolescence. Une cyphose est fréquente et une scoliose survienne chez une minorité de patients. Des signes d'ostéoarthrite avancée comme des formations ostéophytiques et/ou des calcifications périarticulaires sont observés chez les jeunes adultes atteints de PPRD et peuvent être responsables d'une certaine inflammation secondaire. L'analyse moléculaire du gène WISP3 par séquençage de l'ADN génomique permet de confirmer le diagnostic dans la plupart des cas. Néanmoins, des splicings alternatives causés par des mutations introniques peuvent être détectés dans le cDNA des fibroblastes. Dans le cas où l'analyse génomique ne montre aucune mutation chez un individu présentant les signes et symptômes typiques de la maladie, une biopsie de peau est indiquée pour analyse moléculaire du cDNA. La prise en charge de la PPRD est symptomatique et largement insatisfaisante. Le remplacement des articulations les plus atteintes est souvent nécessaire dès l'adolescence afin de diminuer les douleurs et maintenir la mobilité.

Analysis of photoreceptor abnormality in GUCY2D E837D/R838S transgenic pigs

Purpose: We generated genetically engineered pigs expressing the human dominant GUCY2DE837D/R838S allele to modelize cone dystrophy. After a functional follow-up showing reduced photopic ERG responses (ARVO 2011), we analyzed the eyes by immunohistochemistry and revealed retinal modifications in the transgenic group. Methods: Lentiviral vectors encoding the human double mutant GUCY2DE837D/R838S cDNA under the control of a portion of the pig arrestin-3 promoter (Arr3) were produced and used for lentiviral-mediated transgenesis in pigs. Animals were regularly submitted to behavioral and functional investigations and were sacrificed at 4, 7, 15 and 18 months of age for histological and RT-PCR analyses. Retinal markers were used to evaluate the retinal status of eleven transgenic pigs and 6 non-transgenic controls. The expression of the mutant cDNA was also assayed by RT-PCR. Results: A significant increase in the number of displaced nuclei in the outersegment layer is observed in transgenic animals compared to control animals independently of their age. Part of these nuclei originate from cones as demonstrated by colocalization with cone markers. No significant change in the ONL thickness (central and peripheral retina) was measured between 4 and 18 months of age, showing a slow progression of the disease in the transgenic pigs within this time-frame. Conclusions: Arr3-GUCY2DE837D/R838S pigs show signs of retinal abnormality with slow progression which parallels the loss of photopic function. Further characterization of this model should help to elucidate the molecular mechanisms underlying the disease evolution.

Gene expression changes in chronic inflammatory demyelinating polyneuropathy skin biopsies.

Chronic-inflammatory demyelinating polyneuropathy (CIDP) is an immune-mediated disease with no known biomarkers for diagnosing the disease or assessing its prognosis. We performed transcriptional profiling microarray analysis on skin punch biopsies from 20 CIDP patients and 17 healthy controls to identify disease-associated gene expression changes. We demonstrate changes in expression of genes involved in immune and chemokine regulation, growth and repair. We also found a combination of two upregulated genes that can be proposed as a novel biomarker of the disorder.

Definition of clinically distinct molecular subtypes in estrogen receptor-positive breast carcinomas through genomic grade.

PURPOSE: A number of microarray studies have reported distinct molecular profiles of breast cancers (BC), such as basal-like, ErbB2-like, and two to three luminal-like subtypes. These were associated with different clinical outcomes. However, although the basal and the ErbB2 subtypes are repeatedly recognized, identification of estrogen receptor (ER) -positive subtypes has been inconsistent. Therefore, refinement of their molecular definition is needed. MATERIALS AND METHODS: We have previously reported a gene expression grade index (GGI), which defines histologic grade based on gene expression profiles. Using this algorithm, we assigned ER-positive BC to either high-or low-genomic grade subgroups and compared these with previously reported ER-positive molecular classifications. As further validation, we classified 666 ER-positive samples into subtypes and assessed their clinical outcome. RESULTS: Two ER-positive molecular subgroups (high and low genomic grade) could be defined using the GGI. Despite tracking a single biologic pathway, these were highly comparable to the previously described luminal A and B classification and significantly correlated to the risk groups produced using the 21-gene recurrence score. The two subtypes were associated with statistically distinct clinical outcome in both systemically untreated and tamoxifen-treated populations. CONCLUSION: The use of genomic grade can identify two clinically distinct ER-positive molecular subtypes in a simple and highly reproducible manner across multiple data sets. This study emphasizes the important role of proliferation-related genes in predicting prognosis in ER-positive BC.

Database construction and Peptide identification strategies for proteogenomic studies on sequenced genomes.

Since the advent of high-throughput DNA sequencing technologies, the ever-increasing rate at which genomes have been published has generated new challenges notably at the level of genome annotation. Even if gene predictors and annotation softwares are more and more efficient, the ultimate validation is still in the observation of predicted gene product( s). Mass-spectrometry based proteomics provides the necessary high throughput technology to show evidences of protein presence and, from the identified sequences, confirmation or invalidation of predicted annotations. We review here different strategies used to perform a MS-based proteogenomics experiment with a bottom-up approach. We start from the strengths and weaknesses of the different database construction strategies, based on different genomic information (whole genome, ORF, cDNA, EST or RNA-Seq data), which are then used for matching mass spectra to peptides and proteins. We also review the important points to be considered for a correct statistical assessment of the peptide identifications. Finally, we provide references for tools used to map and visualize the peptide identifications back to the original genomic information.

Protection of pancreatic beta-cells by high density lipoproteins

SUMMARYThe incidence of type 2 diabetes (T2D) is increasing worldwide and is linked to the enhancement of obesity. The principal cause of T2D development is insulin resistance, which lead to the increase of insulin production by the pancreatic beta-cells. In a pathological environment, namely dyslipidaemia, hyperglycaemia and inflammation, beta-cell compensation will fail in more vulnerable cells and diabetes will occur. High Density Lipoproteins (HDLs), commonly named "good cholesterol" are known to be atheroprotective. Low levels of HDLs are associated with increased prevalence of cardiovascular disease but are also an independent risk factor for the development of T2D. HDLs were demonstrated to protect pancreatic beta-cells against several stresses. However the molecular mechanisms of the protection are unknown and the objectives of this work were to try to elucidate the way how HDLs protect. The first approach was a broad screening of genes regulated by the stress and HDLs. A microarray analysis was performed on beta-cells stressed by serum deprivation and rescued by HDLs. Among the genes regulated, we focused on 4E-BP1, a cap-dependent translational inhibitor. In addition, HDLs were also found to protect against several other stresses.Endoplasmic reticulum (ER) stress is a mechanism that may play a role in the onset of T2D. The unfolded protein response (UPR) is a physiological process that aims at maintaining ER homeostasis in conditions where the protein folding and secretion is perturbed. Specific signalling pathways are involved in the increase of folding, export and degradation capacity of the ER. However, in case where the stress is prolonged, this mechanism turns to be pathological, by inducing cell death effector pathways, leading to beta-cell apoptosis. In our study, we discovered that HDLs were protective against ER stress induced by drugs and physiological stresses such as saturated free fatty acids. HDLs protected beta-cells by promoting ER homeostasis via the improvement of the folding and trafficking od proteins from the ER to the Golgi apparatus.Altogether our results suggest that HDLs are important for beta-cell function and survival, by protecting them from several stresses and acting on ER homeostasis. This suggests that attempt in keeping normal HDLs levels or function in patients is crucial to lessen the development of T2D.RÉSUMÉL'incidence du diabète de type 2 est en constante augmentation et est fortement liée à l'accroissement du taux d'obésité. La cause principale du diabète de type 2 est la résistance à l'insuline, qui entraîne une surproduction d'insuline par les cellules bêta pancréatiques. Dans un environnement pathologique associé à l'obésité (dyslipidémie, hyperglycémie et inflammation), les cellules bêta les plus vulnérables ne sont plus capables de compenser en augmentant leur production d'insuline, dysfonctionnent, ce qui conduit à leur mort par apoptose. Les lipoprotéines de hautes densités (HDLs), communément appelées (( bon cholestérol », sont connues pour leurs propriétés protectrices contre l'athérosclérose. Des niveaux bas de HDLs sanguins sont associés au risque de développer un diabète de type 2. Les HDLs ont également montré des propriétés protectrices contre divers stresses dans la cellule bêta. Cependant, les mécanismes de protection restent encore inconnus et l'objectif de ce travail a été d'investiguer les mécanismes moléculaires de protection des HDLs. La première approche choisie a été une étude du profil d'expression génique par puce à ADN afin d'identifier les gènes régulés par le stress et les HDLs. Parmi les gènes régulés, notre intérêt s'est porté sur 4E-BP1, un inhibiteur de la traduction coiffe- dépendante, dont l'induction par le stress était corrélée avec une augmentation de l'apoptose. Suite à cette étude, les HDLs ont également montrés un rôle protecteur contre d'autres stresses. Il s'agit particulièrement du stress du réticulum endoplasmique (RE), qui est un mécanisme qui semble jouer un rôle clé dans le développement du diabète. L'UPR (« Unfolded Protein Response ») est un processus physiologique tendant à maintenir l'homéostasie du réticulum endoplasmique, organelle prépondérante pour la fonction des cellules sécrétrices, notamment lorsqu'elle est soumise à des conditions extrêmes telles que des perturbations de la conformation tertiaire des protéines ou de la sécrétion. Dans ces cas, des voies de signalisation moléculaires sont activées, ce qui mène à l'exportation des protéines mal repliées, à leur dégradation et à l'augmentation de l'expression de chaperonnes capables d'améliorer le repliement des protéines mal formées. Toutefois, en cas de stress persistant, ce mécanisme de protection s'avère être pathologique. En induisant des voies de signalisation effectrices de l'apoptose, il conduit finalement au développement du diabète. Dans cette étude, nous avons démontré que les HDLs étaient capables de protéger la cellule bêta contre le stress du RE induits par des inhibiteurs (thapsigargine, tunicamycine) ou des stresses physiologiques tels que les acides gras libres. Les HDLs ont la capacité d'améliorer l'homéostasie du RE, notamment en favorisant le repliement et le transfert des protéines du RE à l'appareil de Golgi.En résumé, ces données suggèrent que les HDLs sont bénéfiques pour la survie des cellules bêta soumises à des stresses impliqués dans le développement du diabète, notamment en restaurant l'homéostasie du RE. Ces résultats conduisent à soutenir que le maintien des taux de cholestérol joue un rôle important dans la limitation de l'incidence du diabète.

Molecular characterization of HLA-C incompatibilities in HLA-ABDR-matched unrelated bone marrow donor-recipient pairs. Sequence of two new Cw alleles (Cw*02023 and Cw*0707) and recognition by cytotoxic T lymphocytes.

While the influence of HLA-AB and -DRB1 matching on the outcome of bone marrow transplantation (BMT) with unrelated donors is clear, the evaluation of HLA-C has been hampered by its poor serological definition. Because the low resolution of standard HLA-C typing could explain the significant number of positive cytotoxic T lymphocyte precursor frequency (CTLpf) tests found among HLA-AB-subtype, DRB1/B3/B5-subtype matched patient/donor pairs, we have identified by sequencing the incompatibilities recognized by CD8+ CTL clones obtained from such positive CTLpf tests. In most cases the target molecules were HLA-C antigens that had escaped detection by serology (e.g. Cw*1601, 1502 or 0702). Direct recognition of HLA-C by a CTL clone was demonstrated by lysis of the HLA class I-negative 721.221 cell line transfected with Cw*1601 cDNA. Because of the functional importance of Cw polymorphism, a PCR-SSO oligotyping procedure was set up allowing the resolution of 29 Cw alleles. Oligotyping of a panel of 382 individuals (including 101 patients and their 272 potential unrelated donors, 5 related donors and 4 platelet donors) allowed to determine HLA-C and HLA A-B-Cw-DRB1 allelic frequencies, as well as a number of A-Cw, B-Cw, and DRB1-Cw associations. Two new HLA-Cw alleles (Cw*02023 and Cw*0707) were identified by DNA sequencing of PCR-amplified exon 2-intron 2-exon 3 amplicons. Furthermore, we determined the degree of HLA-C compatibility in 287 matched pairs that could be formed from 73 patients and their 184 potential unrelated donors compatible for HLA-AB by serology and for HLA-DRB1/ B3/B5 by oligotyping. Cw mismatches were identified in 42.1% of these pairs, and AB-subtype oligotyping showed that 30% of these Cw-incompatible pairs were also mismatched for A or B-locus subtype. The degree of HLA-C incompatibility was strongly influenced by the linkage with B alleles and by the ABDR haplotypes. Cw alleles linked with B*4403, B*5101, B18, and B62 haplotypes were frequently mismatched. Apparently high resolution DNA typing for HLA-AB does not result in full matching at locus C. Since HLA-C polymorphism is recognized by alloreactive CTLs, such incompatibilities might be as relevant as AB-subtype mismatches in clinical transplantation.

Purification and characterization of the human Rad51 protein, an analogue of E. coli RecA.

In bacteria, genetic recombination is catalysed by RecA protein, the product of the recA gene. A human gene that shares homology with Escherichia coli recA (and its yeast homologue RAD51) has been cloned from a testis cDNA library, and its 37 kDa product (hRad51) purified to homogeneity. The human Rad51 protein binds to single- and double-stranded DNA and exhibits DNA-dependent ATPase activity. Using a topological assay, we demonstrate that hRad51 underwinds duplex DNA, in a reaction dependent upon the presence of ATP or its non-hydrolysable analogue ATP gamma S. Complexes formed with single- and double-stranded DNA have been observed by electron microscopy following negative staining. With nicked duplex DNA, hRad51 forms helical nucleoprotein filaments which exhibit the striated appearance characteristic of RecA or yeast Rad51 filaments. Contour length measurements indicate that the DNA is underwound and extended within the nucleoprotein complex. In contrast to yeast Rad51 protein, human Rad51 forms filaments with single-stranded DNA in the presence of ATP/ATP gamma S. These resemble the inactive form of the RecA filament which is observed in the absence of a nucleotide cofactor.

Complete Nucleotide Sequence of Mouse Mammary Tumor Virus from JYG Chinese Wild Mice: Absence of Bacterial Insertion Sequences in the Cloned Viral gag Gene.

Mammary tumors of a newly isolated strain of Chinese wild mouse (JYG mouse) harbor exogenous mouse mammary tumor virus (MMTV). The complete nucleotide sequence of exogenous JYG-MMTV was determined on the proviral 5' long terminal repeat (LTR)(partial)-gag-pol-env-3' LTR (partial) fragment cloned into a plasmid vector and the cDNA sequence from JYG-MMTV producing cells. Similarly to the other MMTV species the LTR of JYG-MMTV contains an open reading frame (ORF). The amino acid sequence of the JYG-MMTV ORF resembles that of SW-MMTV (92% identity) and endogenous Mtv-7 (93% identity) especially at the C-terminal region. Thus, a functional similarity in T-cell receptor V beta recognition as a superantigen is implicated among these MMTV species. Analysis of the viral gag nucleotide sequence revealed that this gene is not disrupted by the bacterial insertion sequence IS1 or IS2, which have been reported to be present in the majority of the plasmids containing the gag region. Comparison of amino acid sequences of JYG-MMTV with those of BR6-MMTV showed that over 96% of the amino acids of gag, pol, protease and env products are identical. These results suggest the intact nature of the nucleotide sequence of the near full-length MMTV genome cloned in the plasmid.

Transepithelial transport of HIV-1 by M cells is receptor-mediated.

Human colon carcinoma Caco-2 cell monolayers undergo conversion into cells that share morphological and functional features of M cells when allowed to interact with B lymphocytes. A lymphotropic (X4) HIV-1 strain crosses M cell monolayers and infects underlying CD4(+) target cells. Transport requires both lactosyl cerebroside and CXCR4 receptors, which are expressed on the apical surface of Caco-2 and M cells. Antibodies specific for each receptor block transport. In contrast, a monotropic (R5) HIV-1 strain is unable to cross M cell monolayers and infect underlying monocytes, despite efficient transport of latex beads. Caco-2 and M cells do not express CCR5, but transfection of these cells with CCR5 cDNA restores transport of R5 virus, which demonstrates that HIV-1 transport across M cells is receptor-mediated. The follicle-associated epithelium covering human gut lymphoid follicles expresses CCR5, but not CXCR4, and lactosyl cerebroside, suggesting that HIV-1 infection may occur through M cells and enterocytes at these sites.

Recurrent deletions and reciprocal duplications of 10q11.21q11.23 including CHAT and SLC18A3 are likely mediated by complex low-copy repeats.

We report 24 unrelated individuals with deletions and 17 additional cases with duplications at 10q11.21q21.1 identified by chromosomal microarray analysis. The rearrangements range in size from 0.3 to 12 Mb. Nineteen of the deletions and eight duplications are flanked by large, directly oriented segmental duplications of >98% sequence identity, suggesting that nonallelic homologous recombination (NAHR) caused these genomic rearrangements. Nine individuals with deletions and five with duplications have additional copy number changes. Detailed clinical evaluation of 20 patients with deletions revealed variable clinical features, with developmental delay (DD) and/or intellectual disability (ID) as the only features common to a majority of individuals. We suggest that some of the other features present in more than one patient with deletion, including hypotonia, sleep apnea, chronic constipation, gastroesophageal and vesicoureteral refluxes, epilepsy, ataxia, dysphagia, nystagmus, and ptosis may result from deletion of the CHAT gene, encoding choline acetyltransferase, and the SLC18A3 gene, mapping in the first intron of CHAT and encoding vesicular acetylcholine transporter. The phenotypic diversity and presence of the deletion in apparently normal carrier parents suggest that subjects carrying 10q11.21q11.23 deletions may exhibit variable phenotypic expressivity and incomplete penetrance influenced by additional genetic and nongenetic modifiers.

Phenotypic switching in Pseudomonas brassicacearum involves GacS- and GacA-dependent Rsm small RNAs.

The plant-beneficial bacterium Pseudomonas brassicacearum forms phenotypic variants in vitro as well as in planta during root colonization under natural conditions. Transcriptome analysis of typical phenotypic variants using microarrays containing coding as well as noncoding DNA fragments showed differential expression of several genes relevant to secondary metabolism and of the small RNA (sRNA) genes rsmX, rsmY, and rsmZ. Naturally occurring mutations in the gacS-gacA system accounted for phenotypic switching, which was characterized by downregulation of antifungal secondary metabolites (2,4-diacetylphloroglucinol and cyanide), indoleacetate, exoenzymes (lipase and protease), and three different N-acyl-homoserine lactone molecules. Moreover, in addition to abrogating these biocontrol traits, gacS and gacA mutations resulted in reduced expression of the type VI secretion machinery, alginate biosynthesis, and biofilm formation. In a gacA mutant, the expression of rsmX was completely abolished, unlike that of rsmY and rsmZ. Overexpression of any of the three sRNAs in the gacA mutant overruled the pleiotropic changes and restored the wild-type phenotypes, suggesting functional redundancy of these sRNAs. In conclusion, our data show that phenotypic switching in P. brassicacearum results from mutations in the gacS-gacA system.