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.

Interleukin-1 receptor antagonist gene (IL-1RN) polymorphism is a predictive factor of clinical pregnancy after IVF.

BACKGROUND: Only 25% of IVF transfer cycles lead to a clinical pregnancy, calling for continued technical progress but also more in depth analysis of patients' individual characteristics. The interleukin-1 (IL-1) system and matrix metalloproteinases (MMPs) are strongly implicated in embryo implantation. The genes coding for IL-1Ra (gene symbol IL-1RN), IL-1beta, MMP2 and MMP9 bear functional polymorphisms. We analysed the maternal genetic profile at these polymorphic sites in IVF patients, to determine possible correlations with IVF outcome. METHODS: One hundred and sixty women undergoing an IVF cycle were enrolled and a buccal smear was obtained. The presence of IL-1RN variable number of tandem repeats and IL-1B + 3953, MMP2-1306 and MMP9-1562 single nucleotide substitutions were determined. Patients were divided into pregnancy failures (119), biochemical pregnancies (8) and clinical pregnancies (33). RESULTS: There was a 40% decrease in IL-1RN*2 allele frequency (P = 0.024) and a 45% decrease in IL-1RN*2 carrier status in the clinical pregnancy group as compared to the pregnancy failure group (P = 0.017). This decrease was still statistically significant after a multivariate logistic regression analysis. The likelihood of a clinical pregnancy was decreased accordingly in IL-1RN*2 carriers: odds ratio = 0.349, 95% confidence interval = 0.2-0.8, P = 0.017. The IL-1B, MMP2 and MMP9 polymorphisms showed no correlation with IVF outcome. CONCLUSIONS: IL-1RN*2 allele carriage is associated with a poor prognosis of achieving a pregnancy after IVF.

Myc's other life: stem cells and beyond.

Over the last three decades genetic and biochemical studies have revealed the pleiotropic effects of the Myc oncoprotein. While cell line studies have defined the intracellular processes regulated by Myc such as proliferation, differentiation, and metabolic growth, in vivo studies have confirmed these functions, and revealed roles in acquisition and maintenance of stem cell properties. These roles may be partially mediated by Myc's capacity to modify the chromatin landscape on a global scale. Myc also regulates numerous protein-coding transcripts, and many noncoding RNAs (rRNAs, tRNAs, and miRNAs). As Myc activity directly correlates with protein expression, further complexity is provided by post-translational modifications that regulate Myc in normal stem cells or deregulate it in malignant stem cells.

Avoiding epigenetic variations in transgene expression with genetic insulators

Gene transfer that relies on integrating vectors often suffers from epigenetic or regulatory effects that influence the expression of the therapeutic gene and=or of cellular genes located near the vector integration site in the chromosome. Insulator elements act to block gene activation by enhancers, while chromatin domain boundary or barrier sequences prevent gene-silencing effects. At present, the modes of action of insulator and barriers are poorly understood, and their use in the context of gene therapies remains to be documented. Using combinations of reporter genes coding for indicator fluorescent proteins, we constructed assay systems that allow the quantification of the insulator or of the barrier activities of genetic elements in individual cells. This presentation will illustrate how these assay systems were used to identify short DNA elements that can insulate nearby genes from activation by viral vector enhancer elements, and=or that can block the propagation of a silent chromatin structure that leads to gene silencing. We will show that small elements of the order of 100-400 nucleotides can be designed to achieve both insulator and boundary function, as needed for safer integrating viral vectors.

Origins, evolution, and phenotypic impact of new genes.

Ever since the pre-molecular era, the birth of new genes with novel functions has been considered to be a major contributor to adaptive evolutionary innovation. Here, I review the origin and evolution of new genes and their functions in eukaryotes, an area of research that has made rapid progress in the past decade thanks to the genomics revolution. Indeed, recent work has provided initial whole-genome views of the different types of new genes for a large number of different organisms. The array of mechanisms underlying the origin of new genes is compelling, extending way beyond the traditionally well-studied source of gene duplication. Thus, it was shown that novel genes also regularly arose from messenger RNAs of ancestral genes, protein-coding genes metamorphosed into new RNA genes, genomic parasites were co-opted as new genes, and that both protein and RNA genes were composed from scratch (i.e., from previously nonfunctional sequences). These mechanisms then also contributed to the formation of numerous novel chimeric gene structures. Detailed functional investigations uncovered different evolutionary pathways that led to the emergence of novel functions from these newly minted sequences and, with respect to animals, attributed a potentially important role to one specific tissue--the testis--in the process of gene birth. Remarkably, these studies also demonstrated that novel genes of the various types significantly impacted the evolution of cellular, physiological, morphological, behavioral, and reproductive phenotypic traits. Consequently, it is now firmly established that new genes have indeed been major contributors to the origin of adaptive evolutionary novelties.

Germ cell-specific targeting of DICER or DGCR8 reveals a novel role for endo-siRNAs in the progression of mammalian spermatogenesis and male fertility.

Small non-coding RNAs act as critical regulators of gene expression and are essential for male germ cell development and spermatogenesis. Previously, we showed that germ cell-specific inactivation of Dicer1, an endonuclease essential for the biogenesis of micro-RNAs (miRNAs) and endogenous small interfering RNAs (endo-siRNAs), led to complete male infertility due to alterations in meiotic progression, increased spermatocyte apoptosis and defects in the maturation of spermatozoa. To dissect the distinct physiological roles of miRNAs and endo-siRNAs in spermatogenesis, we compared the testicular phenotype of mice with Dicer1 or Dgcr8 depletion in male germ cells. Dgcr8 mutant mice, which have a defective miRNA pathway while retaining an intact endo-siRNA pathway, were also infertile and displayed similar defects, although less severe, to Dicer1 mutant mice. These included cumulative defects in meiotic and haploid phases of spermatogenesis, resulting in oligo-, terato-, and azoospermia. In addition, we found by RNA sequencing of purified spermatocytes that inactivation of Dicer1 and the resulting absence of miRNAs affected the fine tuning of protein-coding gene expression by increasing low level gene expression. Overall, these results emphasize the essential role of miRNAs in the progression of spermatogenesis, but also indicate a role for endo-siRNAs in this process.

Reduced choroidal neovascularization by AAV-anti-VEGF shRNA delivery

VEGF plays an essential role in ocular angiogenic diseases including the late-stage form of AMD, the primary cause of vision loss in the western world. Over-expression of VEGF leads to development of vasculature emanating from the choroid, invading the subretinal space through breaks in Bruch's membrane. Strategies leading to long-term suppression of inappropriate ocular angiogenesis are required. A panel of 10 shRNAs targeting the coding region of human VEGF165 was tested in HEK293 cells and in the human retinal pigment epithelial cell line, ARPE-19. VEGF knock-down up to 92% was achieved by co-transfecting shRNAexpressing constructs with plasmid encoding the Renilla luciferase gene fused to the VEGF165 sequence. For in vivo delivery of the most potent shRNA cassette, both single-stranded and self-complementary rAAV vectors were packaged in serotype 8 capsids. Intramuscular administration in mice led to localized expression and 96% knock-down of endogenous VEGF. Using eGFP as a marker, efficient gene transfer of retinal pigment epithelial cells, the cells thought to be responsible for the abnormal VEGF production, was obtained by subretinal delivery of rAAV2.8 vectors. The capacity of rAAV-encoded shRNAs to silence endogenous VEGF gene expression was evaluated in the laser-induced murine model of choroidal neovascularization (CNV). In this mouse model of AMD, sizes of the CNV were found to be significantly reduced following rAAV-shRNA subretinal delivery. Thus, our results indicate that gene transfer combining AAV-mediated delivery with triggering of the endogenous RNAi pathway can be used for anti-VEGF therapy and holds great promise for the treatment of AMD.

Efavirenz intoxication due to a new CYP2B6 constellation.

Here, we describe severe neuropsychiatric symptoms in an HIV-positive Asian man with extremely high efavirenz plasma levels while receiving standard treatment with efavirenz/tenofovir/emtricitabine fixed-dose regimen. Genetic examination revealed compound homozygosity for loss-of-function alleles of CYP2B6, including coding for a rare truncated protein. Neuropsychiatric symptoms resolved completely after efavirenz discontinuation.

ICD-11 for quality and safety: overview of the who quality and safety topic advisory group.

This paper outlines the approach that the WHO's Family of International Classifications (WHO-FIC) network is undertaking to create ICD-11. We also outline the more focused work of the Quality and Safety Topic Advisory Group, whose activities include the following: (i) cataloguing existing ICD-9 and ICD-10 quality and safety indicators; (ii) reviewing ICD morbidity coding rules for main condition, diagnosis timing, numbers of diagnosis fields and diagnosis clustering; (iii) substantial restructuring of the health-care related injury concepts coded in the ICD-10 chapters 19/20, (iv) mapping of ICD-11 quality and safety concepts to the information model of the WHO's International Classification for Patient Safety and the AHRQ Common Formats; (v) the review of vertical chapter content in all chapters of the ICD-11 beta version and (vi) downstream field testing of ICD-11 prior to its official 2015 release. The transition from ICD-10 to ICD-11 promises to produce an enhanced classification that will have better potential to capture important concepts relevant to measuring health system safety and quality-an important use case for the classification.

Mutations in the LHX2 gene are not a frequent cause of micro/anophthalmia.

Purpose: Microphthalmia and anophthalmia are at the severe end of the spectrum of abnormalities in ocular development. A few genes (orthodenticle homeobox 2 [OTX2], retina and anterior neural fold homeobox [RAX], SRY-box 2 [SOX2], CEH10 homeodomain-containing homolog [CHX10], and growth differentiation factor 6 [GDF6]) have been implicated mainly in isolated micro/anophthalmia but causative mutations of these genes explain less than a quarter of these developmental defects. The essential role of the LIM homeobox 2 (LHX2) transcription factor in early eye development has recently been documented. We postulated that mutations in this gene could lead to micro/anophthalmia, and thus performed molecular screening of its sequence in patients having micro/anophthalmia. Methods: Seventy patients having non-syndromic forms of colobomatous microphthalmia (n=25), isolated microphthalmia (n=18), or anophthalmia (n=17), and syndromic forms of micro/anophthalmia (n=10) were included in this study after negative molecular screening for OTX2, RAX, SOX2, and CHX10 mutations. Mutation screening of LHX2 was performed by direct sequencing of the coding sequences and intron/exon boundaries. Results: Two heterozygous variants of unknown significance (c.128C > G [p.Pro43Arg]; c.776C > A [p.Pro259Gln]) were identified in LHX2 among the 70 patients. These variations were not identified in a panel of 100 control patients of mixed origins. The variation c.776C > A (p.Pro259Gln) was considered as non pathogenic by in silico analysis, while the variation c.128C > G (p.Pro43Arg) considered as deleterious by in silico analysis and was inherited from the asymptomatic father. Conclusions: Mutations in LHX2 do not represent a frequent cause of micro/anophthalmia.

Expression of genes encoding cytotoxic cell-associated serine proteases in thymocytes.

A family of homologous serine esterases designated granzyme A-H and the pore-forming protein perforin are present in cytoplasmic granules of mature peripheral cytolytic T lymphocytes and natural killer cells. In vivo, the majority of cytotoxic T cells containing these granule-associated proteins are of the CD4-CD8+ phenotype. It is generally assumed that these cells are derived from immature CD4-CD8- thymocytes. However, the precise intrathymic differentiation steps leading to functionally mature cytotoxic T cells are unclear. Thus we decided to analyze the expression of genes in the thymus which are preferentially expressed in mature cytotoxic cells, i.e. granzyme A, granzyme B, and perforin. In situ hybridization on tissue sections revealed the expression of genes coding for granzyme A and granzyme B in the thymus. No evidence was found, however, for thymocytes expressing the perforin gene. Granzyme A and granzyme B mRNA positive cells in the thymus are almost exclusively CD4-CD8- thymocytes, particularly of the CD3- IL2R- phenotype.

Synergistic transcriptional activation by CTF/NF-I and the estrogen receptor involves stabilized interactions with a limiting target factor.

Transcription initiation at eukaryotic protein-coding gene promoters is regulated by a complex interplay of site-specific DNA-binding proteins acting synergistically or antagonistically. Here, we have analyzed the mechanisms of synergistic transcriptional activation between members of the CCAAT-binding transcription factor/nuclear factor I (CTF/NF-I) family and the estrogen receptor. By using cotransfection experiments with HeLa cells, we show that the proline-rich transcriptional activation domain of CTF-1, when fused to the GAL4 DNA-binding domain, synergizes with each of the two estrogen receptor-activating regions. Cooperative DNA binding between the GAL4-CTF-1 fusion and the estrogen receptor does not occur in vitro, and in vivo competition experiments demonstrate that both activators can be specifically inhibited by the overexpression of a proline-rich competitor, indicating that a common limiting factor is mediating their transcriptional activation functions. Furthermore, the two activators functioning synergistically are much more resistant to competition than either factor alone, suggesting that synergism between CTF-1 and the estrogen receptor is the result of a stronger tethering of the limiting target factor(s) to the two promoter-bound activators.

trome, trEST and trGEN: databases of predicted protein sequences.

We previously introduced two new protein databases (trEST and trGEN) of hypothetical protein sequences predicted from EST and HTG sequences, respectively. Here, we present the updates made on these two databases plus a new database (trome), which uses alignments of EST data to HTG or full genomes to generate virtual transcripts and coding sequences. This new database is of higher quality and since it contains the information in a much denser format it is of much smaller size. These new databases are in a Swiss-Prot-like format and are updated on a weekly basis (trEST and trGEN) or every 3 months (trome). They can be downloaded by anonymous ftp from ftp://ftp.isrec.isb-sib.ch/pub/databases.

No evidence for mutations of CTCFL/BORIS in Silver-Russell syndrome patients with IGF2/H19 imprinting control region 1 hypomethylation.

BACKGROUND: Silver-Russell syndrome (SRS) is a genetically and clinically heterogeneous disease. Although no protein coding gene defects have been reported in SRS patients, approximately 50% of SRS patients carry epimutations (hypomethylation) at the IGF2/H19 imprinting control region 1 (ICR1). Proper methylation at ICR1 is crucial for the imprinted expression of IGF2, a fetal growth factor. CTCFL, a testis-specific protein, has recently been proposed to play a role in the establishment of DNA methylation at the murine equivalent of ICR1. A screen was undertaken to assess whether CTCFL is mutated in SRS patients with hypomethylation, to explore a link between the observed epimutations and a genetic cause of the disease. METHODOLOGY/PRINCIPAL FINDINGS: DNA was obtained from 36 SRS patients with hypomethylation at ICR1. All CTCFL coding exons were sequenced and analyzed for duplications/deletions using both multiplex ligation-dependent probe amplification, with a custom CTCFL probe set, and genomic qPCR. Novel SNP alleles were analyzed for potential differential splicing in vitro utilizing a splicing assay. Neither mutations of CTCFL nor duplications/deletions were observed. Five novel SNPs were identified and have been submitted to dbSNP. In silico splice prediction suggested one novel SNP, IVS2-66A>C, activated a cryptic splice site, resulting in aberrant splicing and premature termination. In vitro splicing assays did not confirm predicted aberrant splicing. CONCLUSIONS/SIGNIFICANCE: As no mutations were detected at CTCFL in the patients examined, we conclude that genetic alterations of CTCFL are not responsible for the SRS hypomethylation. We suggest that analysis of other genes involved in the establishment of DNA methylation at imprinted genes, such as DNMT3A and DNMT3L, may provide insight into the genetic cause of hypomethylation in SRS patients.

An integrated encyclopedia of DNA elements in the human genome.

The human genome encodes the blueprint of life, but the function of the vast majority of its nearly three billion bases is unknown. The Encyclopedia of DNA Elements (ENCODE) project has systematically mapped regions of transcription, transcription factor association, chromatin structure and histone modification. These data enabled us to assign biochemical functions for 80% of the genome, in particular outside of the well-studied protein-coding regions. Many discovered candidate regulatory elements are physically associated with one another and with expressed genes, providing new insights into the mechanisms of gene regulation. The newly identified elements also show a statistical correspondence to sequence variants linked to human disease, and can thereby guide interpretation of this variation. Overall, the project provides new insights into the organization and regulation of our genes and genome, and is an expansive resource of functional annotations for biomedical research.