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.

Clinical follow-up of women infected with human papillomavirus-16, either alone or with other human papillomavirus types: identification of different risk groups.

OBJECTIVES: Evaluation of the clinical impact of multiple infections of the cervix by human papillomavirus, including human papillomavirus-16, compared with single human papillomavirus-16 infection. STUDY DESIGN: One hundred sixty-nine women were classified in 3 categories depending on their human papillomavirus profile: human papillomavirus-16 only, human papillomavirus-16 and low-risk type(s), and human papillomavirus-16 and other high-risk type(s). Cervical brush samples were analyzed for human papillomavirus DNA by polymerase chain reaction and reverse line blot hybridization. All women were evaluated with colposcopy during 24 months or more. Management was according to the Bethesda recommendations. RESULTS: Women infected with human papillomavirus-16 and other high-risk human papillomavirus type(s) presented more progression or no change in the grade of dysplasia, compared with women of the other groups (relative risk [RR], 1.39; 95% confidence interval [CI], 1.07-1.82; P = .02 at 6 months; RR, 2.10; 95% CI, 1.46-3.02; P < .001 at 12 months; RR, 1.82; 95% CI, 1.21-2.72; P = .004 at 24 months). CONCLUSION: Coinfection of women with human papillomavirus-16 and other high-risk human papillomavirus type(s) increases the risk of unfavorable evolution.

IL-13 induces expression of CD36 in human monocytes through PPARgamma activation.

The class B scavenger receptor CD36 is a component of the pattern recognition receptors on monocytes that recognizes a variety of molecules. CD36 expression in monocytes depends on exposure to soluble mediators. We demonstrate here that CD36 expression is induced in human monocytes following exposure to IL-13, a Th2 cytokine, via the peroxisome proliferator-activated receptor (PPAR)gamma pathway. Induction of CD36 protein was paralleled by an increase in CD36 mRNA. The PPARgamma pathway was demonstrated using transfection of a PPARgamma expression plasmid into the murine macrophage cell line RAW264.7, expressing very low levels of PPARgamma, and in peritoneal macrophages from PPARgamma-conditional null mice. We also show that CD36 induction by IL-13 via PPARgamma is dependent on phospholipase A2 activation and that IL-13 induces the production of endogenous 15-deoxy-Delta12,14-prostaglandin J2, an endogenous PPARgamma ligand, and its nuclear localization in human monocytes. Finally, we demonstrate that CD36 and PPARgamma are involved in IL-13-mediated phagocytosis of Plasmodium falciparum-parasitized erythrocytes. These results reveal a novel role for PPARgamma in the alternative activation of monocytes by IL-13, suggesting that endogenous PPARgamma ligands, produced by phospholipase A2 activation, could contribute to the biochemical and cellular functions of CD36.

Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project.

We report the generation and analysis of functional data from multiple, diverse experiments performed on a targeted 1% of the human genome as part of the pilot phase of the ENCODE Project. These data have been further integrated and augmented by a number of evolutionary and computational analyses. Together, our results advance the collective knowledge about human genome function in several major areas. First, our studies provide convincing evidence that the genome is pervasively transcribed, such that the majority of its bases can be found in primary transcripts, including non-protein-coding transcripts, and those that extensively overlap one another. Second, systematic examination of transcriptional regulation has yielded new understanding about transcription start sites, including their relationship to specific regulatory sequences and features of chromatin accessibility and histone modification. Third, a more sophisticated view of chromatin structure has emerged, including its inter-relationship with DNA replication and transcriptional regulation. Finally, integration of these new sources of information, in particular with respect to mammalian evolution based on inter- and intra-species sequence comparisons, has yielded new mechanistic and evolutionary insights concerning the functional landscape of the human genome. Together, these studies are defining a path for pursuit of a more comprehensive characterization of human genome function.

CHO expression of a novel human recombinant IgG1 anti-RhD antibody isolated by phage display.

Replacement of the hyperimmune anti-Rhesus (Rh) D immunoglobulin, currently used to prevent haemolytic disease of the newborn, by fully recombinant human anti-RhD antibodies would solve the current logistic problems associated with supply and demand. The combination of phage display repertoire cloning with precise selection procedures enables isolation of specific genes that can then be inserted into mammalian expression systems allowing production of large quantities of recombinant human proteins. With the aim of selecting high-affinity anti-RhD antibodies, two human Fab libraries were constructed from a hyperimmune donor. Use of a new phage panning procedure involving bromelin-treated red blood cells enabled the isolation of two high-affinity Fab-expressing phage clones. LD-6-3 and LD-6-33, specific for RhD. These showed a novel reaction pattern by recognizing the D variants D(III), D(IVa), D(IVb), D(Va), D(VI) types I and II. D(VII), Rh33 and DFR. Full-length immunoglobulin molecules were constructed by cloning the variable regions into expression vectors containing genomic DNA encoding the immunoglobulin constant regions. We describe the first, stable, suspension growth-adapted Chinese hamster ovary (CHO) cell line producing a high affinity recombinant human IgG1 anti-RhD antibody adapted to pilot-scale production. Evaluation of the Fc region of this recombinant antibody by either chemiluminescence or antibody-dependent cell cytotoxicity (ADCC) assays demonstrated macrophage activation and lysis of red blood cells by human lymphocytes. A consistent source of recombinant human anti-RhD immunoglobulin produced by CHO cells is expected to meet the stringent safety and regulatory requirements for prophylactic application.

AA-amyloidosis caused by visceral leishmaniasis in a human immunodeficiency virus-infected patient.

AA-amyloidosis in the setting of chronic visceral leishmaniasis (VL) has been reported in animal models but documentation in humans is unavailable. Here, we report on a Portuguese man who in 1996 was diagnosed with both human immunodeficiency virus (HIV)-infection and VL. Antiretroviral treatment led to sustained suppression of HIV viremia but CD4+ lymphocytes rose from 8 to only 160 cells/mL. Several courses of antimony treatment did not prevent VL relapses. Renal failure developed in 2006 and renal biopsy revealed AA-amyloidosis. The patient had cryoglobulinemia and serum immune complexes containing antibodies directed against seven leishmanial antigens. Antimony plus amphotericin B, followed by oral miltefosine resulted in a sustained VL treatment response with elimination of circulating Leishmania infantum DNA and CD4+ recovery. The concomitant reduction of serum AA levels and disappearance of circulating leishmanial immune complexes suggests that prolonged VL may lead to AA-amyloidosis in immunocompromised humans.

Widespread occurrence of non-canonical transcription termination by human RNA polymerase III.

Human RNA polymerase (Pol) III-transcribed genes are thought to share a simple termination signal constituted by four or more consecutive thymidine residues in the coding DNA strand, just downstream of the RNA 3'-end sequence. We found that a large set of human tRNA genes (tDNAs) do not display any T(≥4) stretch within 50 bp of 3'-flanking region. In vitro analysis of tDNAs with a distanced T(≥4) revealed the existence of non-canonical terminators resembling degenerate T(≥5) elements, which ensure significant termination but at the same time allow for the production of Pol III read-through pre-tRNAs with unusually long 3' trailers. A panel of such non-canonical signals was found to direct transcription termination of unusual Pol III-synthesized viral pre-miRNA transcripts in gammaherpesvirus 68-infected cells. Genome-wide location analysis revealed that human Pol III tends to trespass into the 3'-flanking regions of tDNAs, as expected from extensive terminator read-through. The widespread occurrence of partial termination suggests that the Pol III primary transcriptome in mammals is unexpectedly enriched in 3'-trailer sequences with the potential to contribute novel functional ncRNAs.

E2F1 mediates DNA damage and apoptosis through HCF-1 and the MLL family of histone methyltransferases.

E2F1 is a key positive regulator of human cell proliferation and its activity is altered in essentially all human cancers. Deregulation of E2F1 leads to oncogenic DNA damage and anti-oncogenic apoptosis. The molecular mechanisms by which E2F1 mediates these two processes are poorly understood but are important for understanding cancer progression. During the G1-to-S phase transition, E2F1 associates through a short DHQY sequence with the cell-cycle regulator HCF-1 together with the mixed-lineage leukaemia (MLL) family of histone H3 lysine 4 (H3K4) methyltransferases. We show here that the DHQY HCF-1-binding sequence permits E2F1 to stimulate both DNA damage and apoptosis, and that HCF-1 and the MLL family of H3K4 methyltransferases have important functions in these processes. Thus, HCF-1 has a broader role in E2F1 function than appreciated earlier. Indeed, sequence changes in the E2F1 HCF-1-binding site can modulate both up and down the ability of E2F1 to induce apoptosis indicating that HCF-1 association with E2F1 is a regulator of E2F1-induced apoptosis.

Persistent human parvovirus B19 infection in children under maintenance chemotherapy for acute lymphocytic leukemia.

OBJECTIVE: To report on B19 infection management and chemotherapy schedule consequences in five children treated for acute lymphocytic leukemia (ALL). PATIENTS AND METHODS: Between May 2001 and February 2002, five patients between 4 and 12 years of age, receiving maintenance chemotherapy for ALL, presented with symptoms suggesting B19 infection (pallor, fatigue, petechiae and pancytopenia in four patients; generalized rash in two patients; acute hepatitis in one patient). Qualitative polymerase chain reaction (PCR) on peripheral blood was used for diagnosis and follow-up of infection; quantitative PCR was used for viral load measurement. Intravenous nonspecific high-dose immunoglobulin therapy was administered until PCR was negative. RESULTS: Qualitative B19 DNA was found in the peripheral blood of all patients, confirming the infection. Viral load at diagnosis ranged from 10 to 10 particles/mL blood. B19 DNA was detectable in four patients at 45, 21, 40, and 44 weeks, respectively. Chemotherapy was delayed in all patients. No clear benefit of intravenous immunoglobulin was noted. CONCLUSIONS: Infection with B19 is rarely reported in patients with ALL, but it should be suspected when unexplained pancytopenia occurs during chemotherapy. Persistent B19 infection remains a challenge in the management of patients receiving maintenance chemotherapy for ALL, as no specific therapy such as a specific immunoglobulin or vaccine exists. The role of viral load measurement needs to be established in terms of its use in follow-up and evaluation of the therapeutic response.

Discovery of amino acid variants in the human glucose-dependent insulinotropic polypeptide (GIP) receptor: the impact on the pancreatic beta cell responses and functional expression studies in Chinese hamster fibroblast cells.

The two incretins, glucose-dependent insulinotropic polypeptide (GIP) and glucagon-like peptide-1 (GLP-1), are insulinotropic factors released from the small intestine to the blood stream in response to oral glucose ingestion. The insulinotropic effect of GLP-1 is maintained in patients with Type II (non-insulin-dependent) diabetes mellitus, whereas, for unknown reasons, the effect of GIP is diminished or lacking. We defined the exon-intron boundaries of the human GIP receptor, made a mutational analysis of the gene and identified two amino acid substitutions, A207 V and E354Q. In an association study of 227 Caucasian Type II diabetic patients and 224 matched glucose tolerant control subjects, the allelic frequency of the A207 V polymorphism was 1.1% in Type II diabetic patients and 0.7% in control subjects (p = 0.48), whereas the allelic frequency of the codon 354 polymorphism was 24.9% in Type II diabetic patients versus 23.2% in control subjects. Interestingly, the glucose tolerant subjects (6% of the population) who were homozygous for the codon 354 variant had on average a 14% decrease in fasting serum C-peptide concentration (p = 0.01) and an 11% decrease in the same variable 30 min after an oral glucose load (p = 0.03) compared with subjects with the wild-type receptor. Investigation of the function of the two GIP receptor variants in Chinese hamster fibroblasts showed, however, that the GIP-induced cAMP formation and the binding of GIP to cells expressing the variant receptors were not different from the findings in cells expressing the wildtype GIP receptor. In conclusion, amino acid variants in the GIP receptor are not associated with random Type II diabetes in patients of Danish Caucasian origin or with altered GIP binding and GIP-induced cAMP production when stably transfected in Chinese hamster fibroblasts. The finding of an association between homozygosity for the codon 354 variant and reduced fasting and post oral glucose tolerance test (OGTT) serum C-peptide concentrations, however, calls for further investigations and could suggest that GIP even in the fasting state regulates the beta-cell secretory response.

Functional and structural basis for a bacteriophage homolog of human RAD52.

In eukaryotes, homologous recombination proteins such as RAD51 and RAD52 play crucial roles in DNA repair and genome stability. Human RAD52 is a member of a large single-strand annealing protein (SSAP) family [1] and stimulates Rad51-dependent recombination [2, 3]. In prokaryotes and phages, it has been difficult to establish the presence of RAD52 homologs with conserved sequences. Putative SSAPs were recently found in several phages that infect strains of Lactococcus lactis[4]. One of these SSAPs was identified as Sak and was found in the virulent L. lactis phage ul36, which belongs to the Siphoviridae family [4, 5]. In this study, we show that Sak is homologous to the N terminus of human RAD52. Purified Sak binds single-stranded DNA (ssDNA) preferentially over double-stranded DNA (dsDNA) and promotes the renaturation of long complementary ssDNAs. Sak also binds RecA and stimulates homologous recombination reactions. Mutations shown to modulate RAD52 DNA binding [6] affect Sak similarly. Remarkably, electron-microscopic reconstruction of Sak reveals an undecameric (11) subunit ring, similar to the crystal structure of the N-terminal fragment of human RAD52 [7, 8]. For the first time, we propose a viral homolog of RAD52 at the amino acid, phylogenic, functional, and structural levels.

Isolation of genetic mutation leading to abnormal phenotype in human through ultra-high-throughput sequencing (UHTS)

The introduction of Next Generation Sequencing (NGS) facilitated the task of localizing DNA variation and identifying the genetic cause of yet unsolved Mendelian disorders. Using Whole Exome Capture method and NGS, we identified the causative genetic aberration responsible for a number of monogenic disorders previously undetermined. Due to the novelty of the NGS method we benchmarked different algorithms to assess their merits and defects. This allowed us to establish a pipeline that we successfully used to pinpoint genes responsible for a form of West's syndrome, a Complex Intellectual Disability syndrome associated with patellar dislocation and celiac disease, and correcting some erroneous molecular diagnosis of Alport's syndrome in a Saudi Arabian family.

The ENCODE (ENCyclopedia Of DNA Elements) Project.

The ENCyclopedia Of DNA Elements (ENCODE) Project aims to identify all functional elements in the human genome sequence. The pilot phase of the Project is focused on a specified 30 megabases (approximately 1%) of the human genome sequence and is organized as an international consortium of computational and laboratory-based scientists working to develop and apply high-throughput approaches for detecting all sequence elements that confer biological function. The results of this pilot phase will guide future efforts to analyze the entire human genome.

Telomerase activity and human telomerase reverse transcriptase mRNA expression in soft tissue tumors: correlation with grade, histology, and proliferative activity.

Telomerase activity (TA) is detected in most human cancers but, with few exceptions, not in normal somatic cells. Little is known about TA in soft tissue tumors. We have examined a series of benign and malignant soft tissue tumors for TA using the telomerase repeat amplification protocol assay. Analysis of the expression of the human telomerase reverse transcriptase was also carried out using RT-PCR. TA was undetectable in benign lesions (15 of 15) and low-grade sarcomas (6 of 6) and was detectable in 50% (19 of 38) of intermediate-/high-grade sarcomas. Although the presence of TA in soft tissue tumors is synonymous with malignancy, it is neither a reliable method in making the distinction between reactive/benign and malignant (especially low-grade) lesions nor a reliable marker of tumor aggressiveness. Leiomyosarcomas and storiform/pleomorphic malignant fibrous histiocytomas rarely showed TA, irrespective of their grade. A strong correlation between human telomerase reverse transcriptase mRNA expression and TA was observed, supporting the close relationship between both parameters. No significant relationship was observed between proliferative activity (as assessed by MIB-1 immunolabeling) and TA. We verified that the absence of telomerase expression was not due to the presence of telomerase inhibitors and therefore alternative mechanism(s) for cell immortalization, yet to be determined, seem to be involved in the development and/or maintenance of some soft tissue sarcomas.

Multiple adenosine 3',5'-cyclic [corrected] monophosphate response element DNA-binding proteins generated by gene diversification and alternative exon splicing.

The protein sequence deduced from the open reading frame of a human placental cDNA encoding a cAMP-responsive enhancer (CRE)-binding protein (CREB-327) has structural features characteristic of several other transcriptional transactivator proteins including jun, fos, C/EBP, myc, and CRE-BP1. Results of Southwestern analysis of nuclear extracts from several different cell lines show that there are multiple CRE-binding proteins, which vary in size in cell lines derived from different tissues and animal species. To examine the molecular diversity of CREB-327 and related proteins at the nucleic acid level, we used labeled cDNAs from human placenta that encode two different CRE-binding proteins (CREB-327 and CRE-BP1) to probe Northern and Southern blots. Both probes hybridized to multiple fragments on Southern blots of genomic DNA from various species. Alternatively, when a human placental c-jun probe was hybridized to the same blot, a single fragment was detected in most cases, consistent with the intronless nature of the human c-jun gene. The CREB-327 probe hybridized to multiple mRNAs, derived from human placenta, ranging in size from 2-9 kilobases. In contrast, the CRE-BP1 probe identified a single 4-kilobase mRNA. Sequence analyses of several overlapping human genomic cosmid clones containing CREB-327 sequences in conjunction with polymerase chain reaction indicates that the CREB-327/341 cDNAs are composed of at least eight or nine exons, and analyses of human placental cDNAs provide direct evidence for at least one alternatively spliced exon. Analyses of mouse/hamster-human hybridoma DNAs by Southern blotting and polymerase chain reaction localizes the CREB-327/341 gene to human chromosome 2. The results indicate that there is a dichotomy of CREB-like proteins, those that are related by overall structure and DNA-binding specificity as well as those that are related by close similarities of primary sequences.