Origins and functional evolution of Y chromosomes across mammals.

Y chromosomes underlie sex determination in mammals, but their repeat-rich nature has hampered sequencing and associated evolutionary studies. Here we trace Y evolution across 15 representative mammals on the basis of high-throughput genome and transcriptome sequencing. We uncover three independent sex chromosome originations in mammals and birds (the outgroup). The original placental and marsupial (therian) Y, containing the sex-determining gene SRY, emerged in the therian ancestor approximately 180 million years ago, in parallel with the first of five monotreme Y chromosomes, carrying the probable sex-determining gene AMH. The avian W chromosome arose approximately 140 million years ago in the bird ancestor. The small Y/W gene repertoires, enriched in regulatory functions, were rapidly defined following stratification (recombination arrest) and erosion events and have remained considerably stable. Despite expression decreases in therians, Y/W genes show notable conservation of proto-sex chromosome expression patterns, although various Y genes evolved testis-specificities through differential regulatory decay. Thus, although some genes evolved novel functions through spatial/temporal expression shifts, most Y genes probably endured, at least initially, because of dosage constraints.

Bgee: integrating and comparing heterogeneous transcriptome data among species

Gene expression patterns are a key feature in understanding gene function, notably in development. Comparing gene expression patterns between animals is a major step in the study of gene function as well as of animal evolution. It also provides a link between genes and phenotypes. Thus we have developed Bgee, a database designed to compare expression patterns between animals, by implementing ontologies describing anatomies and developmental stages of species, and then designing homology relationships between anatomies and comparison criteria between developmental stages. To define homology relationships between anatomical features we have developed the software Homolonto, which uses a modified ontology alignment approach to propose homology relationships between ontologies. Bgee then uses these aligned ontologies, onto which heterogeneous expression data types are mapped. These already include microarrays and ESTs.

NCoR1 is a conserved physiological modulator of muscle mass and oxidative function.

Transcriptional coregulators control the activity of many transcription factors and are thought to have wide-ranging effects on gene expression patterns. We show here that muscle-specific loss of nuclear receptor corepressor 1 (NCoR1) in mice leads to enhanced exercise endurance due to an increase of both muscle mass and of mitochondrial number and activity. The activation of selected transcription factors that control muscle function, such as MEF2, PPARβ/δ, and ERRs, underpins these phenotypic alterations. NCoR1 levels are decreased in conditions that require fat oxidation, resetting transcriptional programs to boost oxidative metabolism. Knockdown of gei-8, the sole C. elegans NCoR homolog, also robustly increased muscle mitochondria and respiration, suggesting conservation of NCoR1 function. Collectively, our data suggest that NCoR1 plays an adaptive role in muscle physiology and that interference with NCoR1 action could be used to improve muscle function.

Bristle pattern formation in the thorax formation of Drosophila: a systems level approach

A fundamental question in developmental biology is how tissues are patterned to give rise to differentiated body structures with distinct morphologies. The Drosophila wing disc offers an accessible model to understand epithelial spatial patterning. It has been studied extensively using genetic and molecular approaches. Bristle patterns on the thorax, which arise from the medial part of the wing disc, are a classical model of pattern formation, dependent on a pre-pattern of trans-activators and –repressors. Despite of decades of molecular studies, we still only know a subset of the factors that determine the pre-pattern. We are applying a novel and interdisciplinary approach to predict regulatory interactions in this system. It is based on the description of expression patterns by simple logical relations (addition, subtraction, intersection and union) between simple shapes (graphical primitives). Similarities and relations between primitives have been shown to be predictive of regulatory relationships between the corresponding regulatory factors in other Systems, such as the Drosophila egg. Furthermore, they provide the basis for dynamical models of the bristle-patterning network, which enable us to make even more detailed predictions on gene regulation and expression dynamics. We have obtained a data-set of wing disc expression patterns which we are now processing to obtain average expression patterns for each gene. Through triangulation of the images we can transform the expression patterns into vectors which can easily be analysed by Standard clustering methods. These analyses will allow us to identify primitives and regulatory interactions. We expect to identify new regulatory interactions and to understand the basic Dynamics of the regulatory network responsible for thorax patterning. These results will provide us with a better understanding of the rules governing gene regulatory networks in general, and provide the basis for future studies of the evolution of the thorax-patterning network in particular.

Molecular changes underlying mammalian phenotypic innovation

Mammals are characterized by specific phenotypic traits that include lactation, hair, and relatively large brains with unique structures. Individual mammalian lineages have, in turn, evolved characteristic traits that distinguish them from others. These include obvious anatom¬ical differences but also differences related to reproduction, life span, cognitive abilities, be¬havior. and disease susceptibility. However, the molecular basis of the diverse mammalian phenotypes and the selective pressures that shaped their evolution remain largely unknown. In the first part of my thesis, I analyzed the genetic factors associated with the origin of a unique mammalian phenotype lactation and I studied the selective pressures that forged the transition from oviparity to viviparity. Using a comparative genomics approach and evolutionary simulations, I showed that the emergence of lactation, as well as the appear¬ance of the casein gene family, significantly reduced selective pressure on the major egg-yolk proteins (the vitellogenin family). This led to a progressive loss of vitellogenins, which - in oviparous species - act as storage proteins for lipids, amino acids, phosphorous and calcium in the isolated egg. The passage to internal fertilization and placentation in therian mam¬mals rendered vitellogenins completely dispensable, which ended in the loss of the whole gene family in this lineage. As illustrated by the vitellogenin study, changes in gene content are one possible underlying factor for the evolution of mammalian-specific phenotypes. However, more subtle genomic changes, such as mutations in protein-coding sequences, can also greatly affect the phenotypes. In particular, it was proposed that changes at the level of gene reg¬ulation could underlie many (or even most) phenotypic differences between species. In the second part of my thesis, I participated in a major comparative study of mammalian tissue transcriptomes, with the goal of understanding how evolutionary forces affected expression patterns in the past 200 million years of mammalian evolution. I showed that, while com¬parisons of gene expressions are in agreement with the known species phylogeny, the rate of expression evolution varies greatly among lineages. Species with low effective population size, such as monotremes and hominoids, showed significantly accelerated rates of gene expression evolution. The most likely explanation for the high rate of gene expression evolution in these lineages is the accumulation of mildly deleterious mutations in regulatory regions, due to the low efficiency of purifying selection. Thus, our observations are in agreement with the nearly neutral theory of molecular evolution. I also describe substantial differences in evolutionary rates between tissues, with brain being the most constrained (especially in primates) and testis significantly accelerated. The rate of gene expression evolution also varies significantly between chromosomes. In particular, I observed an acceleration of gene expression changes on the X chromosome, probably as a result of adaptive processes associated with the origin of therian sex chromosomes. Lastly, I identified several individual genes as well as co-regulated expression modules that have undergone lineage specific expression changes and likely under¬lie various phenotypic innovations in mammals. The methods developed during my thesis, as well as the comprehensive gene content analyses and transcriptomics datasets made available by our group, will likely prove to be useful for further exploratory analyses of the diverse mammalian phenotypes.

Multiple non-collinear TF-map alignments of promoter regions

Background: The analysis of the promoter sequence of genes with similar expression patterns isa basic tool to annotate common regulatory elements. Multiple sequence alignments are on thebasis of most comparative approaches. The characterization of regulatory regions from coexpressedgenes at the sequence level, however, does not yield satisfactory results in manyoccasions as promoter regions of genes sharing similar expression programs often do not shownucleotide sequence conservation.Results: In a recent approach to circumvent this limitation, we proposed to align the maps ofpredicted transcription factors (referred as TF-maps) instead of the nucleotide sequence of tworelated promoters, taking into account the label of the corresponding factor and the position in theprimary sequence. We have now extended the basic algorithm to permit multiple promotercomparisons using the progressive alignment paradigm. In addition, non-collinear conservationblocks might now be identified in the resulting alignments. We have optimized the parameters ofthe algorithm in a small, but well-characterized collection of human-mouse-chicken-zebrafishorthologous gene promoters.Conclusion: Results in this dataset indicate that TF-map alignments are able to detect high-levelregulatory conservation at the promoter and the 3'UTR gene regions, which cannot be detectedby the typical sequence alignments. Three particular examples are introduced here to illustrate thepower of the multiple TF-map alignments to characterize conserved regulatory elements inabsence of sequence similarity. We consider this kind of approach can be extremely useful in thefuture to annotate potential transcription factor binding sites on sets of co-regulated genes fromhigh-throughput expression experiments.

Transcription factor map alignment of promoter regions

We address the problem of comparing and characterizing the promoter regions of genes with similar expression patterns. This remains a challenging problem in sequence analysis, because often the promoter regions of co-expressed genes do not show discernible sequence conservation. In our approach, thus, we have not directly compared the nucleotide sequence of promoters. Instead, we have obtained predictions of transcription factor binding sites, annotated the predicted sites with the labels of the corresponding binding factors, and aligned the resulting sequences of labels—to which we refer here as transcription factor maps (TF-maps). To obtain the global pairwise alignment of two TF-maps, we have adapted an algorithm initially developed to align restriction enzyme maps. We have optimized the parameters of the algorithm in a small, but well-curated, collection of human–mouse orthologous gene pairs. Results in this dataset, as well as in an independent much larger dataset from the CISRED database, indicate that TF-map alignments are able to uncover conserved regulatory elements, which cannot be detected by the typical sequence alignments.

Peroxisome proliferator-activated receptor beta regulates acyl-CoA synthetase 2 in reaggregated rat brain cell cultures.

Peroxisome proliferator-activated receptors (PPARs) are nuclear hormone receptors that regulate the expression of many genes involved in lipid metabolism. The biological roles of PPARalpha and PPARgamma are relatively well understood, but little is known about the function of PPARbeta. To address this question, and because PPARbeta is expressed to a high level in the developing brain, we used reaggregated brain cell cultures prepared from dissociated fetal rat telencephalon as experimental model. In these primary cultures, the fetal cells initially form random aggregates, which progressively acquire a tissue-specific pattern resembling that of the brain. PPARs are differentially expressed in these aggregates, with PPARbeta being the prevalent isotype. PPARalpha is present at a very low level, and PPARgamma is absent. Cell type-specific expression analyses revealed that PPARbeta is ubiquitous and most abundant in some neurons, whereas PPARalpha is predominantly astrocytic. We chose acyl-CoA synthetases (ACSs) 1, 2, and 3 as potential target genes of PPARbeta and first analyzed their temporal and cell type-specific pattern. This analysis indicated that ACS2 and PPARbeta mRNAs have overlapping expression patterns, thus designating the ACS2 gene as a putative target of PPARbeta. Using a selective PPARbeta activator, we found that the ACS2 gene is transcriptionally regulated by PPARbeta, demonstrating a role for PPARbeta in brain lipid metabolism.

Evolutionary origin and functions of retrogene introns.

Retroposed genes (retrogenes) originate via the reverse transcription of mature messenger RNAs from parental source genes and are therefore usually devoid of introns. Here, we characterize a particular set of mammalian retrogenes that acquired introns upon their emergence and thus represent rare cases of intron gain in mammals. We find that although a few retrogenes evolved introns in their coding or 3' untranslated regions (untranslated region, UTR), most introns originated together with untranslated exons in the 5' flanking regions of the retrogene insertion site. They emerged either de novo or through fusions with 5' UTR exons of host genes into which the retrogenes inserted. Generally, retrogenes with introns display high transcription levels and show broader spatial expression patterns than other retrogenes. Our experimental expression analyses of individual intron-containing retrogenes show that 5' UTR introns may indeed promote higher expression levels, at least in part through encoded regulatory elements. By contrast, 3' UTR introns may lead to downregulation of expression levels via nonsense-mediated decay mechanisms. Notably, the majority of retrogenes with introns in their 5' flanks depend on distant, sometimes bidirectional CpG dinucleotide-enriched promoters for their expression that may be recruited from other genes in the genomic vicinity. We thus propose a scenario where the acquisition of new 5' exon-intron structures was directly linked to the recruitment of distant promoters by these retrogenes, a process potentially facilitated by the presence of proto-splice sites in the genomic vicinity of retrogene insertion sites. Thus, the primary role and selective benefit of new 5' introns (and UTR exons) was probably initially to span the often substantial distances to potent CpG promoters driving retrogene transcription. Later in evolution, these introns then obtained additional regulatory roles in fine tuning retrogene expression levels. Our study provides novel insights regarding mechanisms underlying the origin of new introns, the evolutionary relevance of intron gain, and the origin of new gene promoters.

Segmental copy number variation shapes tissue transcriptomes.

Copy number variation (CNV) is a key source of genetic diversity, but a comprehensive understanding of its phenotypic effect is only beginning to emerge. We have generated a CNV map in wild mice and classical inbred strains. Genome-wide expression data from six major organs show not only that expression of genes within CNVs tend to correlate with copy number changes, but also that CNVs influence the expression of genes in their vicinity, an effect that extends up to half a megabase. Genes within CNVs show lower expression and more specific spatial expression patterns than genes mapping elsewhere. Our analyses reveal differential constraint on copy number changes of genes expressed in different tissues. Dosage alterations of brain-expressed genes are less frequent than those of other genes and are buffered by tighter transcriptional regulation. Our study provides initial evidence that CNVs shape tissue transcriptomes on a global scale.

A novel bioinformatics pipeline to discover genes related to arbuscular mycorrhizal symbiosis based on their evolutionary conservation pattern among higher plants.

BACKGROUND: Genes involved in arbuscular mycorrhizal (AM) symbiosis have been identified primarily by mutant screens, followed by identification of the mutated genes (forward genetics). In addition, a number of AM-related genes has been identified by their AM-related expression patterns, and their function has subsequently been elucidated by knock-down or knock-out approaches (reverse genetics). However, genes that are members of functionally redundant gene families, or genes that have a vital function and therefore result in lethal mutant phenotypes, are difficult to identify. If such genes are constitutively expressed and therefore escape differential expression analyses, they remain elusive. The goal of this study was to systematically search for AM-related genes with a bioinformatics strategy that is insensitive to these problems. The central element of our approach is based on the fact that many AM-related genes are conserved only among AM-competent species. RESULTS: Our approach involves genome-wide comparisons at the proteome level of AM-competent host species with non-mycorrhizal species. Using a clustering method we first established orthologous/paralogous relationships and subsequently identified protein clusters that contain members only of the AM-competent species. Proteins of these clusters were then analyzed in an extended set of 16 plant species and ranked based on their relatedness among AM-competent monocot and dicot species, relative to non-mycorrhizal species. In addition, we combined the information on the protein-coding sequence with gene expression data and with promoter analysis. As a result we present a list of yet uncharacterized proteins that show a strongly AM-related pattern of sequence conservation, indicating that the respective genes may have been under selection for a function in AM. Among the top candidates are three genes that encode a small family of similar receptor-like kinases that are related to the S-locus receptor kinases involved in sporophytic self-incompatibility. CONCLUSIONS: We present a new systematic strategy of gene discovery based on conservation of the protein-coding sequence that complements classical forward and reverse genetics. This strategy can be applied to diverse other biological phenomena if species with established genome sequences fall into distinguished groups that differ in a defined functional trait of interest.

Parent-of-origin-specific allelic associations among 106 genomic loci for age at menarche.

Age at menarche is a marker of timing of puberty in females. It varies widely between individuals, is a heritable trait and is associated with risks for obesity, type 2 diabetes, cardiovascular disease, breast cancer and all-cause mortality. Studies of rare human disorders of puberty and animal models point to a complex hypothalamic-pituitary-hormonal regulation, but the mechanisms that determine pubertal timing and underlie its links to disease risk remain unclear. Here, using genome-wide and custom-genotyping arrays in up to 182,416 women of European descent from 57 studies, we found robust evidence (P < 5 × 10(-8)) for 123 signals at 106 genomic loci associated with age at menarche. Many loci were associated with other pubertal traits in both sexes, and there was substantial overlap with genes implicated in body mass index and various diseases, including rare disorders of puberty. Menarche signals were enriched in imprinted regions, with three loci (DLK1-WDR25, MKRN3-MAGEL2 and KCNK9) demonstrating parent-of-origin-specific associations concordant with known parental expression patterns. Pathway analyses implicated nuclear hormone receptors, particularly retinoic acid and γ-aminobutyric acid-B2 receptor signalling, among novel mechanisms that regulate pubertal timing in humans. Our findings suggest a genetic architecture involving at least hundreds of common variants in the coordinated timing of the pubertal transition.

Arbuscular mycorrhiza-specific signaling in rice transcends the common symbiosis signaling pathway.

Knowledge about signaling in arbuscular mycorrhizal (AM) symbioses is currently restricted to the common symbiosis (SYM) signaling pathway discovered in legumes. This pathway includes calcium as a second messenger and regulates both AM and rhizobial symbioses. Both monocotyledons and dicotyledons form symbiotic associations with AM fungi, and although they differ markedly in the organization of their root systems, the morphology of colonization is similar. To identify and dissect AM-specific signaling in rice (Oryza sativa), we developed molecular phenotyping tools based on gene expression patterns that monitor various steps of AM colonization. These tools were used to distinguish common SYM-dependent and -independent signaling by examining rice mutants of selected putative legume signaling orthologs predicted to be perturbed both upstream (CASTOR and POLLUX) and downstream (CCAMK and CYCLOPS) of the central, calcium-spiking signal. All four mutants displayed impaired AM interactions and altered AM-specific gene expression patterns, therefore demonstrating functional conservation of SYM signaling between distant plant species. In addition, differential gene expression patterns in the mutants provided evidence for AM-specific but SYM-independent signaling in rice and furthermore for unexpected deviations from the SYM pathway downstream of calcium spiking.

Predictive value of the MGMT promoter methylation status in metastatic melanoma patients receiving first-line temozolomide plus bevacizumab in the trial SAKK 50/07.

The O6-methylguanine-DNA-methyltransferase (MGMT) promoter methylation status is a predictive parameter for the response of malignant gliomas to alkylating agents such as temozolomide. First clinical trials with temozolomide plus bevacizumab therapy in metastatic melanoma patients are ongoing, although the predictive value of the MGMT promoter methylation status in this setting remains unclear. We assessed MGMT promoter methylation in formalin-fixed, primary tumor tissue of metastatic melanoma patients treated with first-line temozolomide and bevacizumab from the trial SAKK 50/07 by methylation-specific polymerase chain reaction. In addition, the MGMT expression levels were also analyzed by MGMT immunohistochemistry. Eleven of 42 primary melanomas (26%) revealed a methylated MGMT promoter. Promoter methylation was significantly associated with response rates CR + PR versus SD + PD according to RECIST (response evaluation criteria in solid tumors) (p<0.05) with a trend to prolonged median progression-free survival (8.1 versus 3.4 months, p>0.05). Immunohistochemically different protein expression patterns with heterogeneous and homogeneous nuclear MGMT expression were identified. Negative MGMT expression levels were associated with overall disease stabilization CR+PR+SD versus PD (p=0.05). There was only a poor correlation between MGMT methylation and lack of MGMT expression. A significant proportion of melanomas have a methylated MGMT promoter. The MGMT promoter methylation status may be a promising predictive marker for temozolomide therapy in metastatic melanoma patients. Larger sample sizes may help to validate significant differences in survival type endpoints.

Involvement of the CXCR7/CXCR4/CXCL12 Axis in the Malignant Progression of Human Neuroblastoma.

Neuroblastoma (NB) is a typical childhood and heterogeneous neoplasm for which efficient targeted therapies for high-risk tumors are not yet identified. The chemokine CXCL12, and its receptors CXCR4 and CXCR7 have been involved in tumor progression and dissemination. While CXCR4 expression is associated to undifferentiated tumors and poor prognosis, the role of CXCR7, the recently identified second CXCL12 receptor, has not yet been elucidated in NB. In this report, CXCR7 and CXCL12 expressions were evaluated using a tissue micro-array including 156 primary and 56 metastatic NB tissues. CXCL12 was found to be highly associated to NB vascular and stromal structures. In contrast to CXCR4, CXCR7 expression was low in undifferentiated tumors, while its expression was stronger in matured tissues and specifically associated to differentiated neural tumor cells. As determined by RT-PCR, CXCR7 expression was mainly detected in N-and S-type NB cell lines, and was slightly induced upon NB cell differentiation in vitro. The relative roles of the two CXCL12 receptors were further assessed by overexpressing CXCR7 or CXCR4 receptor alone, or in combination, in the IGR-NB8 and the SH-SY5Y NB cell lines. In vitro functional analyses indicated that, in response to their common ligand, both receptors induced activation of ERK1/2 cascade, but not Akt pathway. CXCR7 strongly reduced in vitro growth, in contrast to CXCR4, and impaired CXCR4/CXCL12-mediated chemotaxis. Subcutaneous implantation of CXCR7-expressing NB cells showed that CXCR7 also significantly reduced in vivo growth. Moreover, CXCR7 affected CXCR4-mediated orthotopic growth in a CXCL12-producing environment. In such model, CXCR7, in association with CXCR4, did not induce NB cell metastatic dissemination. In conclusion, the CXCR7 and CXCR4 receptors revealed specific expression patterns and distinct functional roles in NB. Our data suggest that CXCR7 elicits anti-tumorigenic functions, and may act as a regulator of CXCR4/CXCL12-mediated signaling in NB.