A fish-specific transposable element shapes the repertoire of p53 target genes in zebrafish.

Transposable elements, as major components of most eukaryotic organisms' genomes, define their structural organization and plasticity. They supply host genomes with functional elements, for example, binding sites of the pleiotropic master transcription factor p53 were identified in LINE1, Alu and LTR repeats in the human genome. Similarly, in this report we reveal the role of zebrafish (Danio rerio) EnSpmN6_DR non-autonomous DNA transposon in shaping the repertoire of the p53 target genes. The multiple copies of EnSpmN6_DR and their embedded p53 responsive elements drive in several instances p53-dependent transcriptional modulation of the adjacent gene, whose human orthologs were frequently previously annotated as p53 targets. These transposons define predominantly a set of target genes whose human orthologs contribute to neuronal morphogenesis, axonogenesis, synaptic transmission and the regulation of programmed cell death. Consistent with these biological functions the orthologs of the EnSpmN6_DR-colonized loci are enriched for genes expressed in the amygdala, the hippocampus and the brain cortex. Our data pinpoint a remarkable example of convergent evolution: the exaptation of lineage-specific transposons to shape p53-regulated neuronal morphogenesis-related pathways in both a hominid and a teleost fish.

Phylogenetic analyses of C4 photosynthesis evolution in grasses: genetics of convergence and convergence of genes

Summary : During the evolutionary diversification of organisms, similar ecological constraints led to the recurrent appearances of the same traits (phenotypes) in distant lineages, a phenomenon called convergence. In most cases, the genetic origins of the convergent traits remain unknown, but recent studies traced the convergent phenotypes to recurrent alterations of the same gene or, in a few cases, to identical genetic changes. However, these cases remain anecdotal and there is a need for a study system that evolved several times independently and whose genetic determinism is well resolved and straightforward, such as C4 photosynthesis. This adaptation to warm environments, possibly driven by past atmospheric CO2 decreases, consists in a CO2-concentrating pump, created by numerous morphological and biochemical novelties. All genes encoding C4 enzymes already existed in C3 ancestors, and are supposed to have been recruited through gene duplication followed by neo-functionalization, to acquire the cell specific expression pattern and altered kinetic properties that characterize Ca-specific enzymes. These predictions have so far been tested only in species-poor and ecologically marginal C4 dicots. The monocots, and especially the grass family (Poaceae), the most important C4 family in terms of species number, ecological dominance and economical importance, have been largely under-considered as suitable study systems. This thesis aimed at understanding the evolution of the C4 trait in grasses at a molecular level and to use the genetics of C4 photosynthesis to infer the evolutionary history of the C4 phenotype and its driving selective pressures. A molecular phylogeny of grasses and affiliated monocots identified 17 to 18 independent acquisitions of the C4 pathway in the grass family. A relaxed molecular clock was used to date these events and the first C4 evolution was estimated in the Chloridoideae subfamily, between 32-25 million years ago, at a period when atmospheric CO2 abruptly declined. Likelihood models showed that after the COZ decline the probability of evolving the C4 pathway strongly increased, confirming low CO2 as a likely driver of C4 photosynthesis evolution. In order to depict the genetic changes linked to the numerous C4 origins, genes encoding phopshoenolpyruvate carboxylase (PEPC), the key-enzyme responsible for the initial fixation of atmospheric CO2 in the C4 pathway, were isolated from a large sample of C3 and C4 grasses. Phylogenetic analyses were used to reconstruct the evolutionary history of the PEPC multigene family and showed that the evolution of C4-specific PEPC had been driven by positive selection on 21 codons simultaneously in up to eight C4 lineages. These selective pressures led to numerous convergent genetic changes in many different C4 clades, highlighting the repeatability of some evolutionary processes, even at the molecular level. PEPC C4-adaptive changes were traced and used to show multiple appearances of the C, pathway in clades where species tree inferences were unable to differentiate multiple C4 appearances and a single appearance followed by C4 to C3 reversion. Further investigations of genes involved in some of the C4 subtypes only (genes encoding decarboxylating enzymes NADP-malic enzyme and phosphoenolpyruvate carboxykinase) showed that these C4-enzymes also evolved through strong positive selection and underwent parallel genetic changes during the different Ca origins. The adaptive changes on these subtype-specific C4 genes were used to retrace the history of the C4-subtypes phenotypes, which revealed that the evolution of C4-PEPC and C4-decarboxylating enzymes was in several cases disconnected, emphasizing the multiplicity of the C4 trait and the gradual acquisition of the features that create the CO2-pump. Finally, phylogenetic analyses of a gene encoding the Rubisco (the enzyme responsible for the fixation of CO2 into organic compounds in all photosynthetic organisms) showed that C4 evolution switched the selective pressures on this gene. Five codons were recurrently mutated to adapt the enzyme kinetics to the high CO2 concentrations of C4 photosynthetic cells. This knowledge could be used to introgress C4-like Rubisco in C3 crops, which could lead to an increased yield under predicted future high CO2 atmosphere. Globally, the phylogenetic framework adopted during this thesis demonstrated the widespread occurrence of genetic convergence on C4-related enzymes. The genetic traces of C4 photosynthesis evolution allowed reconstructing events that happened during the last 30 million years and proved the usefulness of studying genes directly responsible for phenotype variations when inferring evolutionary history of a given trait. Résumé Durant la diversification évolutive des organismes, des pressions écologiques similaires ont amené à l'apparition récurrente de certains traits (phénotypes) dans des lignées distantes, un phénomène appelé évolution convergente. Dans la plupart des cas, l'origine génétique des traits convergents reste inconnue mais des études récentes ont montré qu'ils étaient dus dans certains cas à des changements répétés du même gène ou, dans de rares cas, à des changements génétiques identiques. Malgré tout, ces cas restent anecdotiques et il y a un réel besoin d'un système d'étude qui ait évolué indépendamment de nombreuses fois et dont le déterminisme génétique soit clairement identifié. La photosynthèse dite en Ça répond à ces critères. Cette adaptation aux environnements chauds, dont l'évolution a pu être encouragé par des baisses passées de la concentration atmosphérique en CO2, est constituée de nombreuses nouveautés morphologiques et biochimiques qui créent une pompe à CO2. La totalité des gènes codant les enzymes Ç4 étaient déjà présents dans les ancêtres C3. Leur recrutement pour la photosynthèse Ç4 est supposé s'être fait par le biais de duplications géniques suivies par une néo-fonctionnalisation pour leur conférer l'expression cellule-spécifique et les propriétés cinétiques qui caractérisent les enzymes C4. Ces prédictions n'ont jusqu'à présent été testées que dans des familles C4 contenant peu d'espèces et ayant un rôle écologique marginal. Les graminées (Poaceae), qui sont la famille C4 la plus importante, tant en termes de nombre d'espèces que de dominance écologique et d'importance économique, ont toujours été considérés comme un système d'étude peu adapté et ont fait le sujet de peu d'investigations évolutives. Le but de cette thèse était de comprendre l'évolution de la photosynthèse en C4 chez les graminées au niveau génétique et d'utiliser les gènes pour inférer l'évolution du phénotype C4 ainsi que les pressions de sélection responsables de son évolution. Une phylogénie moléculaire de la famille des graminées et des monocotylédones apparentés a identifié 17 à 18 acquisitions indépendantes de la photosynthèse chez les graminées. Grâce à une méthode d'horloge moléculaire relâchée, ces évènements ont été datés et la première apparition C4 a été estimée dans la sous-famille des Chloridoideae, il y a 32 à 25 millions d'années, à une période où les concentrations atmosphériques de CO2 ont décliné abruptement. Des modèles de maximum de vraisemblance ont montré qu'à la suite du déclin de CO2, la probabilité d'évoluer la photosynthèse C4 a fortement augmenté, confirmant ainsi qu'une faible concentration de CO2 est une cause potentielle de l'évolution de la photosynthèse C4. Afin d'identifier les mécanismes génétiques responsables des évolutions répétées de la photosynthèse C4, un segment des gènes codant pour la phosphoénolpyruvate carboxylase (PEPC), l'enzyme responsable de la fixation initiale du CO2 atmosphérique chez les plantes C4, ont été séquencés dans une centaine de graminées C3 et C4. Des analyses phylogénétiques ont permis de reconstituer l'histoire évolutive de la famille multigénique des PEPC et ont montré que l'évolution de PEPC spécifiques à la photosynthèse Ça a été causée par de la sélection positive agissant sur 21 codons, et ce simultanément dans huit lignées C4 différentes. Cette sélection positive a conduit à un grand nombre de changements génétiques convergents dans de nombreux clades différents, ce qui illustre la répétabilité de certains phénomènes évolutifs, et ce même au niveau génétique. Les changements sur la PEPC liés au C4 ont été utilisés pour confirmer des évolutions indépendantes du phénotype C4 dans des clades où l'arbre des espèces était incapable de différencier des apparitions indépendantes d'une seule apparition suivie par une réversion de C4 en C3. En considérant des gènes codant des protéines impliquées uniquement dans certains sous-types C4 (deux décarboxylases, l'enzyme malique à NADP et la phosphoénolpyruvate carboxykinase), des études ultérieures ont montré que ces enzymes C4 avaient elles-aussi évolué sous forte sélection positive et subi des changements génétiques parallèles lors des différentes origines de la photosynthèse C4. Les changements adaptatifs sur ces gènes liés seulement à certains sous-types C4 ont été utilisés pour retracer l'histoire des phénotypes de sous-types C4, ce qui a révélé que les caractères formant le trait C4 ont, dans certains cas, évolué de manière déconnectée. Ceci souligne la multiplicité du trait C4 et l'acquisition graduelle de composants participant à la pompe à CO2 qu'est la photosynthèse C4. Finalement, des analyses phylogénétiques des gènes codant pour la Rubisco (l'enzyme responsable de la fixation du CO2 en carbones organiques dans tous les organismes photosynthétiques) ont montré que l'évolution de la photosynthèse Ça a changé les pressions de sélection sur ce gène. Cinq codons ont été mutés de façon répétée afin d'adapter les propriétés cinétiques de la Rubisco aux fortes concentrations de CO2 présentes dans les cellules photosynthétiques des plantes C4. Globalement, l'approche phylogénétique adoptée durant cette thèse de doctorat a permis de démontré des phénomène fréquents de convergence génétique sur les enzymes liées à la photosynthèse C4. Les traces génétiques de l'évolution de la photosynthèse C4 ont permis de reconstituer des évènements qui se sont produits durant les derniers 30 millions d'années et ont prouvé l'utilité d'étudier des gènes directement responsables des variations phénotypiques pour inférer l'histoire évolutive d'un trait donné.

Analyses of six homologous proteins of Protochlamydia amoebophila UWE25 encoded by large GC-rich genes (lgr): a model of evolution and concatenation of leucine-rich repeats.

BACKGROUND: Along the chromosome of the obligate intracellular bacteria Protochlamydia amoebophila UWE25, we recently described a genomic island Pam100G. It contains a tra unit likely involved in conjugative DNA transfer and lgrE, a 5.6-kb gene similar to five others of P. amoebophila: lgrA to lgrD, lgrF. We describe here the structure, regulation and evolution of these proteins termed LGRs since encoded by "Large G+C-Rich" genes. RESULTS: No homologs to the whole protein sequence of LGRs were found in other organisms. Phylogenetic analyses suggest that serial duplications producing the six LGRs occurred relatively recently and nucleotide usage analyses show that lgrB, lgrE and lgrF were relocated on the chromosome. The C-terminal part of LGRs is homologous to Leucine-Rich Repeats domains (LRRs). Defined by a cumulative alignment score, the 5 to 18 concatenated octacosapeptidic (28-meric) LRRs of LGRs present all a predicted alpha-helix conformation. Their closest homologs are the 28-residue RI-like LRRs of mammalian NODs and the 24-meres of some Ralstonia and Legionella proteins. Interestingly, lgrE, which is present on Pam100G like the tra operon, exhibits Pfam domains related to DNA metabolism. CONCLUSION: Comparison of the LRRs, enable us to propose a parsimonious evolutionary scenario of these domains driven by adjacent concatenations of LRRs. Our model established on bacterial LRRs can be challenged in eucaryotic proteins carrying less conserved LRRs, such as NOD proteins and Toll-like receptors.

Gain of function mutations in CgPDR1 of Candida glabrata not only mediate antifungal resistance but also enhance virulence.

CgPdr1p is a Candida glabrata Zn(2)-Cys(6) transcription factor involved in the regulation of the ABC-transporter genes CgCDR1, CgCDR2, and CgSNQ2, which are mediators of azole resistance. Single-point mutations in CgPDR1 are known to increase the expression of at least CgCDR1 and CgCDR2 and thus to contribute to azole resistance of clinical isolates. In this study, we investigated the incidence of CgPDR1 mutations in a large collection of clinical isolates and tested their relevance, not only to azole resistance in vitro and in vivo, but also to virulence. The comparison of CgPDR1 alleles from azole-susceptible and azole-resistant matched isolates enabled the identification of 57 amino acid substitutions, each positioned in distinct CgPDR1 alleles. These substitutions, which could be grouped into three different "hot spots," were gain of function (GOF) mutations since they conferred hyperactivity to CgPdr1p revealed by constitutive high expression of ABC-transporter genes. Interestingly, the major transporters involved in azole resistance (CgCDR1, CgCDR2, and CgSNQ2) were not always coordinately expressed in presence of specific CgPDR1 GOF mutations, thus suggesting that these are rather trans-acting elements (GOF in CgPDR1) than cis-acting elements (promoters) that lead to azole resistance by upregulating specific combinations of ABC-transporter genes. Moreover, C. glabrata isolates complemented with CgPDR1 hyperactive alleles were not only more virulent in mice than those with wild type alleles, but they also gained fitness in the same animal model. The presence of CgPDR1 hyperactive alleles also contributed to fluconazole treatment failure in the mouse model. In conclusion, this study shows for the first time that CgPDR1 mutations are not only responsible for in vitro/in vivo azole resistance but that they can also confer a selective advantage under host conditions.

Meta-analysis of estrogen response in MCF-7 distinguishes early target genes involved in signaling and cell proliferation from later target genes involved in cell cycle and DNA repair.

ABSTRACT: BACKGROUND: Many studies have been published outlining the global effects of 17 beta-estradiol (E2) on gene expression in human epithelial breast cancer derived MCF-7 cells. These studies show large variation in results, reporting between ~100 and ~1500 genes regulated by E2, with poor overlap. RESULTS: We performed a meta-analysis of these expression studies, using the Rank product method to obtain a more accurate and stable list of the differentially expressed genes, and of pathways regulated by E2. We analyzed 9 time-series data sets, concentrating on response at 3-4 hrs (early) and at 24 hrs (late). We found >1000 statistically significant probe sets after correction for multiple testing at 3-4 hrs, and >2000 significant probe sets at 24 hrs. Differentially expressed genes were examined by pathway analysis. This revealed 15 early response pathways, mostly related to cell signaling and proliferation, and 20 late response pathways, mostly related to breast cancer, cell division, DNA repair and recombination. CONCLUSIONS: Our results show that meta-analysis identified more differentially expressed genes than the individual studies, and that these genes act together in networks. These results provide new insight into E2 regulated mechanisms, especially in the context of breast cancer.

Validated prediction of clinical outcome in sarcomas and multiple types of cancer on the basis of a gene expression signature related to genome complexity.

Sarcomas are heterogeneous and aggressive mesenchymal tumors. Histological grading has so far been the best predictor for metastasis-free survival, but it has several limitations, such as moderate reproducibility and poor prognostic value for some histological types. To improve patient grading, we performed genomic and expression profiling in a training set of 183 sarcomas and established a prognostic gene expression signature, complexity index in sarcomas (CINSARC), composed of 67 genes related to mitosis and chromosome management. In a multivariate analysis, CINSARC predicts metastasis outcome in the training set and in an independent 127 sarcomas validation set. It is superior to the Fédération Francaise des Centres de Lutte Contre le Cancer grading system in determining metastatic outcome for sarcoma patients. Furthermore, it also predicts outcome for gastrointestinal stromal tumors (GISTs), breast carcinomas and lymphomas. Application of the signature will permit more selective use of adjuvant therapies for people with sarcomas, leading to decreased iatrogenic morbidity and improved outcomes for such individuals.

Identification of CT46/HORMAD1, an immunogenic cancer/testis antigen encoding a putative meiosis-related protein.

Transcripts with ESTs derived exclusively or predominantly from testis, and not from other normal tissues, are likely to be products of genes with testis-restricted expression, and are thus potential cancer/testis (CT) antigen genes. A list of 371 genes with such characteristics was compiled by analyzing publicly available EST databases. RT-PCR analysis of normal and tumor tissues was performed to validate an initial selection of 20 of these genes. Several new CT and CT-like genes were identified. One of these, CT46/HORMAD1, is expressed strongly in testis and weakly in placenta; the highest level of expression in other tissues is <1% of testicular expression. The CT46/HORMAD1 gene was expressed in 31% (34/109) of the carcinomas examined, with 11% (12/109) showing expression levels >10% of the testicular level of expression. CT46/HORMAD1 is a single-copy gene on chromosome 1q21.3, encoding a putative protein of 394 aa. Conserved protein domain analysis identified a HORMA domain involved in chromatin binding. The CT46/HORMAD1 protein was found to be homologous to the prototype HORMA domain-containing protein, Hop1, a yeast meiosis-specific protein, as well as to asy1, a meiotic synaptic mutant protein in Arabidopsis thaliana.

Calcineurin A of Candida albicans: involvement in antifungal tolerance, cell morphogenesis and virulence.

The azole antifungal fluconazole possesses only fungistatic activity in Candida albicans and, therefore, this human pathogen is tolerant to this agent. However, tolerance to fluconazole can be inhibited when C. albicans is exposed to fluconazole combined with the immunosuppressive drug cyclosporin A, which is known to inhibit calcineurin activity in yeast. A mutant lacking both alleles of a gene encoding the calcineurin A subunit (CNA) lost viability in the presence of fluconazole, thus making calcineurin essential for fluconazole tolerance. Consistent with this observation, tolerance to fluconazole was modulated by calcium ions or by the expression of a calcineurin A derivative autoactivated by the removal of its C-terminal inhibitory domain. Interestingly, CNA was also essential for tolerance to other antifungal agents (voriconazole, itraconazole, terbinafine, amorolfine) and to several other metabolic inhibitors (caffeine, brefeldin A, mycophenolic acid, fluphenazine) or cell wall-perturbing agents (SDS, calcofluor white, Congo red), thus indicating that the calcineurin pathway plays an important role in the survival of C. albicans in the presence of external growth inhibitors. Several genes, including PMC1, a vacuolar calcium P-type ATPase, were regulated in a calcineurin- and fluconazole-dependent manner. However, PMC1 did not play a direct role in the survival of C. albicans when exposed to fluconazole. In addition to these different properties, calcineurin was found to affect colony morphology in several media known to modulate the C. albicans dimorphic switch. In particular, calcineurin was found to be essential for C. albicans viability in serum-containing media. Finally, calcineurin was found to be necessary for the virulence of C. albicans in a mice model of infection, thus making calcineurin an important element for adequate adaptation to the conditions of the host environment.

Vitellogenin B2 gene in Xenopus laevis: isolation, in vitro transcription and relation to other vitellogenin genes.

The isolation of the four Xenopus laevis vitellogenin genes has been completed by the purification from a DNA library of the B2 gene together with its flanking sequences. The overlapping DNA fragments analyzed cover 34 kilobases. The B2 gene which has a length of 17.5 kilobases was characterized by heteroduplex and R-loop mapping in the electron microscope and by in vitro transcription in a HeLa whole-cell extract. Its structural organization is compared with that of the closely related B1 gene. The mRNA-coding sequence of about 6 kilobases is interrupted 34 times in the B1 gene and 33 times in the B2 gene. Sequence homology between the two genes was not only found in exons. In addition, 54% of the intron sequences as well as 63% and 48.5% respectively of the 5' and 3' flanking sequences, show enough homology to form stable duplexes. These findings are compared with earlier results obtained with the two other closely related members of the vitellogenin gene family, the A1 and the A2 genes.

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.

Milbemycins: more than efflux inhibitors for fungal pathogens.

Existing antifungal agents are still confronted to activities limited to specific fungal species and to the development of resistance. Several improvements are possible either by tackling and overcoming resistance or exacerbating the activity of existing antifungal agents. In Candida glabrata, azole resistance is almost exclusively mediated by ABC transporters (including C. glabrata CDR1 [CgCDR1] and CgCDR2) via gain-of-function mutations in the transcriptional activator CgPDR1 or by mitochondrial dysfunctions. We also observed that azole resistance was correlating with increasing virulence and fitness of C. glabrata in animal models of infection. This observation motivated the re-exploitation of ABC transporter inhibitors as a possible therapeutic intervention to decrease not only the development of azole resistance but also to interfere with the virulence of C. glabrata. Milbemycins are known ABC transporter inhibitors, and here we used commercially available milbemycin A3/A4 oxim derivatives to verify this effect. As expected, the derivatives were inhibiting C. glabrata efflux with the highest activity for A3 oxim below 1 μg/ml. More surprising was that oxim derivatives had intrinsic fungicidal activity above 3.2 μg/ml, thus highlighting effects additional to the efflux inhibition. Similar values were obtained with C. albicans. Our data show that the fungicidal activity could be related to reactive oxygen species formation in these species. Transcriptional analysis performed both in C. glabrata and C. albicans exposed to A3 oxim highlighted a core of commonly regulated genes involved in stress responses, including genes involved in oxidoreductive processes, protein ubiquitination, and vesicle trafficking, as well as mitogen-activated protein kinases. However, the transcript profiles contained also species-specific signatures. Following these observations, experimental treatments of invasive infections were performed in mice treated with the commercial A3/A4 oxim preparation alone or in combination with fluconazole. Tissue burden analysis revealed that oxims on their own were able to decrease fungal burdens in both Candida species. In azole-resistant isolates, oxims acted synergistically in vivo with fluconazole to reduce fungal burden to levels of azole-susceptible isolates. In conclusion, we show here the potential of milbemycins not only as drug efflux inhibitors but also as effective fungal growth inhibitors in C. glabrata and C. albicans.

Evolution of mRNA expression levels of autosomal and sex-linked genes in amniotes

In addition to differences in protein-coding gene sequences, changes in expression resulting from mutations in regulatory sequences have long been hypothesized to be responsible for phenotypic differences between species. However, unlike comparison of genome sequences, few studies, generally restricted to pairwise comparisons of closely related mammalian species, have assessed between-species differences at the transcriptome level. They reported that gene expression evolves at different rates in various organs and in a pattern that is overall consistent with neutral models of evolution. In the first part of my thesis, I investigated the evolution of gene expression in therian mammals (i.e.7 placental and marsupials), based on microarray data from human, mouse and the gray short-tailed opossum (Monodelphis domestica). In addition to autosomal genes, a special focus was given to the evolution of X-linked genes. The therian X chromosome was recently shown to be younger than previously thought and to harbor a specific gene content (e.g., genes involved in brain or reproductive functions) that is thought to have been shaped by specific sex-related evolutionary forces. Sex chromosomes derive from ordinary autosomes and their differentiation led to the degeneration of the Y chromosome (in mammals) or W chromosome (in birds). Consequently, X- or Z-linked genes differ in gene dose between males and females such that the heterogametic sex has half the X/Z gene dose compared to the ancestral state. To cope with this dosage imbalance, mammals have been reported to have evolved mechanisms of dosage compensation.¦In the first project, I could first show that transcriptomes evolve at different rates in different organs. Out of the five tissues I investigated, the testis is the most rapidly evolving organ at the gene expression level while the brain has the most conserved transcriptome. Second, my analyses revealed that mammalian gene expression evolution is compatible with a neutral model, where the rates of change in gene expression levels is linked to the efficiency of purifying selection in a given lineage, which, in turn, is determined by the long-term effective population size in that lineage. Thus, the rate of DNA sequence evolution, which could be expected to determine the rate of regulatory sequence change, does not seem to be a major determinant of the rate of gene expression evolution. Thus, most gene expression changes seem to be (slightly) deleterious. Finally, X-linked genes seem to have experienced elevated rates of gene expression change during the early stage of X evolution. To further investigate the evolution of mammalian gene expression, we generated an extensive RNA-Seq gene expression dataset for nine mammalian species and a bird. The analyses of this dataset confirmed the patterns previously observed with microarrays and helped to significantly deepen our view on gene expression evolution.¦In a specific project based on these data, I sought to assess in detail patterns of evolution of dosage compensation in amniotes. My analyses revealed the absence of male to female dosage compensation in monotremes and its presence in marsupials and, in addition, confirmed patterns previously described for placental mammals and birds. I then assessed the global level of expression of X/Z chromosomes and contrasted this with its ancestral gene expression levels estimated from orthologous autosomal genes in species with non-homologous sex chromosomes. This analysis revealed a lack of up-regulation for placental mammals, the level of expression of X-linked genes being proportional to gene dose. Interestingly, the ancestral gene expression level was at least partially restored in marsupials as well as in the heterogametic sex of monotremes and birds. Finally, I investigated alternative mechanisms of dosage compensation and found that gene duplication did not seem to be a widespread mechanism to restore the ancestral gene dose. However, I could show that placental mammals have preferentially down-regulated autosomal genes interacting with X-linked genes which underwent gene expression decrease, and thus identified a novel alternative mechanism of dosage compensation.

Bioinformatics-driven identification and examination of candidate genes for non-alcoholic fatty liver disease.

ObjectiveCandidate genes for non-alcoholic fatty liver disease (NAFLD) identified by a bioinformatics approach were examined for variant associations to quantitative traits of NAFLD-related phenotypes.Research Design and MethodsBy integrating public database text mining, trans-organism protein-protein interaction transferal, and information on liver protein expression a protein-protein interaction network was constructed and from this a smaller isolated interactome was identified. Five genes from this interactome were selected for genetic analysis. Twenty-one tag single-nucleotide polymorphisms (SNPs) which captured all common variation in these genes were genotyped in 10,196 Danes, and analyzed for association with NAFLD-related quantitative traits, type 2 diabetes (T2D), central obesity, and WHO-defined metabolic syndrome (MetS).Results273 genes were included in the protein-protein interaction analysis and EHHADH, ECHS1, HADHA, HADHB, and ACADL were selected for further examination. A total of 10 nominal statistical significant associations (P<0.05) to quantitative metabolic traits were identified. Also, the case-control study showed associations between variation in the five genes and T2D, central obesity, and MetS, respectively. Bonferroni adjustments for multiple testing negated all associations.ConclusionsUsing a bioinformatics approach we identified five candidate genes for NAFLD. However, we failed to provide evidence of associations with major effects between SNPs in these five genes and NAFLD-related quantitative traits, T2D, central obesity, and MetS.

Common variants in mendelian kidney disease genes and their association with renal function.

Many common genetic variants identified by genome-wide association studies for complex traits map to genes previously linked to rare inherited Mendelian disorders. A systematic analysis of common single-nucleotide polymorphisms (SNPs) in genes responsible for Mendelian diseases with kidney phenotypes has not been performed. We thus developed a comprehensive database of genes for Mendelian kidney conditions and evaluated the association between common genetic variants within these genes and kidney function in the general population. Using the Online Mendelian Inheritance in Man database, we identified 731 unique disease entries related to specific renal search terms and confirmed a kidney phenotype in 218 of these entries, corresponding to mutations in 258 genes. We interrogated common SNPs (minor allele frequency >5%) within these genes for association with the estimated GFR in 74,354 European-ancestry participants from the CKDGen Consortium. However, the top four candidate SNPs (rs6433115 at LRP2, rs1050700 at TSC1, rs249942 at PALB2, and rs9827843 at ROBO2) did not achieve significance in a stage 2 meta-analysis performed in 56,246 additional independent individuals, indicating that these common SNPs are not associated with estimated GFR. The effect of less common or rare variants in these genes on kidney function in the general population and disease-specific cohorts requires further research.

Lack of hypotriglyceridemic effect of gemfibrozil as a consequence of age-related changes in rat liver PPARalpha.

We have investigated if changes in hepatic lipid metabolism produced by old age are related to changes in liver peroxisome proliferator-activated receptor alpha (PPARalpha). Our results indicate that 18-month-old rats showed a marked decrease in the expression and activity of liver PPARalpha, as shown by significant reductions in PPARalpha mRNA, protein and binding activity, resulting in a reduction in the relative mRNA levels of PPARalpha target genes, such as liver-carnitine-palmitoyl transferase-I (CPT-I) and mitochondrial medium-chain acyl-CoA dehydrogenase (MCAD). Further, in accordance with a liver PPARalpha deficiency in old rats, treatment of old animals with a therapeutic dose of gemfibrozil (GFB) (3mg/kg per day, 21 days) was ineffective in reducing plasma triglyceride concentrations (TG), despite attaining a 50% reduction in TG when GFB was administered to young animals at the same dose and length of treatment. We hypothesize that the decrease in hepatic PPARalpha can be related to a state of leptin resistance present in old animals.