Developmental and Anatomical Constraints on Vertebrate Genome Evolution

Pendant ma thèse de doctorat, j'ai utilisé des espèces modèles, comme la souris et le poisson-zèbre, pour étudier les facteurs qui affectent l'évolution des gènes et leur expression. Plus précisément, j'ai montré que l'anatomie et le développement sont des facteurs clés à prendre en compte, car ils influencent la vitesse d'évolution de la séquence des gènes, l'impact sur eux de mutations (i.e. la délétion du gène est-elle létale ?), et leur tendance à se dupliquer. Où et quand il est exprimé impose à un gène certaines contraintes ou au contraire lui donne des opportunités d'évoluer. J'ai pu comparer ces tendances aux modèles classiques d'évolution de la morphologie, que l'on pensait auparavant refléter directement les contraintes s'appliquant sur le génome. Nous avons montré que les contraintes entre ces deux niveaux d'organisation ne peuvent pas être transférées simplement : il n'y a pas de lien direct entre la conservation du génotype et celle de phénotypes comme la morphologie. Ce travail a été possible grâce au développement d'outils bioinformatiques. Notamment, j'ai travaillé sur le développement de la base de données Bgee, qui a pour but de comparer l'expression des gènes entre différentes espèces de manière automatique et à large échelle. Cela implique une formalisation de l'anatomie, du développement et de concepts liés à l'homologie grâce à l'utilisation d'ontologies. Une intégration cohérente de données d'expression hétérogènes (puces à ADN, marqueurs de séquence exprimée, hybridations in situ) a aussi été nécessaire. Cette base de données est mise à jour régulièrement et disponible librement. Elle devrait contribuer à étendre les possibilités de comparaison de l'expression des gènes entre espèces pour des études d'évo-devo (évolution du développement) et de génomique. During my PhD, I used model species of vertebrates, such as mouse and zebrafish, to study factors affecting the evolution of genes and their expression. More precisely I have shown that anatomy and development are key factors to take into account, influencing the rate of gene sequence evolution, the impact of mutations (i.e. is the deletion of a gene lethal?), and the propensity of a gene to duplicate. Where and when genes are expressed imposes constraints, or on the contrary leaves them some opportunity to evolve. We analyzed these patterns in relation to classical models of morphological evolution in vertebrates, which were previously thought to directly reflect constraints on the genomes. We showed that the patterns of evolution at these two levels of organization do not translate smoothly: there is no direct link between the conservation of genotype and phenotypes such as morphology. This work was made possible by the development of bioinformatics tools. Notably, I worked on the development of the database Bgee, which aims at comparing gene expression between different species in an automated and large-scale way. This involves the formalization of anatomy, development, and concepts related to homology, through the use of ontologies. A coherent integration of heterogeneous expression data (microarray, expressed sequence tags, in situ hybridizations) is also required. This database is regularly updated and freely available. It should contribute to extend the possibilities for comparison of gene expression between species in evo-devo and genomics studies.

Evolution of vitellogenin genes: comparative analysis of the nucleotide sequences downstream of the transcription initiation site of four Xenopus laevis and one chicken gene.

Electron microscopic analysis of heteroduplexes between the most distantly related Xenopus vitellogenin genes (A genes X B genes) has revealed the distribution of homologous regions that have been preferentially conserved after the duplication events that gave rise to the multigene family in Xenopus laevis. DNA sequence analysis was limited to the region downstream of the transcription initiation site of the Xenopus genes A1, B1 and B2 and a comparison with the Xenopus A2 and the major chicken vitellogenin gene is presented. Within the coding regions of the first three exons, nucleotide substitutions resulting in amino acid changes accumulate at a rate similar to that observed in globin genes. This suggests that the duplication event which led to the formation of the A and B ancestral genes in Xenopus laevis occurred about 150 million years ago. Homologous exons of the A1-A2 and B1-B2 gene pairs, which formed about 30 million years ago, show a quite similar sequence divergence. In contrast, A1-A2 homologous introns seem to have evolved much faster than their B1-B2 counterparts.

Interpretation of array comparative genome hybridization data: a major challenge.

The advent and application of high-resolution array-based comparative genome hybridization (array CGH) has led to the detection of large numbers of copy number variants (CNVs) in patients with developmental delay and/or multiple congenital anomalies as well as in healthy individuals. The notion that CNVs are also abundantly present in the normal population challenges the interpretation of the clinical significance of detected CNVs in patients. In this review we will illustrate a general clinical workflow based on our own experience that can be used in routine diagnostics for the interpretation of CNVs.

16p11.2 600 kb Duplications confer risk for typical and atypical Rolandic epilepsy.

Rolandic epilepsy (RE) is the most common idiopathic focal childhood epilepsy. Its molecular basis is largely unknown and a complex genetic etiology is assumed in the majority of affected individuals. The present study tested whether six large recurrent copy number variants at 1q21, 15q11.2, 15q13.3, 16p11.2, 16p13.11 and 22q11.2 previously associated with neurodevelopmental disorders also increase risk of RE. Our association analyses revealed a significant excess of the 600 kb genomic duplication at the 16p11.2 locus (chr16: 29.5-30.1 Mb) in 393 unrelated patients with typical (n = 339) and atypical (ARE; n = 54) RE compared with the prevalence in 65,046 European population controls (5/393 cases versus 32/65,046 controls; Fisher's exact test P = 2.83 × 10(-6), odds ratio = 26.2, 95% confidence interval: 7.9-68.2). In contrast, the 16p11.2 duplication was not detected in 1738 European epilepsy patients with either temporal lobe epilepsy (n = 330) and genetic generalized epilepsies (n = 1408), suggesting a selective enrichment of the 16p11.2 duplication in idiopathic focal childhood epilepsies (Fisher's exact test P = 2.1 × 10(-4)). In a subsequent screen among children carrying the 16p11.2 600 kb rearrangement we identified three patients with RE-spectrum epilepsies in 117 duplication carriers (2.6%) but none in 202 carriers of the reciprocal deletion. Our results suggest that the 16p11.2 duplication represents a significant genetic risk factor for typical and atypical RE.

Marked hemiatrophy in carriers of Duchenne muscular dystrophy.

OBJECTIVE: To describe the clinical and molecular genetic findings in 2 carriers of Duchenne muscular dystrophy (DMD) who exhibited marked hemiatrophy. Duchenne muscular dystrophy is an X-linked disorder in which affected male patients harbor mutations in the dystrophin gene. Female patients with heterozygous mutations may be manifesting carriers. DESIGN: Case study. SETTING: Neurology clinic. PATIENTS: Two manifesting carriers of DMD. INTERVENTIONS: Clinical and radiologic examinations along with histologic and molecular investigations. RESULTS: Both patients had marked right-sided hemiatrophy on examination with radiologic evidence of muscle atrophy and fatty replacement on the affected side. In each case, histologic analysis revealed a reduction in dystrophin staining on the right side. Genetic analysis of the dystrophin gene revealed a tandem exonic duplication in patient 1 and a multiexonic deletion in patient 2 with no further point mutations identified on the other chromosome. CONCLUSIONS: Marked hemiatrophy can occur in DMD manifesting carriers. This is likely to result from a combination of skewed X-inactivation and somatic mosaicism.

The serine repeat antigen (SERA) gene family phylogeny in Plasmodium: the impact of GC content and reconciliation of gene and species trees.

Plasmodium falciparum is the parasite responsible for the most acute form of malaria in humans. Recently, the serine repeat antigen (SERA) in P. falciparum has attracted attention as a potential vaccine and drug target, and it has been shown to be a member of a large gene family. To clarify the relationships among the numerous P. falciparum SERAs and to identify orthologs to SERA5 and SERA6 in Plasmodium species affecting rodents, gene trees were inferred from nucleotide and amino acid sequence data for 33 putative SERA homologs in seven different species. (A distance method for nucleotide sequences that is specifically designed to accommodate differing GC content yielded results that were largely compatible with the amino acid tree. Standard-distance and maximum-likelihood methods for nucleotide sequences, on the other hand, yielded gene trees that differed in important respects.) To infer the pattern of duplication, speciation, and gene loss events in the SERA gene family history, the resulting gene trees were then "reconciled" with two competing Plasmodium species tree topologies that have been identified by previous phylogenetic studies. Parsimony of reconciliation was used as a criterion for selecting a gene tree/species tree pair and provided (1) support for one of the two species trees and for the core topology of the amino acid-derived gene tree, (2) a basis for critiquing fine detail in a poorly resolved region of the gene tree, (3) a set of predicted "missing genes" in some species, (4) clarification of the relationship among the P. falciparum SERA, and (5) some information about SERA5 and SERA6 orthologs in the rodent malaria parasites. Parsimony of reconciliation and a second criterion--implied mutational pattern at two key active sites in the SERA proteins-were also seen to be useful supplements to standard "bootstrap" analysis for inferred topologies.

Estructura de datos avanzadas : clases disjuntas, montículos, árboles de búsqueda

En aquest treball s'amplia la implementació en Java de les estructures de dades iniciada per Esteve Mariné, utilitzant el seu disseny bàsic. Concretament, s'ha fet la programació de les estructures de a) classes disjuntes, utilitzant els algorismes de llistes encadenades i amb estructura d'arbre, b) monticles, amb els algorismes binari, binomial i de Fibonacci, i c) arbres de recerca basats en l'algorisme d'arbre binari vermell-negre, el qual complementa els dos ja existents amb algorismes d'encadenaments i AVL. Per a examinar l'evolució de les estructures, s'ha preparat un visualitzador gràfic interactiu amb l'usuari que permet fer les operacions bàsiques de l'estructura. Amb aquest entorn és possible desar les estructures, tornar a reproduir-les i desfer i tornar a repetir les operacions fetes sobre l'estructura. Finalment, aporta una metodologia, amb visualització mitjançant gràfics, de l'avaluació comparativa dels algorismes implementats, que permet modificar els paràmetres d'avaluació com ara nombre d'elements que s'han de tractar, algorismes que s'han de comparar i nombre de repeticions. Les dades obtingudes es poden exportar per a analitzar-les posteriorment.

The amphioxus genome and the evolution of the chordate karyotype.

Lancelets ('amphioxus') are the modern survivors of an ancient chordate lineage, with a fossil record dating back to the Cambrian period. Here we describe the structure and gene content of the highly polymorphic approximately 520-megabase genome of the Florida lancelet Branchiostoma floridae, and analyse it in the context of chordate evolution. Whole-genome comparisons illuminate the murky relationships among the three chordate groups (tunicates, lancelets and vertebrates), and allow not only reconstruction of the gene complement of the last common chordate ancestor but also partial reconstruction of its genomic organization, as well as a description of two genome-wide duplications and subsequent reorganizations in the vertebrate lineage. These genome-scale events shaped the vertebrate genome and provided additional genetic variation for exploitation during vertebrate evolution.

Transcriptional regulation of RNA polymerase III in mammalian cells

L'ARN Polymérase III (Pol III) transcrit un ensemble de petits ARN non traduits impliqués dans des processus cellulaires tels que la biosynthèse des protéines, la maturation des ARNs ou le contrôle transcriptionnel. De ce fait, la Pol III joue un rôle important dans la régulation de la croissance et la prolifération cellulaire. L'initiation de la transcription par la Pol III nécessite l'interaction entre des facteurs de transcription et le complexe de la Pol III lui-même. Un sous- complexe de la Pol III, composé de 3 sous-unités, HsRPC3, HsRPC6 et HsRPC7 sert d'intermédiaire dans cette interaction. Dans cette étude, nous avons caractérisé une nouvelle sous-unité de la Pol III, HsRPC7-Like, homologue à HsRPC7. Nous avons montré que ces deux homologues se trouvent spécifiquement chez les vertébrés. Ils proviennent d'un ancêtre commun qui, après duplication il y a 600 millions d'années, a donné naissance à ces deux paralogues. Dans les cellules humaines, deux formes de Pol III coexistent : l'une contientt HsRPC7, l'autre HsRPC7-Like. Nous avons localisé, à l'échelle du génome entier, la présence de ces deux formes de Pol III dans des cellules humaines et dans le foie de souris. Les deux sous-unités ont démontré des caractéristiques identiques, suggérant qu'elles possèdent des fonctions similaires. Cependant, nous avons analysé les motifs d'expression des gènes codant pour RPC7 et RPC7-Like dans des lignées cellulaires dans des conditions variées telles que la concentration de sérum et la densité cellulaire, ainsi que les motifs d'expression dans le foie de souris et des cellules d'hépatocarcinome de souris. Nos résultats suggèrent que l'expression de ces deux sous-untiés varie en fonction de l'activité de prolifération de la cellule. - RNA polymerase III (Pol III) transcribes a set of genes coding for short untranslated RNAs involved in essential cellular processes as for example protein biosynthesis, RNA maturation, and transcriptional control. Thereby Pol III plays an important role in regulating cell growth and proliferation. Initiation of Pol III transcription requires interactions between transcription factors and the Pol III core complex. A Pol III sub-complex composed of three subunits, HsRPC3, HsRPC6, and HsRPC7 mediates this interaction. In this study, we have characterized a new Pol III subunit, HsRPC7-Like, an homologue of HsRPC7. We have shown that these two homologues are specific to vertebrates and originate from an ancestor gene that duplicated 600 mio years ago to give birth to two paralogues. In human cells, two forms of Pol III coexist, one containing HsRPC7 and the other HsRPC7-Like. We have localized, genome-wide, these two Pol III forms in human cells and mouse liver. Both subunits were found on all types of Pol III genes, suggesting that they share similar function. However, we analysed the expression patterns of the RPC7 and RPC7-Like coding genes under various conditions of serum concentration and cell density in different cell lines, as well as expression patterns in mouse liver and mouse hepatocarcinoma cells. Our results suggest that the expression of these two subunits varies with the proliferation rate of the cell.

Mutations in CNNM4 cause recessive cone-rod dystrophy with amelogenesis imperfecta.

Cone-rod dystrophies are inherited dystrophies of the retina characterized by the accumulation of deposits mainly localized to the cone-rich macular region of the eye. Dystrophy can be limited to the retina or be part of a syndrome. Unlike nonsyndromic cone-rod dystrophies, syndromic cone-rod dystrophies are genetically heterogeneous with mutations in genes encoding structural, cell-adhesion, and transporter proteins. Using a genome-wide single-nucleotide polymorphism (SNP) haplotype analysis to fine map the locus and a gene-candidate approach, we identified homozygous mutations in the ancient conserved domain protein 4 gene (CNNM4) that either generate a truncated protein or occur in highly conserved regions of the protein. Given that CNNM4 is implicated in metal ion transport, cone-rod dystrophy and amelogenesis imperfecta may originate from abnormal ion homeostasis.

The Intra-Alpine terrain: A Paleotethyan remnant in the Alpine Variscides

A Cordilleran type evolution is proposed for the Variscan orogen of middle Europe. This orogenesis is regarded as mainly evolving through terrain accretion and subsequent collapse of the overthickened crust. A major terrain accretion took place between late Devonian and early Carboniferous when the Intra-Alpine terrain collided with the Ligerian-Moldanubian active margin. This terrain is regarded as being a segment of the northern margin of Paleotethys. Oblique subduction of Paleotethys under the newly accreted terrain is responsible for the voluminous calc-alkaline magmatism in late Carboniferous. The Paleotethys subduction has generated a lateral displacement of the eastern part of the Intra-Alpine terrain inducing a duplication of its western end. The late Carboniferous closure of Paleotethys in middle Europe is not found eastward where this closure happened only in early-Triassic times, following the simultaneous opening of the Neotethys ocean and the Meliata back-arc. Palinspastic models of the western Tethyan realm are proposed from the Carboniferous to early Jurassic.

Copy number variation and chromatin structure

The functional consequences of structural variation in the human genome range from adaptation, to phenotypic variation, to predisposition to diseases. Copy number variation (CNV) was shown to influence the phenotype by modifying, in a somewhat dose-dependent manner, the expression of genes that map within them, as well as that of genes located on their flanks. To assess the possible mechanism(s) behind this neighboring effect, we compared histone modification status of cell lines from patients affected by Williams-Beuren, Williams-Beuren region duplication, Smith-Magenis or DiGeorge Syndrome and control individuals using a high-throughput version of chromatin immuno-precipitation method (ChIP), called ChlP-seq. We monitored monomethylation of lysine K20 on histone H4 and trimethylation of lysine K27 on histone H3, as proxies for open and condensed chromatin, respectively. Consistent with the changes in expression levels observed for multiple genes mapping on the entire length of chromosomes affected by structural variants, we also detected regions with modified histone status between samples, up- and downstream from the critical regions, up to the end of the rearranged chromosome. We also gauged the intrachromosomal interactions of these cell lines utilizing chromosome conformation capture (4C-seq) technique. We observed that a set of genes flanking the Williams-Beuren Syndrome critical region (WBSCR) were often looping together, possibly forming an interacting cluster with each other and the WBSCR. Deletion of the WBSCR disrupts the expression of this group of flanking genes, as well as long-range interactions between them and the rearranged interval. We conclude, that large genomic rearrangements can lead to changes in the state of the chromatin spreading far away from the critical region, thus possibly affecting expression globally and as a result modifying the phenotype of the patients. - Les conséquences fonctionnelles des variations structurelles dans le génome humain sont vastes, allant de l'adaptation, en passant par les variations phénotypiques, aux prédispositions à certaines maladies. Il a été démontré que les variations du nombre de copies (CNV) influencent le phénotype en modifiant, d'une manière plus ou moins dose-dépendante, l'expression des gènes se situant à l'intérieur de ces régions, mais également celle des gènes se trouvant dans les régions flanquantes. Afin d'étudier les mécanismes possibles sous-jacents à cet effet de voisinage, nous avons comparé les états de modification des histones dans des lignées cellulaires dérivées de patients atteints du syndrome de Williams-Beuren, de la duplication de la région Williams-Beuren, du syndrome de Smith-Magenis ou du syndrome de Di- George et d'individus contrôles en utilisant une version haut-débit de la méthode d'immunoprécipitation de la chromatine (ChIP), appelée ChIP-seq. Nous avons suivi la mono-méthylation de la lysine K20 sur l'histone H4 et la tri-méthylation de la lysine K27 sur l'histone H3, marqueurs respectifs de la chromatine ouverte et fermée. En accord avec les changements de niveaux d'expression observés pour de multiples gènes tout le long des chromosomes affectés par les CNVs, nous avons aussi détecté des régions présentant des modifications d'histones entre les échantillons, situées de part et d'autre des régions critiques, jusqu'aux extrémités du chromosome réarrangé. Nous avons aussi évalué les interactions intra-chromosomiques ayant lieu dans ces cellules par l'utilisation de la technique de capture de conformation des chromosomes (4C-seq). Nous avons observé qu'un groupe de gènes flanquants la région critique du syndrome de Williams-Beuren (WBSCR) forment souvent une boucle, constituant un groupe d'interactions privilégiées entre ces gènes et la WBSCR. La délétion de la WBSCR perturbe l'expression de ce groupe de gènes flanquants, mais également les interactions à grande échelle entre eux et la région réarrangée. Nous en concluons que les larges réarrangements génomiques peuvent aboutir à des changements de l'état de la chromatine pouvant s'étendre bien plus loin que la région critique, affectant donc potentiellement l'expression de manière globale et ainsi modifiant le phénotype des patients.

How confident can we be that orthologs are similar, but paralogs differ?

Homologous genes are classified into orthologs and paralogs, depending on whether they arose by speciation or duplication. It is widely assumed that orthologs share similar functions, whereas paralogs are expected to diverge more from each other. But does this assumption hold up on further examination? We present evidence that orthologs and paralogs are not so different in either their evolutionary rates or their mechanisms of divergence. We emphasize the importance of appropriately designed studies to test models of gene evolution between orthologs and between paralogs. Thus, functional change between orthologs might be as common as between paralogs, and future studies should be designed to test the impact of duplication against this alternative model.

Birth and expression evolution of mammalian microRNA genes.

MicroRNAs (miRNAs) are major post-transcriptional regulators of gene expression, yet their origins and functional evolution in mammals remain little understood due to the lack of appropriate comparative data. Using RNA sequencing, we have generated extensive and comparable miRNA data for five organs in six species that represent all main mammalian lineages and birds (the evolutionary outgroup) with the aim to unravel the evolution of mammalian miRNAs. Our analyses reveal an overall expansion of miRNA repertoires in mammals, with threefold accelerated birth rates of miRNA families in placentals and marsupials, facilitated by the de novo emergence of miRNAs in host gene introns. Generally, our analyses suggest a high rate of miRNA family turnover in mammals with many newly emerged miRNA families being lost soon after their formation. Selectively preserved mammalian miRNA families gradually evolved higher expression levels, as well as altered mature sequences and target gene repertoires, and were apparently mainly recruited to exert regulatory functions in nervous tissues. However, miRNAs that originated on the X chromosome evolved high expression levels and potentially diverse functions during spermatogenesis, including meiosis, through selectively driven duplication-divergence processes. Overall, our study thus provides detailed insights into the birth and evolution of mammalian miRNA genes and the associated selective forces.