A male with unilateral microphthalmia reveals a role for TMX3 in eye development.

Anophthalmia and microphthalmia are important birth defects, but their pathogenesis remains incompletely understood. We studied a patient with severe unilateral microphthalmia who had a 2.7 Mb deletion at chromosome 18q22.1 that was inherited from his mother. In-situ hybridization showed that one of the deleted genes, TMX3, was expressed in the retinal neuroepithelium and lens epithelium in the developing murine eye. We re-sequenced TMX3 in 162 patients with anophthalmia or microphthalmia, and found two missense substitutions in unrelated patients: c.116G>A, predicting p.Arg39Gln, in a male with unilateral microphthalmia and retinal coloboma, and c.322G>A, predicting p.Asp108Asn, in a female with unilateral microphthalmia and severe micrognathia. We used two antisense morpholinos targeted against the zebrafish TMX3 orthologue, zgc:110025, to examine the effects of reduced gene expression in eye development. We noted that the morphant larvae resulting from both morpholinos had significantly smaller eye sizes and reduced labeling with islet-1 antibody directed against retinal ganglion cells at 2 days post fertilization. Co-injection of human wild type TMX3 mRNA rescued the small eye phenotype obtained with both morpholinos, whereas co-injection of human TMX3(p.Arg39Gln) mutant mRNA, analogous to the mutation in the patient with microphthalmia and coloboma, did not rescue the small eye phenotype. Our results show that haploinsufficiency for TMX3 results in a small eye phenotype and represents a novel genetic cause of microphthalmia and coloboma. Future experiments to determine if other thioredoxins are important in eye morphogenesis and to clarify the mechanism of function of TMX3 in eye development are warranted.

VAX1 mutation associated with microphthalmia, corpus callosum agenesis and orofacial clefting-the first description of a VAX1 phenotype in humans.

Vax1 and Vax2 have been implicated in eye development and the closure of the choroid fissure in mice and zebrafish. We sequenced the coding exons of VAX1 and VAX2 in 70 patients with anophthalmia/microphthalmia. In VAX1, we observed homozygosity for two successive nucleotide substitutions c.453G>A and c.454C>A, predicting p.Arg152Ser, in a proband of Egyptian origin with microphthalmia, small optic nerves, cleft lip/palate and corpus callosum agenesis. This mutation affects an invariant residue in the homeodomain of VAX1 and was absent from 96 Egyptian controls. It is likely that the mutation results in a loss of function, as the mutation results in a phenotype similar to the Vax1 homozygous null mouse. We did not identify any mutations in VAX2. This is the first description of a phenotype associated with a VAX1 mutation in humans and establishes VAX1 as a new causative gene for anophthalmia/microphthalmia. ©2011 Wiley Periodicals, Inc.

Towards a new combination therapy for tuberculosis with next generation benzothiazinones.

The benzothiazinone lead compound, BTZ043, kills Mycobacterium tuberculosis by inhibiting the essential flavo-enzyme DprE1, decaprenylphosphoryl-beta-D-ribose 2-epimerase. Here, we synthesized a new series of piperazine-containing benzothiazinones (PBTZ) and show that, like BTZ043, the preclinical candidate PBTZ169 binds covalently to DprE1. The crystal structure of the DprE1-PBTZ169 complex reveals formation of a semimercaptal adduct with Cys387 in the active site and explains the irreversible inactivation of the enzyme. Compared to BTZ043, PBTZ169 has improved potency, safety and efficacy in zebrafish and mouse models of tuberculosis (TB). When combined with other TB drugs, PBTZ169 showed additive activity against M. tuberculosis in vitro except with bedaquiline (BDQ) where synergy was observed. A new regimen comprising PBTZ169, BDQ and pyrazinamide was found to be more efficacious than the standard three drug treatment in a murine model of chronic disease. PBTZ169 is thus an attractive drug candidate to treat TB in humans.

A general definition and nomenclature for alternative splicing events

Understanding the molecular mechanisms responsible for the regulation of the transcriptome present in eukaryotic cells isone of the most challenging tasks in the postgenomic era. In this regard, alternative splicing (AS) is a key phenomenoncontributing to the production of different mature transcripts from the same primary RNA sequence. As a plethora ofdifferent transcript forms is available in databases, a first step to uncover the biology that drives AS is to identify thedifferent types of reflected splicing variation. In this work, we present a general definition of the AS event along with anotation system that involves the relative positions of the splice sites. This nomenclature univocally and dynamically assignsa specific ‘‘AS code’’ to every possible pattern of splicing variation. On the basis of this definition and the correspondingcodes, we have developed a computational tool (AStalavista) that automatically characterizes the complete landscape of ASevents in a given transcript annotation of a genome, thus providing a platform to investigate the transcriptome diversityacross genes, chromosomes, and species. Our analysis reveals that a substantial part—in human more than a quarter—ofthe observed splicing variations are ignored in common classification pipelines. We have used AStalavista to investigate andto compare the AS landscape of different reference annotation sets in human and in other metazoan species and found thatproportions of AS events change substantially depending on the annotation protocol, species-specific attributes, andcoding constraints acting on the transcripts. The AStalavista system therefore provides a general framework to conductspecific studies investigating the occurrence, impact, and regulation of AS.

Whole genome sequencing in patients with retinitis pigmentosa reveals pathogenic DNA structural changes and NEK2 as a new disease gene.

We performed whole genome sequencing in 16 unrelated patients with autosomal recessive retinitis pigmentosa (ARRP), a disease characterized by progressive retinal degeneration and caused by mutations in over 50 genes, in search of pathogenic DNA variants. Eight patients were from North America, whereas eight were Japanese, a population for which ARRP seems to have different genetic drivers. Using a specific workflow, we assessed both the coding and noncoding regions of the human genome, including the evaluation of highly polymorphic SNPs, structural and copy number variations, as well as 69 control genomes sequenced by the same procedures. We detected homozygous or compound heterozygous mutations in 7 genes associated with ARRP (USH2A, RDH12, CNGB1, EYS, PDE6B, DFNB31, and CERKL) in eight patients, three Japanese and five Americans. Fourteen of the 16 mutant alleles identified were previously unknown. Among these, there was a 2.3-kb deletion in USH2A and an inverted duplication of ∼446 kb in EYS, which would have likely escaped conventional screening techniques or exome sequencing. Moreover, in another Japanese patient, we identified a homozygous frameshift (p.L206fs), absent in more than 2,500 chromosomes from ethnically matched controls, in the ciliary gene NEK2, encoding a serine/threonine-protein kinase. Inactivation of this gene in zebrafish induced retinal photoreceptor defects that were rescued by human NEK2 mRNA. In addition to identifying a previously undescribed ARRP gene, our study highlights the importance of rare structural DNA variations in Mendelian diseases and advocates the need for screening approaches that transcend the analysis of the coding sequences of the human genome.

The metal ion binding protein Cnnm4 is mutated in rod-cone dystrophy/amelogenesis imperfecta syndrome

Purpose:To identify the gene causing rod-cone dystrophy/amelogenesis imperfecta Methods:Homozygosity mapping was performed using the Affymetrix 50K XbaI array in one family and candidate genes in the linked interval were sequenced with ABI Dye Terminator, vers. 1 in the index patient of 3 families. The identified mutations were screened in normal control individuals. Expression analyses were performed on RNA extracted from the brain, various parts of the eye and teeth; immunostaining was done on mouse eyes and jaw and knock-down experiments were carried out in zebrafish embroys. Results:Sequencing the coding regions of ancient conserved domain protein 4 (CNNM4), a metal ions transporter, revealed a 1-base pair duplication (p.L438fs) in family A, a p.R236Q mutation in family B and a p.L324P in family C. All these mutations were homozygous and involved very conserved amino acids in paralogs and orthologs. Immunostaining and RT-PCR confirmed that CNNM4 was strongly expressed in various parts of the eye and in the teeth. Morpholino experiments in zebrafish showed a loss of ganglion cells at 5 days post fertilization. Conclusions:The rod-cone dystrophy/amelogenesis imperfecta syndrome is caused by mutation in CNNM4 and is due to aberrant metal ion homeostasis.

Comparative modular analysis of gene expression in vertebrate organs.

ABSTRACT: BACKGROUND: The degree of conservation of gene expression between homologous organs largely remains an open question. Several recent studies reported some evidence in favor of such conservation. Most studies compute organs' similarity across all orthologous genes, whereas the expression level of many genes are not informative about organ specificity. RESULTS: Here, we use a modularization algorithm to overcome this limitation through the identification of inter-species co-modules of organs and genes. We identify such co-modules using mouse and human microarray expression data. They are functionally coherent both in terms of genes and of organs from both organisms. We show that a large proportion of genes belonging to the same co-module are orthologous between mouse and human. Moreover, their zebrafish orthologs also tend to be expressed in the corresponding homologous organs. Notable exceptions to the general pattern of conservation are the testis and the olfactory bulb. Interestingly, some co-modules consist of single organs, while others combine several functionally related organs. For instance, amygdala, cerebral cortex, hypothalamus and spinal cord form a clearly discernible unit of expression, both in mouse and human. CONCLUSIONS: Our study provides a new framework for comparative analysis which will be applicable also to other sets of large-scale phenotypic data collected across different species.

Identification of Novel Craniofacial Regulatory Domains Located far Upstream of SOX9 and Disrupted in Pierre Robin Sequence.

Mutations in the coding sequence of SOX9 cause campomelic dysplasia (CD), a disorder of skeletal development associated with 46,XY disorders of sex development (DSDs). Translocations, deletions, and duplications within a ∼2 Mb region upstream of SOX9 can recapitulate the CD-DSD phenotype fully or partially, suggesting the existence of an unusually large cis-regulatory control region. Pierre Robin sequence (PRS) is a craniofacial disorder that is frequently an endophenotype of CD and a locus for isolated PRS at ∼1.2-1.5 Mb upstream of SOX9 has been previously reported. The craniofacial regulatory potential within this locus, and within the greater genomic domain surrounding SOX9, remains poorly defined. We report two novel deletions upstream of SOX9 in families with PRS, allowing refinement of the regions harboring candidate craniofacial regulatory elements. In parallel, ChIP-Seq for p300 binding sites in mouse craniofacial tissue led to the identification of several novel craniofacial enhancers at the SOX9 locus, which were validated in transgenic reporter mice and zebrafish. Notably, some of the functionally validated elements fall within the PRS deletions. These studies suggest that multiple noncoding elements contribute to the craniofacial regulation of SOX9 expression, and that their disruption results in PRS.

The hourglass and the early conservation models--co-existing patterns of developmental constraints in vertebrates.

Developmental constraints have been postulated to limit the space of feasible phenotypes and thus shape animal evolution. These constraints have been suggested to be the strongest during either early or mid-embryogenesis, which corresponds to the early conservation model or the hourglass model, respectively. Conflicting results have been reported, but in recent studies of animal transcriptomes the hourglass model has been favored. Studies usually report descriptive statistics calculated for all genes over all developmental time points. This introduces dependencies between the sets of compared genes and may lead to biased results. Here we overcome this problem using an alternative modular analysis. We used the Iterative Signature Algorithm to identify distinct modules of genes co-expressed specifically in consecutive stages of zebrafish development. We then performed a detailed comparison of several gene properties between modules, allowing for a less biased and more powerful analysis. Notably, our analysis corroborated the hourglass pattern at the regulatory level, with sequences of regulatory regions being most conserved for genes expressed in mid-development but not at the level of gene sequence, age, or expression, in contrast to some previous studies. The early conservation model was supported with gene duplication and birth that were the most rare for genes expressed in early development. Finally, for all gene properties, we observed the least conservation for genes expressed in late development or adult, consistent with both models. Overall, with the modular approach, we showed that different levels of molecular evolution follow different patterns of developmental constraints. Thus both models are valid, but with respect to different genomic features.

Analysis of gene expression patterns in animals

During my PhD, my aim was to provide new tools to increase our capacity to analyse gene expression patterns, and to study on a large-scale basis the evolution of gene expression in animals. Gene expression patterns (when and where a gene is expressed) are a key feature in understanding gene function, notably in development. It appears clear now that the evolution of developmental processes and of phenotypes is shaped both by evolution at the coding sequence level, and at the gene expression level.Studying gene expression evolution in animals, with complex expression patterns over tissues and developmental time, is still challenging. No tools are available to routinely compare expression patterns between different species, with precision, and on a large-scale basis. Studies on gene expression evolution are therefore performed only on small genes datasets, or using imprecise descriptions of expression patterns.The aim of my PhD was thus to develop and use novel bioinformatics resources, to study the evolution of gene expression. To this end, I developed the database Bgee (Base for Gene Expression Evolution). The approach of Bgee is to transform heterogeneous expression data (ESTs, microarrays, and in-situ hybridizations) into present/absent calls, and to annotate them to standard representations of anatomy and development of different species (anatomical ontologies). An extensive mapping between anatomies of species is then developed based on hypothesis of homology. These precise annotations to anatomies, and this extensive mapping between species, are the major assets of Bgee, and have required the involvement of many co-workers over the years. My main personal contribution is the development and the management of both the Bgee database and the web-application.Bgee is now on its ninth release, and includes an important gene expression dataset for 5 species (human, mouse, drosophila, zebrafish, Xenopus), with the most data from mouse, human and zebrafish. Using these three species, I have conducted an analysis of gene expression evolution after duplication in vertebrates.Gene duplication is thought to be a major source of novelty in evolution, and to participate to speciation. It has been suggested that the evolution of gene expression patterns might participate in the retention of duplicate genes. I performed a large-scale comparison of expression patterns of hundreds of duplicated genes to their singleton ortholog in an outgroup, including both small and large-scale duplicates, in three vertebrate species (human, mouse and zebrafish), and using highly accurate descriptions of expression patterns. My results showed unexpectedly high rates of de novo acquisition of expression domains after duplication (neofunctionalization), at least as high or higher than rates of partitioning of expression domains (subfunctionalization). I found differences in the evolution of expression of small- and large-scale duplicates, with small-scale duplicates more prone to neofunctionalization. Duplicates with neofunctionalization seemed to evolve under more relaxed selective pressure on the coding sequence. Finally, even with abundant and precise expression data, the majority fate I recovered was neither neo- nor subfunctionalization of expression domains, suggesting a major role for other mechanisms in duplicate gene retention.

Genome-wide association and functional follow-up reveals new loci for kidney function.

Chronic kidney disease (CKD) is an important public health problem with a genetic component. We performed genome-wide association studies in up to 130,600 European ancestry participants overall, and stratified for key CKD risk factors. We uncovered 6 new loci in association with estimated glomerular filtration rate (eGFR), the primary clinical measure of CKD, in or near MPPED2, DDX1, SLC47A1, CDK12, CASP9, and INO80. Morpholino knockdown of mpped2 and casp9 in zebrafish embryos revealed podocyte and tubular abnormalities with altered dextran clearance, suggesting a role for these genes in renal function. By providing new insights into genes that regulate renal function, these results could further our understanding of the pathogenesis of CKD.

BRF1 mutations alter RNA polymerase III-dependent transcription and cause neurodevelopmental anomalies.

RNA polymerase III (Pol III) synthesizes tRNAs and other small noncoding RNAs to regulate protein synthesis. Dysregulation of Pol III transcription has been linked to cancer, and germline mutations in genes encoding Pol III subunits or tRNA processing factors cause neurogenetic disorders in humans, such as hypomyelinating leukodystrophies and pontocerebellar hypoplasia. Here we describe an autosomal recessive disorder characterized by cerebellar hypoplasia and intellectual disability, as well as facial dysmorphic features, short stature, microcephaly, and dental anomalies. Whole-exome sequencing revealed biallelic missense alterations of BRF1 in three families. In support of the pathogenic potential of the discovered alleles, suppression or CRISPR-mediated deletion of brf1 in zebrafish embryos recapitulated key neurodevelopmental phenotypes; in vivo complementation showed all four candidate mutations to be pathogenic in an apparent isoform-specific context. BRF1 associates with BDP1 and TBP to form the transcription factor IIIB (TFIIIB), which recruits Pol III to target genes. We show that disease-causing mutations reduce Brf1 occupancy at tRNA target genes in Saccharomyces cerevisiae and impair cell growth. Moreover, BRF1 mutations reduce Pol III-related transcription activity in vitro. Taken together, our data show that BRF1 mutations that reduce protein activity cause neurodevelopmental anomalies, suggesting that BRF1-mediated Pol III transcription is required for normal cerebellar and cognitive development.

Comparative modular analysis of gene expression in vertebrate development

The focus of my PhD research was the concept of modularity. In the last 15 years, modularity has become a classic term in different fields of biology. On the conceptual level, a module is a set of interacting elements that remain mostly independent from the elements outside of the module. I used modular analysis techniques to study gene expression evolution in vertebrates. In particular, I identified ``natural'' modules of gene expression in mouse and human, and I showed that expression of organ-specific and system-specific genes tends to be conserved between such distance vertebrates as mammals and fishes. Also with a modular approach, I studied patterns of developmental constraints on transcriptome evolution. I showed that none of the two commonly accepted models of the evolution of embryonic development (``evo-devo'') are exclusively valid. In particular, I found that the conservation of the sequences of regulatory regions is highest during mid-development of zebrafish, and thus it supports the ``hourglass model''. In contrast, events of gene duplication and new gene introduction are most rare in early development, which supports the ``early conservation model''. In addition to the biological insights on transcriptome evolution, I have also discussed in detail the advantages of modular approaches in large-scale data analysis. Moreover, I re-analyzed several studies (published in high-ranking journals), and showed that their conclusions do not hold out under a detailed analysis. This demonstrates that complex analysis of high-throughput data requires a co-operation between biologists, bioinformaticians, and statisticians.

Dissecting the role of low-complexity regions in the evolution of vertebrate proteins

Low-complexity regions (LCRs) in proteins are tracts that are highly enriched in one or a few aminoacids. Given their high abundance, and their capacity to expand in relatively short periods of time through replication slippage, they can greatly contribute to increase protein sequence space and generate novel protein functions. However, little is known about the global impact of LCRs on protein evolution. We have traced back the evolutionary history of 2,802 LCRs from a large set of homologous protein families from H.sapiens, M.musculus, G.gallus, D.rerio and C.intestinalis. Transcriptional factors and other regulatory functions are overrepresented in proteins containing LCRs. We have found that the gain of novel LCRs is frequently associated with repeat expansion whereas the loss of LCRs is more often due to accumulation of amino acid substitutions as opposed to deletions. This dichotomy results in net protein sequence gain over time. We have detected a significant increase in the rate of accumulation of novel LCRs in the ancestral Amniota and mammalian branches, and a reduction in the chicken branch. Alanine and/or glycine-rich LCRs are overrepresented in recently emerged LCR sets from all branches, suggesting that their expansion is better tolerated than for other LCR types. LCRs enriched in positively charged amino acids show the contrary pattern, indicating an important effect of purifying selection in their maintenance. We have performed the first large-scale study on the evolutionary dynamics of LCRs in protein families. The study has shown that the composition of an LCR is an important determinant of its evolutionary pattern.

Study of inner ear and lateral line hair cell regeneration

Death of sensory hair cells in the inner ear results in two global health problems that millions of people around the world suffer: hearing loss and balance disorders. Hair cells convert sound vibrations and head movements into electrical signals that are conveyed to the brain, and as a result of aging, exposure to noise, modern drugs or genetic predisposition, hair cells die. In mammals, the great majority of hair cells are produced during embryogenesis, and hair cells that are lost after birth are not replaceable. However, in the last decades, researches have shown some model organisms that retain the ability to regenerate hair cells damaged after embryogenesis, such as Zebrafish and chicken, providing clues as to the cellular and molecular mechanisms that may block hair cell regeneration in mammals. This discovery initiated a search for methods to stimulate regeneration or replacement of hair cells in mammals, a search that, if fruitful, will revolutionize the treatment of hearing loss and balance disorders. One aim of my project is to study the role of retinoic acid in adult Zebrafish and in mice, which is a metabolite of vitamin A known as an essential molecule to activate hair cell regeneration after cells damaged in Zebrafish embryo. We want to study important genes involved in retinoic acid pathway, such as Aldh1a3 and RARs genes, to check what their role is in the inner ear of adult Zebrafish and compare result obtained in the inner ear of mice. On the other hand, Zebrafish lateral line contains neuromast, which are formed by the same structure than the inner ear: hair cells surrounded by supporting cells and neurons. The lateral line is a structure below the skin's surface that makes easier to damage hair cells to study their regeneration. For that reason, another aim of my project is to study how Sox2 and Atoh1, essential genes during the inner ear development, change their expression during hair cell regeneration in the lateral line. In my project, the most important concepts related to Zebrafish world are explained in order to understand why we have studied this animal and these essential genes. Then, techniques that we used are explained, with their protocol attached in the annexes. Finally, results of my project are shown, but many of them were not expected and they would be needed to follow studying.