Functional diversification of duplicate genes through subcellular adaptation of encoded proteins.

BACKGROUND: Gene duplication is the primary source of new genes with novel or altered functions. It is known that duplicates may obtain these new functional roles by evolving divergent expression patterns and/or protein functions after the duplication event. Here, using yeast (Saccharomyces cerevisiae) as a model organism, we investigate a previously little considered mode for the functional diversification of duplicate genes: subcellular adaptation of encoded proteins. RESULTS: We show that for 24-37% of duplicate gene pairs derived from the S. cerevisiae whole-genome duplication event, the two members of the pair encode proteins that localize to distinct subcellular compartments. The propensity of yeast duplicate genes to evolve new localization patterns depends to a large extent on the biological function of their progenitor genes. Proteins involved in processes with a wider subcellular distribution (for example, catabolism) frequently evolved new protein localization patterns after duplication, whereas duplicate proteins limited to a smaller number of organelles (for example, highly expressed biosynthesis/housekeeping proteins with a slow rate of evolution) rarely relocate within the cell. Paralogous proteins evolved divergent localization patterns by partitioning of ancestral localizations ('sublocalization'), but probably more frequently by relocalization to new compartments ('neolocalization'). We show that such subcellular reprogramming may occur through selectively driven substitutions in protein targeting sequences. Notably, our data also reveal that relocated proteins functionally adapted to their new subcellular environments and evolved new functional roles through changes of their physico-chemical properties, expression levels, and interaction partners. CONCLUSION: We conclude that protein subcellular adaptation represents a common mechanism for the functional diversification of duplicate genes.

Birth and adaptive evolution of a hominoid gene that supports high neurotransmitter flux.

The enzyme glutamate dehydrogenase (GDH) is important for recycling the chief excitatory neurotransmitter, glutamate, during neurotransmission. Human GDH exists in housekeeping and brain-specific isotypes encoded by the genes GLUD1 and GLUD2, respectively. Here we show that GLUD2 originated by retroposition from GLUD1 in the hominoid ancestor less than 23 million years ago. The amino acid changes responsible for the unique brain-specific properties of the enzyme derived from GLUD2 occurred during a period of positive selection after the duplication event.

Fourteen new cases contribute to the characterization of the 7q11.23 microduplication syndrome.

Interstitial deletions of 7q11.23 cause Williams-Beuren syndrome, one of the best characterized microdeletion syndromes. The clinical phenotype associated with the reciprocal duplication however is not well defined, though speech delay is often mentioned. We present 14 new 7q11.23 patients with the reciprocal duplication of the Williams-Beuren syndrome critical region, nine familial and five de novo. These were identified by either array-based MLPA or by array-CGH/oligonucleotide analysis in a series of patients with idiopathic mental retardation with an estimated population frequency of 1:13,000-1:20,000. Variable speech delay is a constant finding in our patient group, confirming previous reports. Cognitive abilities range from normal to moderate mental retardation. The association with autism is present in five patients and in one father who also carries the duplication. There is an increased incidence of hypotonia and congenital anomalies: heart defects (PDA), diaphragmatic hernia, cryptorchidism and non-specific brain abnormalities on MRI. Specific dysmorphic features were noted in our patients, including a short philtrum, thin lips and straight eyebrows. Our patient collection demonstrates that the 7q11.23 microduplication not only causes language delay, but is also associated with congenital anomalies and a recognizable face.

Molecular evolutionary patterns in the Glomeromycota;

Abstract Arbuscular mycorrhizal fungi (AMF) form symbiosis with roots of approximately 80% of known land plants. These fungi play a key role in the ecology and adaptation of plants to various ecosystems.by increasing the plant resources for various nutrients. Despite their important ecological role, we still have poor understanding of their genetic structure and their molecular evolution. The work presented in this thesis aims to isolate and analyse AMF genes with various molecular techniques, in order to obtain new insights about their genetics, phylogeny and molecular evolution. Some AMF genes were shown through phylogenetic analyses to be more related with plants or mycoparasites than with other fungal organisms. These results led to the prediction that lateral gene transfers (LGT) occurred between AMF and plants during their long-term co-évolution. By phylogenetic and molecular analyses, in the chapter 2 I demonstrate that the hypothesis of LGT is most likely a consequence of analyses carried out on contaminant non AMF-DNA. In addition, various features characteristic of AMF genes have been determined, allowing researchers to scan their own sequence databases for potential non-AMF contaminants. Phylogenetic relationships of AMF with other fungi has been mostly analysed using molecular markers of ribosomal origin. In chapter 2 I successfully isolated gene encoding α- and ß-tubulins from several AMF genera. Consequently, phylogenetic analyses showed that AMF possess an unexpected relationship with ancestral aquatic fungi (chytrids). These results are consistent with the prediction stating that AMF may have played an important role in the colonisation of land by green plants through the establishment of a symbiosis and after the divergence of AMF from aquatic ancestors. In Chapter 4 I tried to isolate the entire AMF gene family encoding P-Type II ATPases, in order to determine their molecular evolution with the fungal kingdom. These genes were further analysed to detect the level of sequence polymorphism that is present within an AMF population. The results obtained show that mutational events previously thought as occurring only among divergent evolutionary lineages (gene duplications, indel mutations in coding regions) can occur within a single population of AMF. These results have far reaching consequences for our understanding of the genetics and ecology of AMF. Résumé Les champignons endomycorrhiziens arbusculaires (CEA) forment une symbiose racinaire avec environ 80% des plantes vasculaires connues. Ces champignons possèdent un rôle important dans l'écologie et l'adaptation des plantes au sein de différents écosystèmes en .augmentant leurs ressources en nutriments. Le travail présenté dans cette thèse se propose d'isoler et d'analyser certains gènes de CEA avec différentes techniques moléculaires à fin d'obtenir de pÌus amples informations concernant l'évolution moléculaire, la phylogénie et leur diversité génétique à diverses échelles taxonomiques. Certaines analyses phylogénétiques des CEA ont conduit à l'hypothèse que des transferts horizontaux de gènes (THG) ont pu avoir lieu durant leur longue co-évolution avec les plantes vasculaires. Dans le chapitre 2 de cette thèse nous démontrons par analyses moléculaire et phylogénétique que l'hypothèse de THG est une conséquence de contaminations à partir d'ADN de plante ou d'autres micro-organismes. De plus, de nombreuses caractéristiques moléculaires de CEA ont pu être déterminées, permettant la mise en place d'un plan à suivre lors de l'analyse de gènes de CEA dans les études futures. Les relations évolutives des. CEA avec d'autres champignons ont été analysées majoritairement à l'aide de marqueurs moléculaires d'origine ribosomiale. Dans les chapitres 2 et 3 j'ai isolé des gènes codant pour l'a- et la ß-tubuline chez différents genres, de CEA. Les analyses phylogénétiques ont démontré une parenté entre les CEA et des champignons aquatiques ancestraux (chytrides). Ces résultats sont en accord avec l'hypothèse selon laquelle les CEA ont probablement joué un rôle primordial dans l'établissement des plantes sur terre à travers une symbiose et suite à leur évolution à partir d'ancêtres vivant dans des milieux aquatiques: Dans le chapitre 4 j'ai isolé une entière famille de gènes chez les CEA codant des ATPases de la membrane plasmique, et étudié leur évolution moléculaire dans le règne des champignons. Ces mêmes gènes ont été analysés ultérieurement à fin de déterminer le degré de polymorphisme de séquence qui peut être présent au sein d'une population de CEA. Les résultats obtenus montrent que des évènements mutationnels considérés comme apparaissant exclusivement dans des lignées évolutives très divergentes (duplication de gènes, insertions/délétions dans des régions transcrites du génome) ont lieu sein d'une même population de CEA. Cette découverte a un impact important sur nos connaissances concernant la génétique des populations des CEA et leur écologie.

The mechanisms of centriole biogenesis

Résumé : Le centrosome contient une paire de centrioles entourée par du matériel péricentriolaire (PCM) et cet ensemble constitue le centre organisateur des microtubules de la majorité des cellules animales. Tout comme l'ADN, 1'unique centrosome présent au début du cycle cellulaire est dupliqué une et une seule fois pour former deux centrosomes qui vont orchestrer la mise en place du fuseau mitotique. La duplication du centrosome doit être soumise à une régulation précise car la présence d'un seul ou de plus de deux centrosomes peut entraîner la formation d'un fuseau mitotique aberrant, la mauvaise ségrégation des chromosomes et l'aneuploïdie. Bien que la duplication des centrioles soit un phénomène clé pour la duplication du centrosome lui-même, les mécanismes impliqués dans la formation des centrioles sont peu connus et constituent une importante question de biologie cellulaire. Dans cette thèse, nous nous sommes concentrés sur l'analyse de HsSAS-6. Nous avons trouvé que cette protéine est nécessaire pour la formation d'un centriole et qu'elle est localisée spécifiquement à la base des nouveaux centrioles formés. Les niveaux de HsSAS-6 oscillent pendant le cycle cellulaire : la protéine est absente en G1, commence à s'accumuler au niveau du centriole et dans le cytoplasme dès le début de la phase S de synthèse et disparaît abruptement pendant l'anaphase, où probablement APC/CCdlh1 la dirige vers une dégradation par le protéasome 26S. Il est important de noter que la surexpression de HsSAS-6 entraîne la formation de multiples centrioles au lieu d'un seul, ce qui indique que les niveaux de HsSAS-6 déterminent le nombre de centrioles formés. En plus de HsSAS-6, nous avons aussi étudié la lignée mutante sas-2 de C. elegans qui quelques fois assemble un fuseau multi-polaire dans l'embryon à une cellule. Nous avons montré que ce phénotype est la conséquence de la présence de multiples centrioles dans les cellules du sperme. Enfin, nous avons aussi préparé une palette de vecteurs compatibles avec le système Gateway pour permettre la génération rapide de lignées cellulaires humaines exprimant des protéines de manière inductible. De plus, nous avons commencé à développer une méthode pour évaluer la duplication des centrioles par le biais d'une plateforme de criblage d'une librairie de siRNA humains. Dans l'ensemble, notre travail a pu apporter une nouvelle compréhension du processus de duplication des centrioles et a contribué au développement de nouveaux outils de recherche de ce processus. Summary : Centrosomes contain a pair of centrioles surrounded by pericentriolar material (PCM) and serve as the main microtubule organizing centers (MTOCs) of most animal cells. Just like the DNA, the single centrosome present early in the cell cycle duplicates once and only once to give rise to two centrosomes which will then direct assembly of a bipolar spindle. Centrosome duplication must be precisely regulated because the presence of either one or more than two centrosomes can lead to the assembly of an aberrant spindle, chromosome missegregation and aneuploidy. Although duplication of centrioles is key for that of the entire centrosome, the mechanisms underlying centriole formation are poorly understood and represent an important question in cell biology. In this thesis, we focused on the analysis of HsSAS-6. We found that this protein is required for centriole formation and that it is localized specifically at the base of newly forming centrioles. The levels of HsSAS-6 oscillate across the cell cycle. The protein is absent during G1, starts to accumulate at the centriole and in the cytoplasm at the onset of S phase and disappears abruptly during anaphase when it is targeted for 26S proteasome dependent degradation probably by the APC/CCdh1. Importantly, overexpression of HsSAS-6 leads to the formation of multiple centrioles instead of just one, indicating that levels of HsSAS-6 determine the number of centrioles at each cell cycle. Besides HsSAS-6 that is the main focus of this thesis, we have also investigated the C. elegans mutant strain sas-2, which sometimes assembles a multipolar spindle in the one cell stage embryo. We have shown that this phenotype derives from the presence of multiple centrioles in sperm cells. Moreover, we prepared a set of Gateway compatible vectors for fast generation of human cell lines with inducible protein expression. Finally, we started to develop an assay for centriole duplication that can be used in a high throughput setting for screening of human siRNA libraries. Taken together, our work brought novel insights into the process of centriole duplication and lead to the development of new tools for further investigation of this process.

Adaptive divergence of ancient gene duplicates in the avian MHC class II beta.

Gene duplication and neofunctionalization are known to be important processes in the evolution of phenotypic complexity. They account for important evolutionary novelties that confer ecological adaptation, such as the major histocompatibility complex (MHC), a multigene family crucial to the vertebrate immune system. In birds, two MHC class II β (MHCIIβ) exon 3 lineages have been recently characterized, and two hypotheses for the evolutionary history of MHCIIβ lineages were proposed. These lineages could have arisen either by 1) an ancient duplication and subsequent divergence of one paralog or by 2) recent parallel duplications followed by functional convergence. Here, we compiled a data set consisting of 63 MHCIIβ exon 3 sequences from six avian orders to distinguish between these hypotheses and to understand the role of selection in the divergent evolution of the two avian MHCIIβ lineages. Based on phylogenetic reconstructions and simulations, we show that a unique duplication event preceding the major avian radiations gave rise to two ancestral MHCIIβ lineages that were each likely lost once later during avian evolution. Maximum likelihood estimation shows that following the ancestral duplication, positive selection drove a radical shift from basic to acidic amino acid composition of a protein domain facing the α-chain in the MHCII α β-heterodimer. Structural analyses of the MHCII α β-heterodimer highlight that three of these residues are potentially involved in direct interactions with the α-chain, suggesting that the shift following duplication may have been accompanied by coevolution of the interacting α- and β-chains. These results provide new insights into the long-term evolutionary relationships among avian MHC genes and open interesting perspectives for comparative and population genomic studies of avian MHC evolution.

Membrane association of the RNA-dependent RNA polymerase is essential for hepatitis C virus RNA replication.

The hepatitis C virus (HCV) RNA-dependent RNA polymerase (RdRp), represented by nonstructural protein 5B (NS5B), belongs to a class of integral membrane proteins termed tail-anchored proteins. Its membrane association is mediated by the C-terminal 21 amino acid residues, which are dispensable for RdRp activity in vitro. For this study, we investigated the role of this domain, termed the insertion sequence, in HCV RNA replication in cells. Based on a structural model and the amino acid conservation among different HCV isolates, we designed a panel of insertion sequence mutants and analyzed their membrane association and RNA replication. Subgenomic replicons with a duplication of an essential cis-acting replication element overlapping the sequence that encodes the C-terminal domain of NS5B were used to unequivocally distinguish RNA versus protein effects of these mutations. Our results demonstrate that the membrane association of the RdRp is essential for HCV RNA replication. Interestingly, certain amino acid substitutions within the insertion sequence abolished RNA replication without affecting membrane association, indicating that the C-terminal domain of NS5B has functions beyond serving as a membrane anchor and that it may be involved in critical intramembrane protein-protein interactions. These results have implications for the functional architecture of the HCV replication complex and provide new insights into the expanding spectrum of tail-anchored proteins.

Evolutionary patterns of MHC class II B in owls and their implications for the understanding of avian MHC evolution

Owing to its special mode of evolution and central role in the adaptive immune system, the major histocompatibility complex (MHC) has become the focus of diverse disciplines such as immunology, evolutionary ecology, and molecular evolution. MHC evolution has been studied extensively in diverse vertebrate lineages over the last few decades, and it has been suggested that birds differ from the established mammalian norm. Mammalian MHC genes evolve independently, and duplication history (i.e., orthology) can usually be traced back within lineages. In birds, this has been observed in only 3 pairs of closely related species. Here we report strong evidence for the persistence of orthology of MHC genes throughout an entire avian order. Phylogenetic reconstructions of MHC class II B genes in 14 species of owls trace back orthology over tens of thousands of years in exon 3. Moreover, exon 2 sequences from several species show closer relationships than sequences within species, resembling transspecies evolution typically observed in mammals. Thus, although previous studies suggested that long-term evolutionary dynamics of the avian MHC was characterized by high rates of concerted evolution, resulting in rapid masking of orthology, our results question the generality of this conclusion. The owl MHC thus opens new perspectives for a more comprehensive understanding of avian MHC evolution.

The evolutionary history of the CD209 (DC-SIGN) family in humans and non-human primates

The CD209 gene family that encodes C-type lectins in primates includes CD209 (DC-SIGN), CD209L (L-SIGN) and CD209L2. Understanding the evolution of these genes can help understand the duplication events generating this family, the process leading to the repeated neck region and identify protein domains under selective pressure. We compiled sequences from 14 primates representing 40 million years of evolution and from three non-primate mammal species. Phylogenetic analyses used Bayesian inference, and nucleotide substitutional patterns were assessed by codon-based maximum likelihood. Analyses suggest that CD209 genes emerged from a first duplication event in the common ancestor of anthropoids, yielding CD209L2 and an ancestral CD209 gene, which, in turn, duplicated in the common Old World primate ancestor, giving rise to CD209L and CD209. K(A)/K(S) values averaged over the entire tree were 0.43 (CD209), 0.52 (CD209L) and 0.35 (CD209L2), consistent with overall signatures of purifying selection. We also assessed the Toll-like receptor (TLR) gene family, which shares with CD209 genes a common profile of evolutionary constraint. The general feature of purifying selection of CD209 genes, despite an apparent redundancy (gene absence and gene loss), may reflect the need to faithfully recognize a multiplicity of pathogen motifs, commensals and a number of self-antigens

Gene organization of the mating type regions in the ectomycorrhizal fungus Laccaria bicolor reveals distinct evolution between the two mating type loci

In natural conditions, basidiomycete ectomycorrhizal fungi such as Laccaria bicolor are typically in the dikaryotic state when forming symbioses with trees, meaning that two genetically different individuals have to fuse or 'mate'. Nevertheless, nothing is known about the molecular mechanisms of mating in these ecologically important fungi. Here, advantage was taken of the first sequenced genome of the ectomycorrhizal fungus, Laccaria bicolor, to determine the genes that govern the establishment of cell-type identity and orchestrate mating. The L. bicolor mating type loci were identified through genomic screening. The evolutionary history of the genomic regions that contained them was determined by genome-wide comparison of L. bicolor sequences with those of known tetrapolar and bipolar basidiomycete species, and by phylogenetic reconstruction of gene family history. It is shown that the genes of the two mating type loci, A and B, are conserved across the Agaricales, but they are contained in regions of the genome with different evolutionary histories. The A locus is in a region where the gene order is under strong selection across the Agaricales. By contrast, the B locus is in a region where the gene order is likely under a low selection pressure but where gene duplication, translocation and transposon insertion are frequent.

Resistance to nucleoside analog reverse transcriptase inhibitors mediated by human immunodeficiency virus type 1 p6 protein.

Resistance of human immunodeficiency virus type 1 (HIV-1) to antiretroviral agents results from target gene mutation within the pol gene, which encodes the viral protease, reverse transcriptase (RT), and integrase. We speculated that mutations in genes other that the drug target could lead to drug resistance. For this purpose, the p1-p6(gag)-p6(pol) region of HIV-1, placed immediately upstream of pol, was analyzed. This region has the potential to alter Pol through frameshift regulation (p1), through improved packaging of viral enzymes (p6(Gag)), or by changes in activation of the viral protease (p6(Pol)). Duplication of the proline-rich p6(Gag) PTAP motif, necessary for late viral cycle activities, was identified in plasma virus from 47 of 222 (21.2%) patients treated with nucleoside analog RT inhibitor (NRTI) antiretroviral therapy but was identified very rarely from drug-naïve individuals. Molecular clones carrying a 3-amino-acid duplication, APPAPP (transframe duplication SPTSPT in p6(Pol)), displayed a delay in protein maturation; however, they packaged a 34% excess of RT and exhibited a marked competitive growth advantage in the presence of NRTIs. This phenotype is reminiscent of the inoculum effect described in bacteriology, where a larger input, or a greater infectivity of an organism with a wild-type antimicrobial target, leads to escape from drug pressure and a higher MIC in vitro. Though the mechanism by which the PTAP region participates in viral maturation is not known, duplication of this proline-rich motif could improve assembly and packaging at membrane locations, resulting in the observed phenotype of increased infectivity and drug resistance.

Whole genome sequencing in patients with retinitis pigmentosa reveals pathogenic DNA structural changes and NEK2 as a new disease gene.

We performed whole genome sequencing in 16 unrelated patients with autosomal recessive retinitis pigmentosa (ARRP), a disease characterized by progressive retinal degeneration and caused by mutations in over 50 genes, in search of pathogenic DNA variants. Eight patients were from North America, whereas eight were Japanese, a population for which ARRP seems to have different genetic drivers. Using a specific workflow, we assessed both the coding and noncoding regions of the human genome, including the evaluation of highly polymorphic SNPs, structural and copy number variations, as well as 69 control genomes sequenced by the same procedures. We detected homozygous or compound heterozygous mutations in 7 genes associated with ARRP (USH2A, RDH12, CNGB1, EYS, PDE6B, DFNB31, and CERKL) in eight patients, three Japanese and five Americans. Fourteen of the 16 mutant alleles identified were previously unknown. Among these, there was a 2.3-kb deletion in USH2A and an inverted duplication of ∼446 kb in EYS, which would have likely escaped conventional screening techniques or exome sequencing. Moreover, in another Japanese patient, we identified a homozygous frameshift (p.L206fs), absent in more than 2,500 chromosomes from ethnically matched controls, in the ciliary gene NEK2, encoding a serine/threonine-protein kinase. Inactivation of this gene in zebrafish induced retinal photoreceptor defects that were rescued by human NEK2 mRNA. In addition to identifying a previously undescribed ARRP gene, our study highlights the importance of rare structural DNA variations in Mendelian diseases and advocates the need for screening approaches that transcend the analysis of the coding sequences of the human genome.

Secreted subtilisin gene family in Trichophyton rubrum.

Secreted proteases constitute potential virulence factors of dermatophytes. A total of seven genes encoding putative serine proteases of the subtilisin family (SUB) were isolated in Trichophyton rubrum. Based on sequence data and intron-exon structure, a phylogenetic analysis of subtilisins from T. rubrum and other fungi revealed a presumed ancestral lineage comprising T. rubrum SUB2 and Aspergillus SUBs. All other SUBs (SUB1, SUB3-7) are dermatophyte-specific and have apparently emerged more recently, through successive gene duplication events. We showed that two subtilisins, Sub3 and Sub4, were detected in culture supernatants of T. rubrum grown in a medium containing soy protein as a sole nitrogen source. Both recombinant enzymes produced in Pichia pastoris are highly active on keratin azure suggesting that these proteases play an important role in invasion of keratinised tissues by the fungus. The set of deduced amino acid sequences of T. rubrum SUB ORFs allowed the identification of orthologous Subs secreted by other dermatophyte species using proteolysis and mass spectrometry.

Sex chromosomes and male functions: where do new genes go?

The position of a gene in the genome may have important consequences for its function. Therefore, when a new duplicate gene arises, its location may be critical in determining its fate. Our recent work in humans, mouse, and Drosophila provided a test by studying the patterns of duplication in sex chromosome evolution. We revealed a bias in the generation and recruitment of new gene copies involving the X chromosome that has been shaped largely by selection for male germline functions. The gene movement patterns we observed reflect an ongoing process as some of the new genes are very young while others were present before the divergence of humans and mouse. This suggests a continuing redistribution of male-related genes to achieve a more efficient allocation of male functions. This notion should be further tested in organisms employing other sex determination systems or in organisms differing in germline sex chromosome inactivation. It is likely that the selective forces that were detected in these studies are also acting on other types of duplicate genes. As a result, future work elucidating sex chromosome differentiation by other mutational mechanisms will shed light on this important process.

Comparative genomics of chemosensory protein genes reveals rapid evolution and positive selection in ant-specific duplicates.

Gene duplications can have a major role in adaptation, and gene families underlying chemosensation are particularly interesting due to their essential role in chemical recognition of mates, predators and food resources. Social insects add yet another dimension to the study of chemosensory genomics, as the key components of their social life rely on chemical communication. Still, chemosensory gene families are little studied in social insects. Here we annotated chemosensory protein (CSP) genes from seven ant genomes and studied their evolution. The number of functional CSP genes ranges from 11 to 21 depending on species, and the estimated rates of gene birth and death indicate high turnover of genes. Ant CSP genes include seven conservative orthologous groups present in all the ants, and a group of genes that has expanded independently in different ant lineages. Interestingly, the expanded group of genes has a differing mode of evolution from the orthologous groups. The expanded group shows rapid evolution as indicated by a high dN/dS (nonsynonymous to synonymous changes) ratio, several sites under positive selection and many pseudogenes, whereas the genes in the seven orthologous groups evolve slowly under purifying selection and include only one pseudogene. These results show that adaptive changes have played a role in ant CSP evolution. The expanded group of ant-specific genes is phylogenetically close to a conservative orthologous group CSP7, which includes genes known to be involved in ant nestmate recognition, raising an interesting possibility that the expanded CSPs function in ant chemical communication.