Characterization of the expression, regulation and function of Sorting Nexin 9 and Sorting Nexin 18

Summary : Sorting nexin (SNX) family members play important roles in intracellular protein and membrane trafficking, The membrane-tubulating SNX9 protein has been shown to interact with multiple components of the endocytic machinery and to participate in clathrin-mediated endocytosis of cell surface receptors. It has not been investigated if SNX9 may also participate in other protein sorting pathways that involve vesicular transport, specifically the biogenesis of lysosome-related organelles (LROs). Closely related to SNX9 is SNXl8, whose function is largely unknown. In this work, we have characterized the expression of SNX9 and SNXl8 in LRO-containing cells and investigated their role in protein trafficking during the formation of LROs. Our results indicate that SNX9 and SNXl8 are not essential for the formation of LROs, nor for the sorting of melanosomal proteins. We investigated how the level of intracellular SNX9 protein is regulated and found that it is a substrate of the ubiquitin ligase Itch, a member of the NEDD4 family of E3 ubiquitin ligases. Itch ubiquitylates SNX9 and regulates SNX9 levels by enhancing its degradation. Using ? truncated proteins we found that the interaction with SNX9 is mediated by the proline-rich domain of Itch, a domain distinct from the conventional WW recognition domain, and the SH3 domain of SNX9. Interaction with the PRD of Itch is essential for SNX9 ubiquitylation and degradation. We further showed that Itch binding is not affected by tyrosine phosphorylation of SNX9. Using lentivector-mediated siRNA techniques, we found that Itch regulates the level of melanosomal proteins, while knock-down of SNX9 does not alter their level. Interestingly, we revealed that silencing of SNXIS affects the amount of the melanosomal protein Melan-A, but also of SNX9, and that SNXl8 can interact with SNX9. Taken together, our results highlight that the pool of substrates of NEDD4 family E3 ligases extends to proteins containing SH3 domains and provide insight into the potential functions of SNXI8. Résumé : Les membres de la famille des Sorting Nexins (SNX) jouent des rôles importants dans le trafic intracellulaire de protéines et membranes. Il a été démontré que la protéine SNX9, qui génère les tubules membranaires, interagit avec plusieurs composants de la machinerie d'endocytose et participe à l'endocytose des récepteurs de surface mediée par la clathrine. Aucune étude n'a investigué si SNX9 pourrait aussi participer à d'autres voies de trafic de protéines tel que le transport vésiculaire, et plus particulièrement la biogenèse des organites lysosomaux ("lysosome-related organelles", LR©s). SNXl8 est similaire à SNX9, mais sa fonction est largement inconnue. Dans ce travail, nous avons caractérisé l'expression de SNX9 et SNX18 dans des cellules contenants des LROs et investigué leur rôle dans le trafic de protéines pendant la formation des LROS. Nos résultats indiquent que SNX9 et SNXI8 ne sont essentiels ni pour la formation des LR©s, ni pour le trafic de protéines mélanosomales. Nous avons examiné la régulation du niveau intracellulaire de la protéine SNX9 et avons trouvé qu'elle est un substrat de l'ubiquitine ligase Itch, un membre de la famille NEDD4 des ubiquitine ligases E3. Itch ubiquitine SNX9 et régule les niveaux de SNX9 en augmentant sa dégradation. En utilisant des protéines mutées nous avons découvert que l'interaction avec SNX9 est médiée par le domaine riche en proline de Itch, qui est différent du domaine conventionnel de reconnaissance WW, et par le domaine SH3 de SNX9. L'interaction avec le domaine riche en proline de Itch est essentielle pour l'ubiquitination et la dégradation de SNX9. De plus, nous avons montré que cette liaison n'est pas affectée par la phosphorylation des résidus tyrosine de SNX9. En utilisant des vecteurs lentiviraux exprimant des siARN, nous avons trouvé que Itch régule les niveaux de protéines mélanosomales, alors que l'extinction de l'expression de SNX9 ne change pas leurs niveaux. En autre, nous avons révélé que la diminution de SNXl8 affecte le niveau de la protéine mélanosomale Melan-A et de SNX9, et aussi que SNXl8 peut interagir avec SNX9. En résumé, nos résultats démontrent que l'ensemble des substrats de la famille NEDD4 des ubiquitine ligases E3 s'élargit aux protéines contenant des domaines SH3 et ouvrent des perspectives sur les fonctions potentielles de SNXl8.

Pattern and clinical significance of cancer-testis gene expression in head and neck squamous cell carcinoma

Les cancer-testis antigènes appartiennent à la famille des antigènes tumoraux spécifiques. Ils ont montré un pouvoir immunogène chez les patients porteurs de différents cancers. En effet, ils stimulent sélectivement les lymphocytes cytotoxiques, et leur expression spécifique dans les tissus tumoraux en fait une cible idéale pour une vaccination antitumorale. Le but de cette étude est d'identifier l'expression de certains de ces antigènes, d'analyser leur valeur pronostique et de déterminer la meilleure cible antigénique pour permettre une immunothérapie spécifique dans les carcinomes épidermoïdes des voies aérodigestives supérieures. Le profil et le taux d'expression de 12 cancer-testis antigènes (MAGE-A1, MAGE-A3, MAGE-A4, MAGEA10, MAGE-C2, NY-ESO-1, LAGE-1, SSX-2, SSX-4, BAGE, GAGE-1/2, GAGE-3/4) et de 3 autres antigènes tumoraux spécifiques (PRAME, HERV-K-MEL, NA-17A) ont été évalués par RT-PCR sur 57 échantillons de cancers ORL primaires. Les paramètres tumoraux et cliniques ont été prospectivement collectés afin de corréler ces données avec le résultat de nos investigations immunobiologiques. Quatre-vingt-huit pour cent des tumeurs expriment au moins 1 antigène. Une co-expression de 3 gènes ou plus est détectée chez 59% des patients. MAGE-A4 (60%), MAGE-A3 (51%), PRAME (49%) et HERV-K-MEL (42%) sont les gènes le plus fréquemment exprimés. Ils sont totalement absents des muqueuses saines avoisinantes. La présence de MAGE-A et NY-ESO-1 à la surface des cellules a été vérifiée par immunohistochimie. Nos analyses statistiques ont permis d'identifier une diminution de la survie liée au cancer chez les patients porteurs d'une tumeur exprimant de multiples cancer-testis antigènes et notamment MAGE-A4 dont l'expression indépendante d'autres éléments cliniques s'associe statistiquement à un taux de survie diminué. Nos résultats ont permis d'identifier un rôle pronostique de l'expression des gènes associés aux tumeurs dont l'expression est apparemment liée à un phénotype de malignité plus élevé. Cette constatation, corroborée par l'identification parallèle d'un infiltrat lymphocytaire spécifique confirme l'utilité potentielle de certains cancer-testis antigènes comme cible pour une immunothérapie ciblée dans les carcinomes des voies aérodigestives supérieures

A gene expression signature that distinguishes desmoid tumours from nodular fasciitis.

Nodular fasciitis (NF) is a rapidly growing cellular mass composed of fibroblasts/myofibroblasts, usually localized in subcutaneous tissues, that typically undergoes fibrosis and almost never recurs. Desmoid tumours (DTs) are rare forms of fibroblastic/myofibroblastic growth that arise in deep soft tissues, display a propensity for local infiltration and recurrence, but fail to metastasize. Given that both entities are primarily fibroblastic/myofibroblastic lesions with overlapping histological features, their gene expression profiles were compared to identify differentially expressed genes that may provide not only potential diagnostic markers, but also clues as to the pathogenesis of each disorder. Differentially expressed transcripts (89 clones displaying increased expression in DTs and 246 clones displaying increased expression in NF) included genes encoding several receptor and non-receptor tyrosine kinases (EPHB3, PTPRF, GNAZ, SYK, LYN, EPHA4, BIRC3), transcription factors (TWIST1, PITX2, EYA2, OAS1, MITF, TCF20), and members of the Wnt signalling pathway (AXIN2, WISP1, SFRP). Remarkably, almost one-quarter of the differentially expressed genes encode proteins associated with inflammation and tissue remodelling, including members of the interferon (IFN), tumour necrosis factor (TNF), and transforming growth factor beta (TGF-beta) signalling pathways as well as metalloproteinases (MMP1, 9, 13, 23), urokinase plasminogen activator (PLAU), and cathepsins. The observations provide the first comparative molecular characterization of desmoid tumours and nodular fasciitis and suggest that selected tyrosine kinases, transcription factors, and members of the Wnt, TGF-beta, IFN, and TNF signalling pathways may be implicated in influencing and distinguishing their fate.

Early apoptosis of rod photoreceptors in Rpe65(-/-) mice is associated with the upregulated expression of lysosomal-mediated autophagic genes.

RPE65-related Leber's congenital amaurosis (LCA) is a rod-cone dystrophy whose clinical outcome is mainly attributed to the loss of rod photoreceptors followed by cone degeneration. Pathogenesis in Rpe65(-/-) mice is characterized by a slow and progressive degeneration of rods dependent on the constitutive activation of unliganded opsin. We previously reported that this opsin-mediated apoptosis of rods was dependent on Bcl-2-apoptotic pathway and Bax-induced pro-death activity. In this study, we report early initial apoptosis in the newly differentiated retina of Rpe65(-/-) mice. Apoptotic photoreceptors were identified as rods and resulted from pathological phototransduction signaling. This wave of early apoptosis triggered Bcl-2-related pathway and Bax apoptotic activity, while activation of the caspases was not induced. Following cellular stress, multiple signaling pathways are initiated which either commit cells to death or trigger pro-survival responses including autophagy. We report that Bcl-2-related early rod apoptosis was associated with the upregulation of autophagy markers including chaperone-mediated autophagy (CMA) substrate receptor LAMP-2 and lysosomal hydrolases Cathepsin S and Lysozyme. This suggests that lysosomal-mediated autophagy may be triggered in response to early rod apoptosis in Rpe65-LCA disease. These results highlight that Rpe65-related primary stress induces early signaling events, which trigger Bax-induced-apoptotic pathway and autophagy-mediated cellular response. These events may determine retinal cell fate, progression and severity of the disease.

Tumor cell budding and laminin-5 expression in colorectal carcinoma can be modulated by the tissue micro-environment.

Expression of laminin-5 alpha3, beta3 and gamma2 protein subunits was investigated in colorectal adenocarcinomas using immunostaining and confocal microscopy. The laminin-5 heterotrimer was found in basement membranes and as extracellular deposits in tumor stroma. In contrast to the alpha3 subunit, which was under-expressed, the gamma2 and beta3 subunits were detected in the cytoplasm of carcinoma cells dissociating (budding) from neoplastic tubules, suggestive of focal alterations in laminin-5 assembly and secretion. Laminin-5 gamma2 or beta3 subunit-reactive budding carcinoma cells expressed cytokeratins but not vimentin; they did not proliferate and were not apoptotic. Furthermore, expression of laminin-5 gamma2 and beta3 subunits in budding cells was associated with focal under-expression of the E-cadherin-beta-catenin complex. Results from xenograft experiments showed that budding activity in colorectal adenocarcinomas could be suppressed when these tumors grew at ectopic s.c. sites in nude mice. In vitro, cultured colon carcinoma cells, but not adenoma-derived tumor cells, shared the laminin-5 phenotype expressed by carcinoma cells in vivo. Using colon carcinoma cell lines implanted orthotopically and invading the cecum of nude mice, the laminin-5-associated budding was restored, indicating that this phenotype is not only determined by tumor cell properties but also dependent on the tissue micro-environment. Our results indicate that both laminin-5 alpha3 subunit expression and cell-cell cohesiveness are altered in budding carcinoma cells, which we consider to be actively invading. We propose that the local tissue micro-environment contributes to these events.

Homeodomain binding motifs modulate the probability of odorant receptor gene choice in transgenic mice.

Odorant receptor (OR) genes constitute with 1200 members the largest gene family in the mouse genome. A mature olfactory sensory neuron (OSN) is thought to express just one OR gene, and from one allele. The cell bodies of OSNs that express a given OR gene display a mosaic pattern within a particular region of the main olfactory epithelium. The mechanisms and cis-acting DNA elements that regulate the expression of one OR gene per OSN - OR gene choice - remain poorly understood. Here, we describe a reporter assay to identify minimal promoters for OR genes in transgenic mice, which are produced by the conventional method of pronuclear injection of DNA. The promoter transgenes are devoid of an OR coding sequence, and instead drive expression of the axonal marker tau-β-galactosidase. For four mouse OR genes (M71, M72, MOR23, and P3) and one human OR gene (hM72), a mosaic, OSN-specific pattern of reporter expression can be obtained in transgenic mice with contiguous DNA segments of only ~300 bp that are centered around the transcription start site (TSS). The ~150bp region upstream of the TSS contains three conserved sequence motifs, including homeodomain (HD) binding sites. Such HD binding sites are also present in the H and P elements, DNA sequences that are known to strongly influence OR gene expression. When a 19mer encompassing a HD binding site from the P element is multimerized nine times and added upstream of a MOR23 minigene that contains the MOR23 coding region, we observe a dramatic increase in the number of transgene-expressing founders and lines and in the number of labeled OSNs. By contrast, a nine times multimerized 19mer with a mutant HD binding site does not have these effects. We hypothesize that HD binding sites in the H and P elements and in OR promoters modulate the probability of OR gene choice.

Single cell gene expression profiling reveals qualitatively distinct CD8 T cells elicited by different gene-based vaccines

The CD8 T cell response generatedby gene-based vaccines is importantfor protective immunity againstmany infectious diseases but its complexityis incompletely understood.Here, we report that different vaccinesencoding HIV Env elicit qualitativelydistinct CD8 T cells that wereidentified by patterns of gene expressionin individual cells. Three alternativeprime-boost vector combinationsstimulated antigen-specific CD8 Tcell populations of similar magnitudeand function by intracellular cytokinestaining; however, single cell geneexpression profiling enabled the discriminationof distinct CM and EMCD8 cells elicited by the three vaccines.Two previously unrecognizedCD8 T cell subsets have been definedby their coexpression of Eomes,Cxcr3 and Ccr7; or Klrk1, Klrg1 andCcr5 in CM and EM cells respectively.

Retinoid X receptor and peroxisome proliferator-activated receptor activate an estrogen responsive gene independent of the estrogen receptor.

Estrogen receptors regulate transcription of genes essential for sexual development and reproductive function. Since the retinoid X receptor (RXR) is able to modulate estrogen responsive genes and both 9-cis RA and fatty acids influenced development of estrogen responsive tumors, we hypothesized that estrogen responsive genes might be modulated by RXR and the fatty acid receptor (peroxisome proliferator-activated receptor, PPAR). To test this hypothesis, transfection assays in CV-1 cells were performed with an estrogen response element (ERE) coupled to a luciferase reporter construct. Addition of expression vectors for RXR and PPAR resulted in an 11-fold increase in luciferase activity in the presence of 9-cis RA. Furthermore, mobility shift assays demonstrated binding of RXR and PPAR to the vitellogenin A2-ERE and an ERE in the oxytocin promoter. Methylation interference assays demonstrated that specific guanine residues required for RXR/PPAR binding to the ERE were similar to residues required for ER binding. Moreover, RXR domain-deleted constructs in transfection assays showed that activation required RXR since an RXR delta AF-2 mutant completely abrogated reporter activity. Oligoprecipitation binding studies with biotinylated ERE and (35)S-labeled in vitro translated RXR constructs confirmed binding of delta AF-2 RXR mutant to the ERE in the presence of baculovirus-expressed PPAR. Finally, in situ hybridization confirmed RXR and PPAR mRNA expression in estrogen responsive tissues. Collectively, these data suggest that RXR and PPAR are present in reproductive tissues, are capable of activating estrogen responsive genes and suggest that the mechanism of activation may involve direct binding of the receptors to estrogen response elements.

Evolution under monogamy feminizes gene expression in Drosophila melanogaster.

Many genes have evolved sexually dimorphic expression as a consequence of divergent selection on males and females. However, because the sexes share a genome, the extent to which evolution can shape gene expression independently in each sex is controversial. Here, we use experimental evolution to reveal suboptimal sex-specific expression for much of the genome. By enforcing a monogamous mating system in populations of Drosophila melanogaster for over 100 generations, we eliminated major components of selection on males: female choice and male-male competition. If gene expression is subject to sexually antagonistic selection, relaxed selection on males should cause evolution towards female optima. Monogamous males and females show this pattern of feminization in both the whole-body and head transcriptomes. Genes with male-biased expression patterns evolved decreased expression under monogamy, while genes with female-biased expression evolved increased expression, relative to polygamous populations. Our results demonstrate persistent and widespread evolutionary tension between male and female adaptation.

Gene expression profiles in rat mesenteric lymph nodes upon supplementation with Conjugated Linoleic Acid during gestation and suckling

Background Diet plays a role on the development of the immune system, and polyunsaturated fatty acids can modulate the expression of a variety of genes. Human milk contains conjugated linoleic acid (CLA), a fatty acid that seems to contribute to immune development. Indeed, recent studies carried out in our group in suckling animals have shown that the immune function is enhanced after feeding them with an 80:20 isomer mix composed of c9,t11 and t10,c12 CLA. However, little work has been done on the effects of CLA on gene expression, and even less regarding immune system development in early life. Results The expression profile of mesenteric lymph nodes from animals supplemented with CLA during gestation and suckling through dam's milk (Group A) or by oral gavage (Group B), supplemented just during suckling (Group C) and control animals (Group D) was determined with the aid of the specific GeneChip® Rat Genome 230 2.0 (Affymettrix). Bioinformatics analyses were performed using the GeneSpring GX software package v10.0.2 and lead to the identification of 89 genes differentially expressed in all three dietary approaches. Generation of a biological association network evidenced several genes, such as connective tissue growth factor (Ctgf), tissue inhibitor of metalloproteinase 1 (Timp1), galanin (Gal), synaptotagmin 1 (Syt1), growth factor receptor bound protein 2 (Grb2), actin gamma 2 (Actg2) and smooth muscle alpha actin (Acta2), as highly interconnected nodes of the resulting network. Gene underexpression was confirmed by Real-Time RT-PCR. Conclusions Ctgf, Timp1, Gal and Syt1, among others, are genes modulated by CLA supplementation that may have a role on mucosal immune responses in early life.

Mapping the transcription start points of the Staphylococcus aureus eap, emp, and vwb promoters reveals a conserved octanucleotide sequence that is essential for expression of these genes.

Mapping the transcription start points of the eap, emp, and vwb promoters revealed a conserved octanucleotide sequence (COS). Deleting this sequence abolished the expression of eap, emp, and vwb. However, electrophoretic mobility shift assays gave no evidence that this sequence was a binding site for SarA or SaeR, known regulators of eap and emp.

Interleukins (IL)-1 and IL-2 control IL-2 receptor alpha and beta expression in immature thymocytes.

Functional high-affinity interleukin-2 receptors (IL-2R) contain three transmembrane proteins, IL-2R alpha, beta and gamma. We have investigated the expression of IL-2R alpha and beta genes in immature mouse thymocytes. Previous work has shown that during differentiation these cells transiently express IL-2R alpha on their surface. Stimulation of IL-2R alpha+ and IL-2R alpha- immature thymocytes with phorbol 12-myristate 13-acetate and calcium ionophore induces synthesis of IL-2R alpha and IL-2R beta mRNA. Most of this response depends on autocrine stimulation by IL-2. IL-1 synergizes with IL-2 to induce a 120-fold increase in IL-2R alpha mRNA and a 14-fold increase in IL-2R beta mRNA levels. A large proportion of the stimulated cells contains both transcripts. These interleukins do not induce any differentiation to more mature phenotypes. Collectively, these results show that IL-2 plays a major role in the regulation of IL-2R expression in normal immature thymocyte. We suggest that this response to interleukins may be part of a homeostatic mechanism to increase the production of immature thymocytes during stress.

Semiquantitative RT-PCR measurement of gene expression in rat tissues including a correction for varying cell size and number

Background: Current methodology of gene expression analysis limits the possibilities of comparison between cells/tissues of organs in which cell size and/or number changes as a consequence of the study (e.g. starvation). A method relating the abundance of specific mRNA copies per cell may allow direct comparison or different organs and/or changing physiological conditions. Methods: With a number of selected genes, we analysed the relationship of the number of bases and the fluorescence recorded at a present level using cDNA standards. A lineal relationship was found between the final number of bases and the length of the transcript. The constants of this equation and those of the relationship between fluorescence and number of bases in cDNA were determined and a general equation linking the length of the transcript and the initial number of copies of mRNA was deduced for a given pre-established fluorescence setting. This allowed the calculation of the concentration of the corresponding mRNAs per g of tissue. The inclusion of tissue RNA and the DNA content per cell, allowed the calculation of the mRNA copies per cell. Results: The application of this procedure to six genes: Arbp, cyclophilin, ChREBP, T4 deiodinase 2, acetyl-CoA carboxylase 1 and IRS-1, in liver and retroperitoneal adipose tissue of food-restricted rats allowed precise measures of their changes irrespective of the shrinking of the tissue, the loss of cells or changes in cell size, factors that deeply complicate the comparison between changing tissue conditions. The percentage results obtained with the present methods were essentially the same obtained with the delta-delta procedure and with individual cDNA standard curve quantitative RT-PCR estimation. Conclusion: The method presented allows the comparison (i.e. as copies of mRNA per cell) between different genes and tissues, establishing the degree of abundance of the different molecular species tested.

Regional and cellular gene expression changes in human Huntington's disease brain.

Huntington's disease (HD) pathology is well understood at a histological level but a comprehensive molecular analysis of the effect of the disease in the human brain has not previously been available. To elucidate the molecular phenotype of HD on a genome-wide scale, we compared mRNA profiles from 44 human HD brains with those from 36 unaffected controls using microarray analysis. Four brain regions were analyzed: caudate nucleus, cerebellum, prefrontal association cortex [Brodmann's area 9 (BA9)] and motor cortex [Brodmann's area 4 (BA4)]. The greatest number and magnitude of differentially expressed mRNAs were detected in the caudate nucleus, followed by motor cortex, then cerebellum. Thus, the molecular phenotype of HD generally parallels established neuropathology. Surprisingly, no mRNA changes were detected in prefrontal association cortex, thereby revealing subtleties of pathology not previously disclosed by histological methods. To establish that the observed changes were not simply the result of cell loss, we examined mRNA levels in laser-capture microdissected neurons from Grade 1 HD caudate compared to control. These analyses confirmed changes in expression seen in tissue homogenates; we thus conclude that mRNA changes are not attributable to cell loss alone. These data from bona fide HD brains comprise an important reference for hypotheses related to HD and other neurodegenerative diseases.

Analysis of gene expression patterns in animals

During my PhD, my aim was to provide new tools to increase our capacity to analyse gene expression patterns, and to study on a large-scale basis the evolution of gene expression in animals. Gene expression patterns (when and where a gene is expressed) are a key feature in understanding gene function, notably in development. It appears clear now that the evolution of developmental processes and of phenotypes is shaped both by evolution at the coding sequence level, and at the gene expression level.Studying gene expression evolution in animals, with complex expression patterns over tissues and developmental time, is still challenging. No tools are available to routinely compare expression patterns between different species, with precision, and on a large-scale basis. Studies on gene expression evolution are therefore performed only on small genes datasets, or using imprecise descriptions of expression patterns.The aim of my PhD was thus to develop and use novel bioinformatics resources, to study the evolution of gene expression. To this end, I developed the database Bgee (Base for Gene Expression Evolution). The approach of Bgee is to transform heterogeneous expression data (ESTs, microarrays, and in-situ hybridizations) into present/absent calls, and to annotate them to standard representations of anatomy and development of different species (anatomical ontologies). An extensive mapping between anatomies of species is then developed based on hypothesis of homology. These precise annotations to anatomies, and this extensive mapping between species, are the major assets of Bgee, and have required the involvement of many co-workers over the years. My main personal contribution is the development and the management of both the Bgee database and the web-application.Bgee is now on its ninth release, and includes an important gene expression dataset for 5 species (human, mouse, drosophila, zebrafish, Xenopus), with the most data from mouse, human and zebrafish. Using these three species, I have conducted an analysis of gene expression evolution after duplication in vertebrates.Gene duplication is thought to be a major source of novelty in evolution, and to participate to speciation. It has been suggested that the evolution of gene expression patterns might participate in the retention of duplicate genes. I performed a large-scale comparison of expression patterns of hundreds of duplicated genes to their singleton ortholog in an outgroup, including both small and large-scale duplicates, in three vertebrate species (human, mouse and zebrafish), and using highly accurate descriptions of expression patterns. My results showed unexpectedly high rates of de novo acquisition of expression domains after duplication (neofunctionalization), at least as high or higher than rates of partitioning of expression domains (subfunctionalization). I found differences in the evolution of expression of small- and large-scale duplicates, with small-scale duplicates more prone to neofunctionalization. Duplicates with neofunctionalization seemed to evolve under more relaxed selective pressure on the coding sequence. Finally, even with abundant and precise expression data, the majority fate I recovered was neither neo- nor subfunctionalization of expression domains, suggesting a major role for other mechanisms in duplicate gene retention.