Genes, environment and sex : perspectives from ectothermic vertebrates

RESUMELes modèles classiques sur l'évolution des chromosomes sexuels supposent que des gènes sexe- antagonistes s'accumulent sur les chromosomes sexuels, entraînant ainsi l'apparition d'une région non- recombinante, qui se répand progressivement en favorisant l'accumulation de mutations délétères. En accord avec cette théorie, les chromosomes sexuels que l'on observe aujourd'hui chez les mammifères et les oiseaux sont considérablement différenciés. En revanche, chez la plupart des vertébrés ectothermes, les chromosomes sexuels sont indifférenciés et il existe une impressionnante diversité de mécanismes de détermination du sexe. Au cours de cette thèse, j'ai étudié l'évolution des chromosomes sexuels chez les vertébrés ectothermes, en outre pour mieux comprendre ce contraste avec les vertébrés endothermes. L'hypothèse « high-turnover » postule que les chromosomes sexuels sont remplacés régulièrement à partir d'autosomes afin d'éviter leur dégénérescence. L'hypothèse « fountain-of-youth » propose que la recombinaison entre le chromosome X et le chromosome Y au sein de femelles XY empêche la dégénérescence. Les résultats de ma thèse, basés sur des études théoriques et empiriques, suggèrent que les deux processus peuvent être entraînés par l'environnement et ainsi jouent un rôle important dans l'évolution des chromosomes sexuels chez les vertébrés ectothermes.SUMMARYClassical models of sex-chromosome evolution assume that sexually antagonistic genes accumulate on sex chromosomes leading to a non-recombining region, which progressively expands and favors the accumulation of deleterious mutations. Concordant with this theory, sex chromosomes in extant mammals and birds are considerably differentiated. In most ectothermic vertebrates, such as frogs, however, sex chromosomes are undifferentiated and a striking diversity of sex determination systems is observed. This thesis was aimed to investigate this apparent contrast of sex chromosome evolution between endothermic and ectothermic vertebrates. The "high-turnover" hypothesis holds that sex chromosomes arose regularly from autosomes preventing decay. The "fountain-of-youth" hypothesis posits that sex chromosomes undergo episodic X-Y recombination in sex-reversed XY females, thereby purging ("rejuvenating") the Y chromosome. We suggest that both processes likely played an important role in sex chromosome evolution of ectothermic vertebrates. The literature largely views sex determination as a dichotomous process: individual sex is assumed to be determined either by genetic (genotypic sex determination, GSD) or by environmental factors (environmental sex determination, ESD), most often temperature (temperature sex determination, TSD). We endorsed an alternative view, which sees GSD and TSD as the ends of a continuum. The conservatism of molecular processes among different systems of sex determination strongly supports the continuum view. We proposed to define sex as a threshold trait underlain by a liability factor, and reaction norms allowing modeling interactions between genotypic and temperature effects. We showed that temperature changes (due to e.g., climatic changes or range expansions) are expected to provoke turnovers in sex-determination mechanisms maintaining homomorphic sex chromosomes. The balanced lethal system of crested newts might be the result of such a sex determination turnover, originating from two variants of ancient Y-chromosomes. Observations from a group of tree frogs, on the other hand, supported the 'fountain of youth' hypothesis. We then showed that low rates of sex- reversals in species with GSD might actually be adaptive considering joint effects of deleterious mutation purging and sexually antagonistic selection. Ongoing climatic changes are expected to threaten species with TSD by biasing population sex ratios. In contrast, species with GSD are implicitly assumed immune against such changes, because genetic systems are thought to necessarily produce even sex ratios. We showed that this assumption may be wrong and that sex-ratio biases by climatic changes may represent a previously unrecognized extinction threat for some GSD species.

Developmental regulation of expression of schizophrenia susceptibility genes in the primate hippocampal formation

The hippocampal formation is essential for normal memory function and is implicated in many neurodevelopmental, neurodegenerative and neuropsychiatric disorders. In particular, abnormalities in hippocampal structure and function have been identified in schizophrenic subjects. Schizophrenia has a strong polygenic component, but the role of numerous susceptibility genes in normal brain development and function has yet to be investigated. Here we described the expression of schizophrenia susceptibility genes in distinct regions of the monkey hippocampal formation during early postnatal development. We found that, as compared with other genes, schizophrenia susceptibility genes exhibit a differential regulation of expression in the dentate gyrus, CA3 and CA1, over the course of postnatal development. A number of these genes involved in synaptic transmission and dendritic morphology exhibit a developmental decrease of expression in CA3. Abnormal CA3 synaptic organization observed in schizophrenics might be related to some specific symptoms, such as loosening of association. Interestingly, changes in gene expression in CA3 might occur at a time possibly corresponding to the late appearance of the first clinical symptoms. We also found earlier changes in expression of schizophrenia susceptibility genes in CA1, which might be linked to prodromal psychotic symptoms. A number of schizophrenia susceptibility genes including APOE, BDNF, MTHFR and SLC6A4 are involved in other disorders, and thus likely contribute to nonspecific changes in hippocampal structure and function that must be combined with the dysregulation of other genes in order to lead to schizophrenia pathogenesis.

Apparent versus true gene expression changes of three hypoxia-related genes in autopsy derived tissue and the importance of normalisation.

The aim of our work was to show how a chosen normal-isation strategy can affect the outcome of quantitative gene expression studies. As an example, we analysed the expression of three genes known to be upregulated under hypoxic conditions: HIF1A, VEGF and SLC2A1 (GLUT1). Raw RT-qPCR data were normalised using two different strategies: a straightforward normalisation against a single reference gene, GAPDH, using the 2(-ΔΔCt) algorithm and a more complex normalisation against a normalisation factor calculated from the quantitative raw data from four previously validated reference genes. We found that the two different normalisation strategies revealed contradicting results: normalising against a validated set of reference genes revealed an upregulation of the three genes of interest in three post-mortem tissue samples (cardiac muscle, skeletal muscle and brain) under hypoxic conditions. Interestingly, we found a statistically significant difference in the relative transcript abundance of VEGF in cardiac muscle between donors who died of asphyxia versus donors who died from cardiac death. Normalisation against GAPDH alone revealed no upregulation but, in some instances, a downregulation of the genes of interest. To further analyse this discrepancy, the stability of all reference genes used were reassessed and the very low expression stability of GAPDH was found to originate from the co-regulation of this gene under hypoxic conditions. We concluded that GAPDH is not a suitable reference gene for the quantitative analysis of gene expression in hypoxia and that validation of reference genes is a crucial step for generating biologically meaningful data.

Hyperammonemia: regulation of argininosuccinate synthetase and argininosuccinate lyase genes in aggregating cell cultures of fetal rat brain.

Hyperammonemia in the brain leads to poorly understood alterations of nitric oxide (NO) synthesis. Arginine, the substrate of nitric oxide synthases, might be recycled from the citrulline produced with NO by argininosuccinate synthetase (AS) and argininosuccinate lyase (AL). The regulation of AS and AL genes during hyperammonemia is unknown in the brain. We used brain cell aggregates cultured from dissociated telencephalic cortex of rat embryos to analyze the regulation of AS and AL genes in hyperammonemia. Using RNase protection assay and non-radioactive in situ hybridization on aggregate cryosections, we show that both AS and AL genes are induced in astrocytes but not in neurons of aggregates exposed to 5 mM NH4Cl. Our work suggests that the hyperammonemic brain might increase its recycling of citrulline to arginine.

Evolutionary footprint of coevolving positions in genes.

MOTIVATION: The analysis of molecular coevolution provides information on the potential functional and structural implication of positions along DNA sequences, and several methods are available to identify coevolving positions using probabilistic or combinatorial approaches. The specific nucleotide or amino acid profile associated with the coevolution process is, however, not estimated, but only known profiles, such as the Watson-Crick constraint, are usually considered a priori in current measures of coevolution. RESULTS: Here, we propose a new probabilistic model, Coev, to identify coevolving positions and their associated profile in DNA sequences while incorporating the underlying phylogenetic relationships. The process of coevolution is modeled by a 16 × 16 instantaneous rate matrix that includes rates of transition as well as a profile of coevolution. We used simulated, empirical and illustrative data to evaluate our model and to compare it with a model of 'independent' evolution using Akaike Information Criterion. We showed that the Coev model is able to discriminate between coevolving and non-coevolving positions and provides better specificity and specificity than other available approaches. We further demonstrate that the identification of the profile of coevolution can shed new light on the process of dependent substitution during lineage evolution.

Transcriptional regulatory patterns of the myelin basic protein and malic enzyme genes by the thyroid hormone receptors alpha1 and beta1.

While there is evidence that the two ubiquitously expressed thyroid hormone (T3) receptors, TRalpha1 and TRbeta1, have distinct functional specificities, the mechanism by which they discriminate potential target genes remains largely unexplained. In this study, we demonstrate that the thyroid hormone response elements (TRE) from the malic enzyme and myelin basic protein genes (METRE and MBPTRE) respectively, are not functionally equivalent. The METRE, which is a direct repeat motif with a 4-base pair gap between the two half-site hexamers binds thyroid hormone receptor as a heterodimer with 9-cis-retinoic acid receptor (RXR) and mediates a high T3-dependent activation in response to TRalpha1 or TRbeta1 in NIH3T3 cells. In contrast, the MBPTRE, which consists of an inverted palindrome formed by two hexamers spaced by 6 base pairs, confers an efficient transactivation by TRbeta1 but a poor transactivation by TRalpha1. While both receptors form heterodimers with RXR on MBPTRE, the poor transactivation by TRalpha1 correlates also with its ability to bind efficiently as a monomer. This monomer, which is only observed with TRalpha1 bound to MBPTRE, interacts neither with N-CoR nor with SRC-1, explaining its functional inefficacy. However, in Xenopus oocytes, in which RXR proteins are not detectable, the transactivation mediated by TRalpha1 and TRbeta1 is equivalent and independent of a RXR supply, raising the question of the identity of the thyroid hormone receptor partner in these cells. Thus, in mammalian cells, the binding characteristics of TRalpha1 to MBPTRE (i.e. high monomer binding efficiency and low transactivation activity) might explain the particular pattern of T3 responsiveness of MBP gene expression during central nervous system development.

Identification of cancer/testis-antigen genes by massively parallel signature sequencing.

Massively parallel signature sequencing (MPSS) generates millions of short sequence tags corresponding to transcripts from a single RNA preparation. Most MPSS tags can be unambiguously assigned to genes, thereby generating a comprehensive expression profile of the tissue of origin. From the comparison of MPSS data from 32 normal human tissues, we identified 1,056 genes that are predominantly expressed in the testis. Further evaluation by using MPSS tags from cancer cell lines and EST data from a wide variety of tumors identified 202 of these genes as candidates for encoding cancer/testis (CT) antigens. Of these genes, the expression in normal tissues was assessed by RT-PCR in a subset of 166 intron-containing genes, and those with confirmed testis-predominant expression were further evaluated for their expression in 21 cancer cell lines. Thus, 20 CT or CT-like genes were identified, with several exhibiting expression in five or more of the cancer cell lines examined. One of these genes is a member of a CT gene family that we designated as CT45. The CT45 family comprises six highly similar (>98% cDNA identity) genes that are clustered in tandem within a 125-kb region on Xq26.3. CT45 was found to be frequently expressed in both cancer cell lines and lung cancer specimens. Thus, MPSS analysis has resulted in a significant extension of our knowledge of CT antigens, leading to the discovery of a distinctive X-linked CT-antigen gene family.

Epigenetic regulatory elements associate with specific histone modifications to prevent silencing of telomeric genes.

In eukaryotic cells, transgene expression levels may be limited by an unfavourable chromatin structure at the integration site. Epigenetic regulators are DNA sequences which may protect transgenes from such position effect. We evaluated different epigenetic regulators for their ability to protect transgene expression at telomeres, which are commonly associated to low or inconsistent expression because of their repressive chromatin environment. Although to variable extents, matrix attachment regions (MARs), ubiquitous chromatin opening element (UCOE) and the chicken cHS4 insulator acted as barrier elements, protecting a telomeric-distal transgene from silencing. MARs also increased the probability of silent gene reactivation in time-course experiments. Additionally, all MARs improved the level of expression in non-silenced cells, unlike other elements. MARs were associated to histone marks usually linked to actively expressed genes, especially acetylation of histone H3 and H4, suggesting that they may prevent the spread of silencing chromatin by imposing acetylation marks on nearby nucleosomes. Alternatively, an UCOE was found to act by preventing deposition of repressive chromatin marks. We conclude that epigenetic DNA elements used to enhance and stabilize transgene expression all have specific epigenetic signature that might be at the basis of their mode of action.

Revisit of multiple epiphyseal dysplasia: ethnic difference in genotypes and comparison of radiographic features linked to the COMP and MATN3 genes.

Multiple epiphyseal dysplasia (MED) is a genetically heterogeneous group of diseases characterized by variable degrees of epiphyseal abnormality primarily involving the hip and knee joints. The purpose of this study was to investigate the frequency of mutations in individuals with a clinical and radiographic diagnosis of MED and to test the hypothesis that characteristic radiological findings may be helpful in predicting the gene responsible. The radiographs of 74 Korean patients were evaluated by a panel of skeletal dysplasia experts. Six genes known to be associated with MED (COMP, MATN3, COL9A1, COL9A2, COL9A3, and DTDST) were screened by sequencing. Mutations were found in 55 of the 63 patients (87%). MATN3 mutations were found in 30 patients (55%), followed by COMP mutations in 23 (41%), and COL9A2 and DTDST mutations in one patient (2%) each. Comparisons of radiographic findings in patients with COMP and MATN3 mutations showed that albeit marked abnormalities in hip and knee joints were observed in both groups, the degree of involvement and the morphology of dysplastic epiphyses differed markedly. The contour of the pelvic acetabulum, the presence of metaphyseal vertical striations, and/or the brachydactyly of the hand were also found to be highly correlated with the genotypes. The study confirms that MATN3 and COMP are the genes most frequently responsible for MED and that subtle radiographic signs may give precious indications on which gene(s) should be prioritized for mutational screening in a given individual.

Etude des gènes codant les protéines de gouttelettes lipidiques chez les patients avec la stéatose hépatique[Study of genes encoding proteins of lipid droplets in patients with fatty liver]

Introduction: La prévalence de la «non-alcoholic fatty liver disease (NAFLD)» dans les pays industrialisés augment de manière exponentielle. La NAFLD se développe d'une simple stéatose hépatique jusqu'à l'hépatite, puis à la cirrhose. De plus, la stéatose hépatique est fréquemment accompagnée par une résistance à l'insuline, une des causes principales du diabète. Les lipides intermédiaires, tels que céramides et diacylglycérols, ont été décrits comme induisant la résistance à l'insuline. Cependant, nous avons démontré dans notre modèle de stéatose hépatique, que les souris présentant une invalidation de la protéine «microsomal triglyceride transfer protein» (Mtpp) au niveau hépatique, ne développent pas de résistance à l'insuline. Ceci suggère fortement l'existence d'autres mécanismes susceptibles d'induire la résistance à l'insuline. Résultats: Grâce à une analyse de Microarray, nous avons observé une augmentation de l'expression des gènes «cell-death inducing DFFA-like effector c (CIDEC)», «lipid storage droplet protein 5 (LSDP5)» et «Bernardinelli-Seip congenital lipodystrophy 2 homolog (Seipin)» dans le foie des souris Mttp. Ces gènes ont récemment été identifiés comme des protéines localisées autour des gouttelettes lipidiques. Nous avons également constaté que la souris Mttp développe plutôt une microstéatose (petites gouttelettes lipidiques) qu'une macrostéatose qui est normalement observée chez les patients avec NAFLD. Nous avons étudié l'expression des gènes associés aux gouttelettes lipidiques chez les patients obèses avec stéatose hépatique, avec ou sans résistance à l'insuline. Comparés aux sujets sains sans stéatose hépatique, les patients avec la stéatose ont une expression significativement plus élevée. De manière intéressante, les patients avec résistance à l'insuline ont une diminution de ces expressions. Conclusion : Ces données suggèrent que les gènes des gouttelettes lipidiques sont impliqués dans le développement de la stéatose hépatique chez l'homme et peut-être contribue à la mise en place de la résistance à l'insuline.

Mutational screening of splicing factor genes in cases with autosomal dominant retinitis pigmentosa.

PURPOSE: Mutations in genes encoding proteins from the tri-snRNP complex of the spliceosome account for more than 12% of cases of autosomal dominant retinitis pigmentosa (adRP). Although the exact mechanism by which splicing factor defects trigger photoreceptor death is not completely clear, their role in retinitis pigmentosa has been demonstrated by several genetic and functional studies. To test for possible novel associations between splicing factors and adRP, we screened four tri-snRNP splicing factor genes (EFTUD2, PRPF4, NHP2L1, and AAR2) as candidate disease genes. METHODS: We screened up to 303 patients with adRP from Europe and North America who did not carry known RP mutations. Exon-PCR and Sanger methods were used to sequence the NHP2L1 and AAR2 genes, while the sequences of EFTUD2 and PRPF4 were obtained by using long-range PCRs spanning coding and non-coding regions followed by next-generation sequencing. RESULTS: We detected novel missense changes in individual patients in the sequence of the genes PRPF4 and EFTUD2, but the role of these changes in relationship to disease could not be verified. In one other patient we identified a novel nucleotide substitution in the 5' untranslated region (UTR) of NHP2L1, which did not segregate with the disease in the family. CONCLUSIONS: The absence of clearly pathogenic mutations in the candidate genes screened in our cohort suggests that EFTUD2, PRPF4, NHP2L1, and AAR2 are either not involved in adRP or are associated with the disease in rare instances, at least as observed in this study in patients of European and North American origin.

Prospective study of insulin-like growth factor-I, insulin-like growth factor-binding protein 3, genetic variants in the IGF1 and IGFBP3 genes and risk of coronary artery disease.

Although experimental studies have suggested that insulin-like growth factor I (IGF-I) and its binding protein IGFBP-3 might have a role in the aetiology of coronary artery disease (CAD), the relevance of circulating IGFs and their binding proteins in the development of CAD in human populations is unclear. We conducted a nested case-control study, with a mean follow-up of six years, within the EPIC-Norfolk cohort to assess the association between circulating levels of IGF-I and IGFBP-3 and risk of CAD in up to 1,013 cases and 2,055 controls matched for age, sex and study enrolment date. After adjustment for cardiovascular risk factors, we found no association between circulating levels of IGF-I or IGFBP-3 and risk of CAD (odds ratio: 0.98 (95% Cl 0.90-1.06) per 1 SD increase in circulating IGF-I; odds ratio: 1.02 (95% Cl 0.94-1.12) for IGFBP-3). We examined associations between tagging single nucleotide polymorphisms (tSNPs) at the IGF1 and IGFBP3 loci and circulating IGF-I and IGFBP-3 levels in up to 1,133 cases and 2,223 controls and identified three tSNPs (rs1520220, rs3730204, rs2132571) that showed independent association with either circulating IGF-I or IGFBP-3 levels. In an assessment of 31 SNPs spanning the IGF1 or IGFBP3 loci, none were associated with risk of CAD in a meta-analysis that included EPIC-Norfolk and eight additional studies comprising up to 9,319 cases and 19,964 controls. Our results indicate that IGF-I and IGFBP-3 are unlikely to be importantly involved in the aetiology of CAD in human populations.

Loss of egg yolk genes in mammals and the origin of lactation and placentation

Embryonic development in nonmammalian vertebrates depends entirely on nutritional reserves that are predominantly derived from vitellogenin proteins and stored in egg yolk. Mammals have evolved new resources, such as lactation and placentation, to nourish their developing and early offspring. However, the evolutionary timing and molecular events associated with this major phenotypic transition are not known. By means of sensitive comparative genomics analyses and evolutionary simulations, we here show that the three ancestral vitellogenin-encoding genes were progressively lost during mammalian evolution (until around 30-70 million years ago, Mya) in all but the egg-laying monotremes, which have retained a functional vitellogenin gene. Our analyses also provide evidence that the major milk resource genes, caseins, which have similar functional properties as vitellogenins, appeared in the common mammalian ancestor approximately 200-310 Mya. Together, our data are compatible with the hypothesis that the emergence of lactation in the common mammalian ancestor and the development of placentation in eutherian and marsupial mammals allowed for the gradual loss of yolk-dependent nourishment during mammalian evolution

Use of constant denaturant capillary electrophoresis of pooled blood samples to identify single-nucleotide polymorphisms in the genes (Scnn1a and Scnn1b) encoding the alpha and beta subunits of the epithelial sodium channel.

BACKGROUND: The epithelial sodium channel (ENaC) is composed of three homologous subunits: alpha, beta, and gamma. Mutations in the Scnn1b and Scnn1g genes, which encode the beta and the gamma subunits of ENaC, cause a severe form of hypertension (Liddle syndrome). The contribution of genetic variants within the Scnn1a gene, which codes for the alpha subunit, has not been investigated. METHODS: We screened for mutations in the COOH termini of the alpha and beta subunits of ENaC. Blood from 184 individuals from 31 families participating in a study on the genetics of hypertension were analyzed. Exons 13 of Scnn1a and Scnn1b, which encode the second transmembrane segment and the COOH termini of alpha- and beta-ENaC, respectively, were amplified from pooled DNA samples of members of each family by PCR. Constant denaturant capillary electrophoresis (CDCE) was used to detect mutations in PCR products of the pooled DNA samples. RESULTS: The detection limit of CDCE for ENaC variants was 1%, indicating that all members of any family or up to 100 individuals can be analyzed in one CDCE run. CDCE profiles of the COOH terminus of alpha-ENaC in pooled family members showed that the 31 families belonged to four groups and identified families with genetic variants. Using this approach, we analyzed 31 rather than 184 samples. Individual CDCE analysis of members from families with different pooled CDCE profiles revealed five genotypes containing 1853G-->T and 1987A-->G polymorphisms. The presence of the mutations was confirmed by DNA sequencing. For the COOH terminus of beta-ENaC, only one family showed a different CDCE profile. Two members of this family (n = 5) were heterozygous at 1781C-->T (T594M). CONCLUSION: CDCE rapidly detects point mutations in these candidate disease genes.

Identification of genes potentially involved in solute stress response in Sphingomonas wittichii RW1 by transposon mutant recovery.

The term water stress refers to the effects of low water availability on microbial growth and physiology. Water availability has been proposed as a major constraint for the use of microorganisms in contaminated sites with the purpose of bioremediation. Sphingomonas wittichii RW1 is a bacterium capable of degrading the xenobiotic compounds dibenzofuran and dibenzo-p-dioxin, and has potential to be used for targeted bioremediation. The aim of the current work was to identify genes implicated in water stress in RW1 by means of transposon mutagenesis and mutant growth experiments. Conditions of low water potential were mimicked by adding NaCl to the growth media. Three different mutant selection or separation method were tested which, however recovered different mutants. Recovered transposon mutants with poorer growth under salt-induced water stress carried insertions in genes involved in proline and glutamate biosynthesis, and further in a gene putatively involved in aromatic compound catabolism. Transposon mutants growing poorer on medium with lowered water potential also included ones that had insertions in genes involved in more general functions such as transcriptional regulation, elongation factor, cell division protein, RNA polymerase β or an aconitase.