Differences in glucose transporter gene expression between rat pancreatic alpha- and beta-cells are correlated to differences in glucose transport but not in glucose utilization.

Glucose exerts inverse effects upon the secretory function of islet alpha- and beta-cells, suppressing glucagon release and increasing insulin release. This diverse action may result from differences in glucose transport and metabolism between the two cell types. The present study compares glucose transport in rat alpha- and beta-cells. beta-Cells transcribed GLUT2 and, to a lesser extent, GLUT 1; alpha-cells contained GLUT1 but no GLUT2 mRNA. No other GLUT-like sequences were found among cDNAs from alpha- or beta-cells. Both cell types expressed 43-kDa GLUT1 protein which was enhanced by culture. The 62-kDa beta-cell GLUT2 protein was converted to a 58-kDa protein after trypsin treatment of the cells without detectable consequences upon glucose transport kinetics. In beta-cells, the rates of glucose transport were 10-fold higher than in alpha-cells. In both cell types, glucose uptake exceeded the rates of glucose utilization by a factor of 10 or more. Glycolytic flux, measured as D-[5(3)H]glucose utilization, was comparable in alpha- and beta-cells between 1 and 10 mmol/liter substrate. In conclusion, differences in glucose transporter gene expression between alpha- and beta-cells can be correlated with differences in glucose transport kinetics but not with different glucose utilization rates.

Gene expression in cultured endometrium from women with different outcomes following IVF.

Estradiol and progesterone are crucial for the acquisition of receptivity and the change in transcriptional activity of target genes in the implantation window. The aim of this study was to differentiate the regulation of genes in the endometrium of patients with recurrent implantation failure (IF) versus those who became pregnant after in vitro fertilization (IVF) treatment. Moreover, the effect of embryo-derived factors on endometrial transcriptional activity was studied. Nine women with known IVF outcome (IF, M, miscarriage, OP, ongoing pregnancy) and undergoing hysteroscopy with endometrial biopsy were enrolled. Biopsies were taken during the midluteal phase. After culture in the presence of embryo-conditioned IVF media, total RNA was extracted and submitted to reverse transcription, target cDNA synthesis, biotin labelling, fragmentation and hybridization using the Affymetrix Human Genome U133A 2.0 Chip. Differential expression of selected genes was re-analysed by quantitative PCR, in which the results were calculated as threshold cycle differences between the groups and normalized to Glyceraldehyde phosphate dehydrogenase and beta-actin. Differences were seen for several genes from endometrial tissue between the IF and the pregnancy groups, and when comparing OP with M, 1875 up- and 1807 down-regulated genes were returned. Real-time PCR analysis confirmed up-regulation for somatostatin, PLAP-2, mucin 4 and CD163, and down-regulation of glycodelin, IL-24, CD69, leukaemia inhibitory factor and prolactin receptor between Op and M. When the different embryo-conditioned media were compared, no significant differential regulation could be demonstrated. Although microarray profiling may currently not be sensitive enough for studying the effects of embryo-derived factors on the endometrium, the observed differences in gene expression between M and OP suggest that it will become an interesting tool for the identification of fertility-relevant markers produced by the endometrium.

Impact of a cis-associated gene expression SNP on chromosome 20q11.22 on bipolar disorder susceptibility, hippocampal structure and cognitive performance.

BackgroundBipolar disorder is a highly heritable polygenic disorder. Recent enrichment analyses suggest that there may be true risk variants for bipolar disorder in the expression quantitative trait loci (eQTL) in the brain.AimsWe sought to assess the impact of eQTL variants on bipolar disorder risk by combining data from both bipolar disorder genome-wide association studies (GWAS) and brain eQTL.MethodTo detect single nucleotide polymorphisms (SNPs) that influence expression levels of genes associated with bipolar disorder, we jointly analysed data from a bipolar disorder GWAS (7481 cases and 9250 controls) and a genome-wide brain (cortical) eQTL (193 healthy controls) using a Bayesian statistical method, with independent follow-up replications. The identified risk SNP was then further tested for association with hippocampal volume (n = 5775) and cognitive performance (n = 342) among healthy individuals.ResultsIntegrative analysis revealed a significant association between a brain eQTL rs6088662 on chromosome 20q11.22 and bipolar disorder (log Bayes factor = 5.48; bipolar disorder P = 5.85×10(-5)). Follow-up studies across multiple independent samples confirmed the association of the risk SNP (rs6088662) with gene expression and bipolar disorder susceptibility (P = 3.54×10(-8)). Further exploratory analysis revealed that rs6088662 is also associated with hippocampal volume and cognitive performance in healthy individuals.ConclusionsOur findings suggest that 20q11.22 is likely a risk region for bipolar disorder; they also highlight the informative value of integrating functional annotation of genetic variants for gene expression in advancing our understanding of the biological basis underlying complex disorders, such as bipolar disorder.

Identification of molecular subtypes and gene expression patterns of breast cancer analysing RNA-seq data

Breast cancer is the most common diagnosed cancer and the leading cause of cancer death among females worldwide. It is considered a highly heterogeneous disease and it must be classified into more homogeneous groups. Hence, the purpose of this study was to classify breast tumors based on variations in gene expression patterns derived from RNA sequencing by using different class discovery methods. 42 breast tumors paired-samples were sequenced by Illumine Genome Analyzer and the data was analyzed and prepared by TopHat2 and htseq-count. As reported previously, breast cancer could be grouped into five main groups known as basal epithelial-like group, HER2 group, normal breast-like group and two Luminal groups with a distinctive expression profile. Classifying breast tumor samples by using PAM50 method, the most common subtype was Luminal B and was significantly associated with ESR1 and ERBB2 high expression. Luminal A subtype had ESR1 and SLC39A6 significant high expression, whereas HER2 subtype had a high expression of ERBB2 and CNNE1 genes and low luminal epithelial gene expression. Basal-like and normal-like subtypes were associated with low expression of ESR1, PgR and HER2, and had significant high expression of cytokeratins 5 and 17. Our results were similar compared with TGCA breast cancer data results and with known studies related with breast cancer classification. Classifying breast tumors could add significant prognostic and predictive information to standard parameters, and moreover, identify marker genes for each subtype to find a better therapy for patients with breast cancer.

Gene expression profiling of human hibernating myocardium: Increased expression of B-type natriuretic peptide and proenkephalin in hypocontractile vs normally-contracting regions of the heart

A greater understanding of the molecular basis of hibernating myocardium may assist in identifying those patients who would most benefit from revascularization. Paired heart biopsies were taken from hypocontractile and normally-contracting myocardium (identified by cardiovascular magnetic resonance) from 6 patients with chronic stable angina scheduled for bypass grafting. Gene expression profiles of hypocontractile and normally-contracting samples were compared using Affymetrix microarrays. The data for patients with confirmed hibernating myocardium were analysed separately and a different, though overlapping, set (up to 380) of genes was identified which may constitute a molecular fingerprint for hibernating myocardium. The expression of B-type natriuretic peptide (BNP) was increased in hypocontractile relative to normally-contracting myocardium. The expression of BNP correlated most closely with the expression of proenkephalin and follistatin 3, which may constitute additional heart failure markers. Our data illustrate differential gene expression in hypocontractile and/hibernating myocardium relative to normally-contracting myocardium within individual human hearts. Changes in expression of these genes, including increased relative expression of natriuretic and other factors, may constitute a molecular signature for hypocontractile and/or hibernating myocardium.

Influence of sex on basal and dickkopf-1 regulated gene expression in the bovine morula

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Global gene expression profiling of oral cavity cancers suggests molecular heterogeneity within anatomic subsites

Abstract Background Oral squamous cell carcinoma (OSCC) is a frequent neoplasm, which is usually aggressive and has unpredictable biological behavior and unfavorable prognosis. The comprehension of the molecular basis of this variability should lead to the development of targeted therapies as well as to improvements in specificity and sensitivity of diagnosis. Results Samples of primary OSCCs and their corresponding surgical margins were obtained from male patients during surgery and their gene expression profiles were screened using whole-genome microarray technology. Hierarchical clustering and Principal Components Analysis were used for data visualization and One-way Analysis of Variance was used to identify differentially expressed genes. Samples clustered mostly according to disease subsite, suggesting molecular heterogeneity within tumor stages. In order to corroborate our results, two publicly available datasets of microarray experiments were assessed. We found significant molecular differences between OSCC anatomic subsites concerning groups of genes presently or potentially important for drug development, including mRNA processing, cytoskeleton organization and biogenesis, metabolic process, cell cycle and apoptosis. Conclusion Our results corroborate literature data on molecular heterogeneity of OSCCs. Differences between disease subsites and among samples belonging to the same TNM class highlight the importance of gene expression-based classification and challenge the development of targeted therapies.

Gene expression in cultured endometrium from women with different outcomes following IVF

Estradiol and progesterone are crucial for the acquisition of receptivity and the change in transcriptional activity of target genes in the implantation window. The aim of this study was to differentiate the regulation of genes in the endometrium of patients with recurrent implantation failure (IF) versus those who became pregnant after in vitro fertilization (IVF) treatment. Moreover, the effect of embryo-derived factors on endometrial transcriptional activity was studied. Nine women with known IVF outcome (IF, M, miscarriage, OP, ongoing pregnancy) and undergoing hysteroscopy with endometrial biopsy were enrolled. Biopsies were taken during the midluteal phase. After culture in the presence of embryo-conditioned IVF media, total RNA was extracted and submitted to reverse transcription, target cDNA synthesis, biotin labelling, fragmentation and hybridization using the Affymetrix Human Genome U133A 2.0 Chip. Differential expression of selected genes was re-analysed by quantitative PCR, in which the results were calculated as threshold cycle differences between the groups and normalized to Glyceraldehyde phosphate dehydrogenase and beta-actin. Differences were seen for several genes from endometrial tissue between the IF and the pregnancy groups, and when comparing OP with M, 1875 up- and 1807 down-regulated genes were returned. Real-time PCR analysis confirmed up-regulation for somatostatin, PLAP-2, mucin 4 and CD163, and down-regulation of glycodelin, IL-24, CD69, leukaemia inhibitory factor and prolactin receptor between Op and M. When the different embryo-conditioned media were compared, no significant differential regulation could be demonstrated. Although microarray profiling may currently not be sensitive enough for studying the effects of embryo-derived factors on the endometrium, the observed differences in gene expression between M and OP suggest that it will become an interesting tool for the identification of fertility-relevant markers produced by the endometrium.

An atlas of human gene expression from massively parallel signature sequencing (MPSS).

We have used massively parallel signature sequencing (MPSS) to sample the transcriptomes of 32 normal human tissues to an unprecedented depth, thus documenting the patterns of expression of almost 20,000 genes with high sensitivity and specificity. The data confirm the widely held belief that differences in gene expression between cell and tissue types are largely determined by transcripts derived from a limited number of tissue-specific genes, rather than by combinations of more promiscuously expressed genes. Expression of a little more than half of all known human genes seems to account for both the common requirements and the specific functions of the tissues sampled. A classification of tissues based on patterns of gene expression largely reproduces classifications based on anatomical and biochemical properties. The unbiased sampling of the human transcriptome achieved by MPSS supports the idea that most human genes have been mapped, if not functionally characterized. This data set should prove useful for the identification of tissue-specific genes, for the study of global changes induced by pathological conditions, and for the definition of a minimal set of genes necessary for basic cell maintenance. The data are available on the Web at http://mpss.licr.org and http://sgb.lynxgen.com.

Novel loci for adiponectin levels and their influence on type 2 diabetes and metabolic traits: a multi-ethnic meta-analysis of 45,891 individuals.

Circulating levels of adiponectin, a hormone produced predominantly by adipocytes, are highly heritable and are inversely associated with type 2 diabetes mellitus (T2D) and other metabolic traits. We conducted a meta-analysis of genome-wide association studies in 39,883 individuals of European ancestry to identify genes associated with metabolic disease. We identified 8 novel loci associated with adiponectin levels and confirmed 2 previously reported loci (P = 4.5×10(-8)-1.2×10(-43)). Using a novel method to combine data across ethnicities (N = 4,232 African Americans, N = 1,776 Asians, and N = 29,347 Europeans), we identified two additional novel loci. Expression analyses of 436 human adipocyte samples revealed that mRNA levels of 18 genes at candidate regions were associated with adiponectin concentrations after accounting for multiple testing (p<3×10(-4)). We next developed a multi-SNP genotypic risk score to test the association of adiponectin decreasing risk alleles on metabolic traits and diseases using consortia-level meta-analytic data. This risk score was associated with increased risk of T2D (p = 4.3×10(-3), n = 22,044), increased triglycerides (p = 2.6×10(-14), n = 93,440), increased waist-to-hip ratio (p = 1.8×10(-5), n = 77,167), increased glucose two hours post oral glucose tolerance testing (p = 4.4×10(-3), n = 15,234), increased fasting insulin (p = 0.015, n = 48,238), but with lower in HDL-cholesterol concentrations (p = 4.5×10(-13), n = 96,748) and decreased BMI (p = 1.4×10(-4), n = 121,335). These findings identify novel genetic determinants of adiponectin levels, which, taken together, influence risk of T2D and markers of insulin resistance.

The MRC1/CD68 Ratio Is Positively Associated with Adipose Tissue Lipogenesis and with Muscle Mitochondrial Gene Expression in Humans.

BACKGROUND Alternative macrophages (M2) express the cluster differentiation (CD) 206 (MCR1) at high levels. Decreased M2 in adipose tissue is known to be associated with obesity and inflammation-related metabolic disturbances. Here we aimed to investigate MCR1 relative to CD68 (total macrophages) gene expression in association with adipogenic and mitochondrial genes, which were measured in human visceral [VWAT, n = 147] and subcutaneous adipose tissue [SWAT, n = 76] and in rectus abdominis muscle (n = 23). The effects of surgery-induced weight loss were also longitudinally evaluated (n = 6). RESULTS MCR1 and CD68 gene expression levels were similar in VWAT and SWAT. A higher proportion of CD206 relative to total CD68 was present in subjects with less body fat and lower fasting glucose concentrations. The ratio MCR1/CD68was positively associated with IRS1gene expression and with the expression of lipogenic genes such as ACACA, FASN and THRSP, even after adjusting for BMI. The ratio MCR1/CD68 in SWAT increased significantly after the surgery-induced weight loss (+44.7%; p = 0.005) in parallel to the expression of adipogenic genes. In addition, SWAT MCR1/CD68ratio was significantly associated with muscle mitochondrial gene expression (PPARGC1A, TFAM and MT-CO3). AT CD206 was confirmed by immunohistochemistry to be specific of macrophages, especially abundant in crown-like structures. CONCLUSION A decreased ratio MCR1/CD68 is linked to adipose tissue and muscle mitochondrial dysfunction at least at the level of expression of adipogenic and mitochondrial genes.

Increased endometrial placenta growth factor (PLGF) gene expression in women with successful implantation.

Objective: To analyze the vascularization of the endometrium via hysteroscopy and to assess its correlation with angiogenic factor gene expression and embryo implantation rate.Design: Cross-sectional study.Setting: Public university hospital.Patient(s): Patients undergoing hysteroscopy for supposed infertility.Intervention(s): Endometrial quality assessment according to Sakumoto-Masamoto, performed in the early secretory phase of the cycle. Collection of an endometrial tissue biopsy.Main Outcome Measure(s): RNA extraction, reverse transcription, and determination of gene expression of angiogenesis- and implantation-relevant factors using quantitative polymerase chain reaction. Retrieval of pregnancy information from the medical records.Result(s): Good quantity/quality RNA with infertility history was obtained from 63 participating women. Those with a "good" endometrium and subsequent pregnancy showed increased gene expression for placenta growth factor when compared with patients with a "bad" endometrium and who did not succeed with pregnancy to date. Nonpregnant women with a "good" endometrium presented an intermediate result. No significant differences were observed for several other genes tested, but trends in the same direction were observed.Conclusion(s): This study demonstrates for the first time that endometrial PLGF expression corresponds to the hysteroscopic appearance of the endometrium, and therefore has potential as a clinically relevant prognosticator for infertility treatment success. (Fertil Steril (R) 2011;96:663-8. (C)2011 by American Society for Reproductive Medicine.)

Long-term cerebral plasticity-related gene expression : a study of the respective role of CBP and CRTC1 in CREB transcriptional activation

1.1 AbstractThe treatment of memory disorders and cognitive deficits in various forms of mental retardation may greatly benefit from a better understanding of the molecular and cellular mechanisms of memory formation. Different forms of memory have distinct molecular requirements.Short-term memory (STM) is thought to be mediated by covalent modifications of existing synaptic molecules, such as phosphorylation or dephosphorylation of enzymes, receptors or ion channels. In contrast, long-term memoiy (LTM) is thought to be mediated by growth of new synapses and restructuring of existing synapses. There is extensive evidence that changes in gene expression and de novo protein synthesis are key processes for LTM formation. In this context, the transcription factor CREB (cAMP-response element-binding protein) was shown to be crucial. Activation of CREB requires phosphorylation of a serine residue (Ser-133), and the subsequent recruitment of a coactivator called CREB-binding protein (CBP). Moreover, we have recently shown that another coactivator called CREB Regulated Transcription Coactivator 1 (CRTC1) functions as a calcium- and cAMP-sensitive coincidence detector in neurons, and is involved in hippocampal long-term synaptic plasticity. Given the importance of cAMP and calcium signaling for plasticity-related gene expression in neurons and in astrocytes, we sought to determine the respective involvement of the CREB coactivators CBP and CRTC1 in CREB-mediated transcription.We developed various strategies to selectively interfere with these CREB coactivators in mouse primary neurons and in astrocytes in vitro. However, despite several pieces of evidence implicating CBP and/or CRTC1 in the regulation of neuronal plasticity genes, we could not clearly determine the respective requirement of these coactivators for the activation of these genes. Nevertheless, we showed that calcineurin activity, which is important for CRTC1 nuclear translocation, is necessary for the expression of some CREB-regulated plasticity genes. We associated this phenomena to physiopathological conditions observed in Down's syndrome. In addition, we demonstrated that in astrocytes, noradrenaline stimulates CREB-target gene expression through β-adrenergic receptor activation, intracellular cAMP pathway activation, and CRTC-induced CREB transactivation.Defining the respective role of CREB and its coactivators CBP and CRTC1 in neuronal and astrocytic cultures in vitro sets the stage for future in vivo studies and for the possible development of new therapeutic strategies to improve the treatment of memoiy and cognitive disorders.1.2 RésuméUne meilleure connaissance des mécanismes moléculaires et cellulaires responsables de la formation de la mémoire pourrait grandement améliorer le traitement des troubles de la mémoire ainsi que des déficits cognitifs observés dans différentes formes de pathologies psychiatriques telles que le retard mental. Les différentes formes de mémoire dépendent de processus moléculaires différents.La mémoire à court terme (STM) semble prendre forme suite à des modifications covalentes de molécules synaptiques préexistantes, telles que la phosphorylation ou la déphosphorylation d'enzymes, de récepteurs ou de canaux ioniques. En revanche, la mémoire à long terme (LTM) semble être due à la génération de nouvelles synapses et à la restructuration des synapses existantes. De nombreuses études ont permis de démontrer que les changements dans l'expression des gènes et la synthèse de protéine de novo sont des processus clés pour la formation de la LTM. Dans ce contexte, le facteur de transcription CREB (cAMP-response element-binding protein) s'est avéré être un élément crucial. L'activation de CREB nécessite la phosphorylation d'un résidu sérine (Ser-133), et le recrutement d'un coactivateur nommé CBP (CREB binding protein). En outre, nous avons récemment démontré qu'un autre coactivateur de CREB nommé CRTC1 (CREB Regulated Transcription Coactivator 1) agit comme un détecteur de coïncidence de l'AMP cyclique (AMPc) et du calcium dans les neurones et qu'il est impliqué dans la formation de la plasticité synaptique à long terme dans l'hippocampe. Etant donné l'importance des voies de l'AMPc et du calcium dans l'expression des gènes impliqués dans la plasticité cérébrale, nous voulions déterminer le rôle respectif des coactivateurs de CREB, CBP et CRTC1.Nous avons développé diverses stratégies pour interférer de façon sélective avec les coactivateurs de CREB dans les neurones et dans les astrocytes chez la souris in vitro. Nos résultats indiquent que CBP et CRTC1 sont tous deux impliqués dans la transcription dépendante de CREB induite par l'AMPc et le calcium dans les neurones. Cependant, malgré plusieurs évidences impliquant CBP et/ou CRTC1 dans l'expression de gènes de plasticité neuronale, nous n'avons pas pu déterminer clairement leur nécessité respective pour l'activation de ces gènes. Toutefois, nous avons montré que l'activité de la calcineurine, dont dépend la translocation nucléaire de CRTC1, est nécessaire à l'expression de certains de ces gènes. Nous avons pu associer ce phénomène à une condition physiopathologique observée dans le syndrome de Down. Nous avons également montré que dans les astrocytes, la noradrénaline stimule l'expression de gènes cibles de CREB par une activation des récepteurs β- adrénergiques, l'activation de la voie de l'AMPc et la transactivation de CREB par les CRTCs.Définir le rôle respectif de CREB et de ses coactivateurs CBP et CRTC1 dans les neurones et dans les astrocytes in vitro permettra d'acquérir les connaissances nécessaires à de futures études in vivo et, à plus long terme d'éventuellement développer des stratégies thérapeutiques pour améliorer les traitements des troubles cognitifs.

p21 as a transcriptional co-repressor of S-phase and mitotic control genes.

It has been previously described that p21 functions not only as a CDK inhibitor but also as a transcriptional co-repressor in some systems. To investigate the roles of p21 in transcriptional control, we studied the gene expression changes in two human cell systems. Using a human leukemia cell line (K562) with inducible p21 expression and human primary keratinocytes with adenoviral-mediated p21 expression, we carried out microarray-based gene expression profiling. We found that p21 rapidly and strongly repressed the mRNA levels of a number of genes involved in cell cycle and mitosis. One of the most strongly down-regulated genes was CCNE2 (cyclin E2 gene). Mutational analysis in K562 cells showed that the N-terminal region of p21 is required for repression of gene expression of CCNE2 and other genes. Chromatin immunoprecipitation assays indicated that p21 was bound to human CCNE2 and other p21-repressed genes gene in the vicinity of the transcription start site. Moreover, p21 repressed human CCNE2 promoter-luciferase constructs in K562 cells. Bioinformatic analysis revealed that the CDE motif is present in most of the promoters of the p21-regulated genes. Altogether, the results suggest that p21 exerts a repressive effect on a relevant number of genes controlling S phase and mitosis. Thus, p21 activity as inhibitor of cell cycle progression would be mediated not only by the inhibition of CDKs but also by the transcriptional down-regulation of key genes.

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.