Evidence for transcript networks composed of chimeric RNAs in human cells.

The classic organization of a gene structure has followed the Jacob and Monod bacterial gene model proposed more than 50 years ago. Since then, empirical determinations of the complexity of the transcriptomes found in yeast to human has blurred the definition and physical boundaries of genes. Using multiple analysis approaches we have characterized individual gene boundaries mapping on human chromosomes 21 and 22. Analyses of the locations of the 5' and 3' transcriptional termini of 492 protein coding genes revealed that for 85% of these genes the boundaries extend beyond the current annotated termini, most often connecting with exons of transcripts from other well annotated genes. The biological and evolutionary importance of these chimeric transcripts is underscored by (1) the non-random interconnections of genes involved, (2) the greater phylogenetic depth of the genes involved in many chimeric interactions, (3) the coordination of the expression of connected genes and (4) the close in vivo and three dimensional proximity of the genomic regions being transcribed and contributing to parts of the chimeric RNAs. The non-random nature of the connection of the genes involved suggest that chimeric transcripts should not be studied in isolation, but together, as an RNA network.

Hominoid fission of chromosome 14/15 and the role of segmental duplications.

Ape chromosomes homologous to human chromosomes 14 and 15 were generated by a fission event of an ancestral submetacentric chromosome, where the two chromosomes were joined head-to-tail. The hominoid ancestral chromosome most closely resembles the macaque chromosome 7. In this work, we provide insights into the evolution of human chromosomes 14 and 15, performing a comparative study between macaque boundary region 14/15 and the orthologous human regions. We construct a 1.6-Mb contig of macaque BAC clones in the region orthologous to the ancestral hominoid fission site and use it to define the structural changes that occurred on human 14q pericentromeric and 15q subtelomeric regions. We characterize the novel euchromatin-heterochromatin transition region (∼20 Mb) acquired during the neocentromere establishment on chromosome 14, and find it was mainly derived through pericentromeric duplications from ancestral hominoid chromosomes homologous to human 2q14-qter and 10. Further, we show a relationship between evolutionary hotspots and low-copy repeat loci for chromosome 15, revealing a possible role of segmental duplications not only in mediating but also in "stitching" together rearrangement breakpoints.

Genome-wide analysis of cancer/testis gene expression.

Cancer/Testis (CT) genes, normally expressed in germ line cells but also activated in a wide range of cancer types, often encode antigens that are immunogenic in cancer patients, and present potential for use as biomarkers and targets for immunotherapy. Using multiple in silico gene expression analysis technologies, including twice the number of expressed sequence tags used in previous studies, we have performed a comprehensive genome-wide survey of expression for a set of 153 previously described CT genes in normal and cancer expression libraries. We find that although they are generally highly expressed in testis, these genes exhibit heterogeneous gene expression profiles, allowing their classification into testis-restricted (39), testis/brain-restricted (14), and a testis-selective (85) group of genes that show additional expression in somatic tissues. The chromosomal distribution of these genes confirmed the previously observed dominance of X chromosome location, with CT-X genes being significantly more testis-restricted than non-X CT. Applying this core classification in a genome-wide survey we identified >30 CT candidate genes; 3 of them, PEPP-2, OTOA, and AKAP4, were confirmed as testis-restricted or testis-selective using RT-PCR, with variable expression frequencies observed in a panel of cancer cell lines. Our classification provides an objective ranking for potential CT genes, which is useful in guiding further identification and characterization of these potentially important diagnostic and therapeutic targets.

The TPTE gene family: cellular expression, subcellular localization and alternative splicing.

The human TPTE (Transmembrane Phosphatase with TEnsin homology) gene family encodes a PTEN-related tyrosine phosphatase with four potential transmembrane domains. Chromosomal mapping revealed multiple copies of the TPTE gene on chromosomes 13, 15, 21, 22 and Y. Human chromosomes 13 and 21 copies encode two functional proteins, TPIP (TPTE and PTEN homologous Inositol lipid Phosphatase) and TPTE, respectively, whereas only one copy of the gene exists in the mouse genome. In the present study, we show that TPTE and TPIP proteins are expressed in secondary spermatocytes and/or prespermatids. In addition, we report the existence of several novel alternatively spliced isoforms of these two proteins with variable number of transmembrane domains. The latter has no influence on the subcellular localization of these different peptides as shown by co-immunofluorescence experiments. Finally, we identify another expressed TPTE copy, mapping to human chromosome 22, whose transcription appears to be under the control of the LTR of human endogenous retrovirus RTVL-H3.

Extensive gene traffic on the mammalian X chromosome.

Mammalian sex chromosomes have undergone profound changes since evolving from ancestral autosomes. By examining retroposed genes in the human and mouse genomes, we demonstrate that, during evolution, the mammalian X chromosome has generated and recruited a disproportionately high number of functional retroposed genes, whereas the autosomes experienced lower gene turnover. Most autosomal copies originating from X-linked genes exhibited testis-biased expression. Such export is incompatible with mutational bias and is likely driven by natural selection to attain male germline function. However, the excess recruitment is consistent with a combination of both natural selection and mutational bias.

DNA methylation profiles of human active and inactive X chromosomes.

X-chromosome inactivation (XCI) is a dosage compensation mechanism that silences the majority of genes on one X chromosome in each female cell. To characterize epigenetic changes that accompany this process, we measured DNA methylation levels in 45,X patients carrying a single active X chromosome (X(a)), and in normal females, who carry one X(a) and one inactive X (X(i)). Methylated DNA was immunoprecipitated and hybridized to high-density oligonucleotide arrays covering the X chromosome, generating epigenetic profiles of active and inactive X chromosomes. We observed that XCI is accompanied by changes in DNA methylation specifically at CpG islands (CGIs). While the majority of CGIs show increased methylation levels on the X(i), XCI actually results in significant reductions in methylation at 7% of CGIs. Both intra- and inter-genic CGIs undergo epigenetic modification, with the biggest increase in methylation occurring at the promoters of genes silenced by XCI. In contrast, genes escaping XCI generally have low levels of promoter methylation, while genes that show inter-individual variation in silencing show intermediate increases in methylation. Thus, promoter methylation and susceptibility to XCI are correlated. We also observed a global correlation between CGI methylation and the evolutionary age of X-chromosome strata, and that genes escaping XCI show increased methylation within gene bodies. We used our epigenetic map to predict 26 novel genes escaping XCI, and searched for parent-of-origin-specific methylation differences, but found no evidence to support imprinting on the human X chromosome. Our study provides a detailed analysis of the epigenetic profile of active and inactive X chromosomes.

Annotating the human genome

2 Abstract2.1 En françaisLe séquençage du génome humain est un pré-requis fondamental à la compréhension de la biologie de l'être humain. Ce projet achevé, les scientifiques ont dû faire face à une tâche aussi importante, comprendre cette suite de 3 milliards de lettres qui compose notre génome. Le consortium ENCODE (ENCyclopedia Of Dna Elements) fût formé comme une suite logique au projet du génome humain. Son rôle est d'identifier tous les éléments fonctionnels de notre génome incluant les régions transcrites, les sites d'attachement des facteurs de transcription, les sites hypersensibles à la DNAse I ainsi que les marqueurs de modification des histones. Dans le cadre de ma thèse doctorale, j'ai participé à 2 sous-projets d'ENCODE. En premier lieu, j'ai eu la tâche de développer et d'optimiser une technique de validation expérimentale à haut rendement de modèles de gènes qui m'a permis d'estimer la qualité de la plus récente annotation manuelle. Ce nouveau processus de validation est bien plus efficace que la technique RNAseq qui est actuellement en train de devenir la norme. Cette technique basée sur la RT-PCR, m'a notamment permis de découvrir de nouveaux exons dans 10% des régions interrogées. En second lieu j'ai participé à une étude ayant pour but d'identifier les extrémités de tous les gènes des chromosomes humains 21 et 22. Cette étude à permis l'identification à large échelle de transcrits chimères comportant des séquences provenant de deux gènes distincts pouvant être à une grande distance l'un de autre.2.2 In EnglishThe completion of the human genome sequence js the prerequisite to fully understand the biology of human beings. This project achieved, scientists had to face another challenging task, understanding the meaning of the 3 billion letters composing this genome. As a logical continuation of the human genome project, the ENCODE (ENCyclopedia Of DNA Elements) consortium was formed with the aim of annotating all its functional elements. These elements include transcribed regions, transcription binding sites, DNAse I hypersensitive sites and histone modification marks. In the frame of my PhD thesis, I was involved in two sub-projects of ENCODE. Firstly I developed and optimized an high throughput method to validate gene models, which allowed me to assess the quality of the most recent manually-curated annotation. This novel experimental validation pipeline is extremely effective, far more so than transcriptome profiling through RNA sequencing, which is becoming the norm. This RT-PCR-seq targeted-approach is likewise particularly efficient in identifying novel exons, as we discovered about 10% of loci with unannotated exons. Secondly, I participated to a study aiming to identify the gene boundaries of all genes in the human chromosome 21 and 22. This study led to the identification of chimeric transcripts that are composed of sequences coming form two distinct genes that can be map far away from each other.

Nineteen additional unpredicted transcripts from human chromosome 21.

The identification of all human chromosome 21 (HC21) genes is a necessary step in understanding the molecular pathogenesis of trisomy 21 (Down syndrome). The first analysis of the sequence of 21q included 127 previously characterized genes and predicted an additional 98 novel anonymous genes. Recently we evaluated the quality of this annotation by characterizing a set of HC21 open reading frames (C21orfs) identified by mapping spliced expressed sequence tags (ESTs) and predicted genes (PREDs), identified only in silico. This study underscored the limitations of in silico-only gene prediction, as many PREDs were incorrectly predicted. To refine the HC21 annotation, we have developed a reliable algorithm to extract and stringently map sequences that contain bona fide 3' transcript ends to the genome. We then created a specific 21q graphical display allowing an integrated view of the data that incorporates new ESTs as well as features such as CpG islands, repeats, and gene predictions. Using these tools we identified 27 new putative genes. To validate these, we sequenced previously cloned cDNAs and carried out RT-PCR, 5'- and 3'-RACE procedures, and comparative mapping. These approaches substantiated 19 new transcripts, thus increasing the HC21 gene count by 9.5%. These transcripts were likely not previously identified because they are small and encode small proteins. We also identified four transcriptional units that are spliced but contain no obvious open reading frame. The HC21 data presented here further emphasize that current gene prediction algorithms miss a substantial number of transcripts that nevertheless can be identified using a combination of experimental approaches and multiple refined algorithms.

Gene expression variation and expression quantitative trait mapping of human chromosome 21 genes.

Inter-individual differences in gene expression are likely to account for an important fraction of phenotypic differences, including susceptibility to common disorders. Recent studies have shown extensive variation in gene expression levels in humans and other organisms, and that a fraction of this variation is under genetic control. We investigated the patterns of gene expression variation in a 25 Mb region of human chromosome 21, which has been associated with many Down syndrome (DS) phenotypes. Taqman real-time PCR was used to measure expression variation of 41 genes in lymphoblastoid cells of 40 unrelated individuals. For 25 genes found to be differentially expressed, additional analysis was performed in 10 CEPH families to determine heritabilities and map loci harboring regulatory variation. Seventy-six percent of the differentially expressed genes had significant heritabilities, and genomewide linkage analysis led to the identification of significant eQTLs for nine genes. Most eQTLs were in trans, with the best result (P=7.46 x 10(-8)) obtained for TMEM1 on chromosome 12q24.33. A cis-eQTL identified for CCT8 was validated by performing an association study in 60 individuals from the HapMap project. SNP rs965951 located within CCT8 was found to be significantly associated with its expression levels (P=2.5 x 10(-5)) confirming cis-regulatory variation. The results of our study provide a representative view of expression variation of chromosome 21 genes, identify loci involved in their regulation and suggest that genes, for which expression differences are significantly larger than 1.5-fold in control samples, are unlikely to be involved in DS-phenotypes present in all affected individuals.

Genomic organization, fine-mapping, and expression of the human islet-brain 1 (IB1)/c-Jun-amino-terminal kinase interacting protein-1 (JIP-1) gene.

Islet-brain 1 (IB1), a regulator of the pancreatic beta-cell function in the rat, is homologous to JIP-1, a murine inhibitor of c-Jun amino-terminal kinase (JNK). Whether IB1 and JIP-1 are present in humans was not known. We report the sequence of the 2133-bp human IB1 cDNA, the expression, structure, and fine-mapping of the human IB1 gene, and the characterization of an IB1 pseudogene. Human IB1 is 94% identical to rat IB1. The tissue-specific expression of IB1 in human is similar to that observed in rodent. The IB1 gene contains 12 exons and maps to chromosome 11 (11p11.2-p12), a region that is deleted in DEFECT-11 syndrome. Apart from an IB1 pseudogene on chromosome 17 (17q21), no additional IB1-related gene was found in the human genome. Our data indicate that the sequence and expression pattern of IB1 are highly conserved between rodent and human and provide the necessary tools to investigate whether IB1 is involved in human diseases.

Biochemical characterization of recombinant human telomerase

Résumé Les télomères sont les structures ADN-protéines des extrémités des chromosomes des eucaryotes. L'ADN télomérique est constitué de courtes séquences répétitives. L'intégrité des télomères est essentielle pour protéger les extrémités des chromosomes contre les systèmes de dégradations et pour les distinguer des cassures de l'ADN double brin. Parce que la machinerie de la réplication de l'ADN n'est pas capable de répliquer l'extrémité des chromosomes, les télomères raccourcissent au fur et à mesure des cycles de réplication. Dès que les télomères atteignent une longueur critique, leur structure protectrice est perdue. Cela induit un signal de dommage de l'ADN et l'arrêt du cycle cellulaire. Pour contrebalancer le raccourcissement des télomères, les cellules qui s'auto régénèrent, dont les cellules de la moelle osseuse, les lymphocytes activés et 80-90% des cellules cancéreuses, expriment la télomérase. C'est une ribonucléoprotéine qui a la capacité de synthétiser des séquences télomériques par transcription inverse d'une courte séquence contenue dans sa propre sous-unité ARN avec laquelle elle est associée. La télomérase humaine est une enzyme processive au niveau de l'addition des nucléotides et aussi des répétitions télomériques. La télomérase de levure et la télomérase humaine sont toutes deux dimériques et il a été montré que la télomérase humaine recombinante contient deux ARN qui coopèrent pour fonctionner ainsi que deux sous-unités catalytiques. Cependant, il n'a pas encore été montré quel est le rôle de la dimérisation dans l'activité de la télomérase. Afin d'élucider ce rôle, nous avons exprimé, reconstitué et purifié la télomérase humaine dimérique recombinante. Et pour étudier l'effet d'ARN mutants sur l'activité de la télomérase, nous avons développé une méthode pour reconstituer et enrichir en hétérodimères de télomérase. Les hétérodimères contiennent une sous-unité ARN sauvage et une sous-unité ARN mutée au niveau de la séquence de la matrice. Sur l'ARN muté nous avons introduit une étiquette aptamer ARN-S1 puis nous avons purifié la télomérase via l'etiquette Si. Nous avons montré que la dimérisation est essentielle pour l'activité de la télomérase. Nos données indiquent que chaque télomérase du dimère allonge leur substrat, l'ADN télomérique, indépendamment l'une de l'autre à chaque cycle d'élongation mais que l'addition itérative de répétitions télomériques nécessite une coopération entre les deux télomérases du dimère. Nous proposons donc un modèle dans lequel les deux télomérases du dimères se lient et allongent deux substrats télomères et que pendant l'élongation processive les deux enzymes subissent un changement de conformation de manière coordonnée, ce changement va permettre le repositionnement des substrats pour d'autres cycles d'additions de répétitions télomériques. Dyskeratosis congenita est une maladie mortelle due majoritairement au disfonctionnement de la moelle osseuse. Dans la forme autosomale de la maladie, l'ARN de la télomérase contient des mutations. En utilisant notre système de reconstitution, nous avons montré que ces ARN mutés, qui ont perdu leur activité enzymatique dans le cas d'un homodimère de mutants, sont dominant négatifs quand ils sont présents dans les hétérodimères sauvage/mutant. Cet effet trans-dominant négatif pourrait contribuer à la progression de la maladie. Abstract Telomeres are protein-DNA structures at the ends of linear eukaryotic chromosomes. The telomeric DNA consists of tandemly repeated sequences. Telomeric integrity is essential to protect chromosomal ends from nucleolytic degradation and to prevent their recognition as DNA double strand breaks. Due to the inability of the conventional DNA replication machinery to replicate terminal DNA stretches, telomeres shorten with continuous rounds of DNA replication. As soon as telomeres reach a critical length, their protective structure is lost and the deprotected telomeres will induce a DNA damage response leading to cell cycle arrest. To counteract telomere shortening, self-renewing cells, including bone marrow cells, activated lymphocytes and 80-90% of cancer cells express the cellular reverse transcriptase telomerase, which has the capacity to synthesize telomeric repeats by reverse transcription of a short template sequence encoded by its stably associated RNA subunit. Human telomerase is a processive enzyme for nucleotide as well as repeat addition. Both yeast and human telomerase are dimeric enzymes and recombinant human telomerase has been shown to contain two functionally cooperating RNAs and most probably also two protein subunits. However, it has remained unclear how dimerization may contribute to telomerase activity. To study the role of dimerization, we expressed, reconstituted and purified recombinant human telomerase. We also developed a new method to reconstitute and enrich for telomerase heterodimers containing wild-type (wt) and mutant telomerase RNA subunits. To this end we introduced an S1-RNA-aptamer tag into telomerase RNA and purified telomerase reconstituted with a mixture of untagged and tagged RNA via the S1-tag. Using this experimental system, we introduced template mutations in the tagged RNA subunit and examined the effect of mutant RNAs on wt telomerase activity in wt/mutant heterodimers. We obtained evidence that dimerization is essential for telomerase activity. Our data indicate that the two subunits elongate telomere substrates independently of each other during single rounds of elongation, but that iterative addition of telomeric repeats requires cooperation between the two subunits. We suggest a model, in which dimeric telomerases bind and elongate two telomere substrates and that the two subunits undergo coordinated conformational changes during processive elongation that enable repositioning the substrates for subsequent rounds of repeat addition. Dyskeratosis congenita is a multisystemic disease with bone marrow failure as the major cause of death. The autosomal form of this disease was found to harbor mutations in the telomerase RNA. Using our reconstitution system, we tested whether mutant dyskeratosis telomerase RNAs behaved in a dominant negative manner. We observed that dyskeratosis telomerase RNA mutants, which lacked enzymatic activity were dominant negative, when present in wt/ mutant heterodimers. The transdominant negative effect of these mutants may contribute to disease progression.

X-chromosome SNP analyses in 11 human Mediterranean populations show a high overall genetic homogeneity except in North-west Africans (Moroccans).

BACKGROUND: Due to its history, with a high number of migration events, the Mediterranean basin represents a challenging area for population genetic studies. A large number of genetic studies have been carried out in the Mediterranean area using different markers but no consensus has been reached on the genetic landscape of the Mediterranean populations. In order to further investigate the genetics of the human Mediterranean populations, we typed 894 individuals from 11 Mediterranean populations with 25 single-nucleotide polymorphisms (SNPs) located on the X-chromosome. RESULTS: A high overall homogeneity was found among the Mediterranean populations except for the population from Morocco, which seemed to differ genetically from the rest of the populations in the Mediterranean area. A very low genetic distance was found between populations in the Middle East and most of the western part of the Mediterranean Sea.A higher migration rate in females versus males was observed by comparing data from X-chromosome, mt-DNA and Y-chromosome SNPs both in the Mediterranean and a wider geographic area.Multilocus association was observed among the 25 SNPs on the X-chromosome in the populations from Ibiza and Cosenza. CONCLUSION: Our results support both the hypothesis of (1) a reduced impact of the Neolithic Wave and more recent migration movements in NW-Africa, and (2) the importance of the Strait of Gibraltar as a geographic barrier. In contrast, the high genetic homogeneity observed in the Mediterranean area could be interpreted as the result of the Neolithic wave caused by a large demic diffusion and/or more recent migration events. A differentiated contribution of males and females to the genetic landscape of the Mediterranean area was observed with a higher migration rate in females than in males. A certain level of background linkage disequilibrium in populations in Ibiza and Cosenza could be attributed to their demographic background.

Mechanisms regulating the proliferative potential of human CD8+ T lymphocytes overexpressing telomerase.

In human somatic cells, including T lymphocytes, telomeres progressively shorten with each cell division, eventually leading to a state of cellular senescence. Ectopic expression of telomerase results in the extension of their replicative life spans without inducing changes associated with transformation. However, it is yet unknown whether somatic cells that overexpress telomerase are physiologically indistinguishable from normal cells. Using CD8+ T lymphocyte clones overexpressing telomerase, we investigated the molecular mechanisms that regulate T cell proliferation. In this study, we show that early passage T cell clones transduced or not with human telomerase reverse transcriptase displayed identical growth rates upon mitogenic stimulation and no marked global changes in gene expression. Surprisingly, reduced proliferative responses were observed in human telomerase reverse transcriptase-transduced cells with extended life spans. These cells, despite maintaining high expression levels of genes involved in the cell cycle progression, also showed increased expression in several genes found in common with normal aging T lymphocytes. Strikingly, late passage T cells overexpressing telomerase accumulated the cyclin-dependent inhibitors p16Ink4a and p21Cip1 that have largely been associated with in vitro growth arrest. We conclude that alternative growth arrest mechanisms such as those mediated by p16Ink4a and p21Cip1 still remained intact and regulated the growth potential of cells independently of their telomere status.

Genometric analyses of the organization of circular chromosomes: a universal pressure determines the direction of ribosomal RNA genes transcription relative to chromosome replication.

Selective pressures related to gene function and chromosomal architecture are acting on genome sequences and can be revealed, for instance, by appropriate genometric methods. Cumulative nucleotide skew analyses, i.e., GC, TA, and ORF orientation skews, predict the location of the origin of DNA replication for 88 out of 100 completely sequenced bacterial chromosomes. These methods appear fully reliable for proteobacteria, Gram-positives, and spirochetes as well as for euryarchaeotes. Based on this genome architecture information, coorientation analyses reveal that in prokaryotes, ribosomal RNA (rRNA) genes encoding the small and large ribosomal subunits are all transcribed in the same direction as DNA replication; that is, they are located along the leading strand. This result offers a simple and reliable method for circumscribing the region containing the origin of the DNA replication and reveals a strong selective pressure acting on the orientation of rRNA genes similar to the weaker one acting on the orientation of ORFs. Rate of coorientation of transfer RNA (tRNA) genes with DNA replication appears to be taxon-specific. Analyzing nucleotide biases such as GC and TA skews of genes and plotting one against the other reveals a taxonomic clusterization of species. All ribosomal RNA genes are enriched in Gs and depleted in Cs, the only so far known exception being the rRNA genes of deuterostomian mitochondria. However, this exception can be explained by the fact that in the chromosome of the human mitochondrion, the model of the deuterostomian organelle genome, DNA replication, and rRNA transcription proceed in opposite directions. A general rule is deduced from prokaryotic and mitochondrial genomes: ribosomal RNA genes that are transcribed in the same direction as the DNA replication are enriched in Gs, and those transcribed in the opposite direction are depleted in Gs.