José Manuel Rodríguez González

The GENCODE v7 catalog of human long noncoding RNAs: Analysis of their gene structure, evolution, and expression

The human genome contains many thousands of long noncoding RNAs (lncRNAs). While several studies have demonstrated compelling biological and disease roles for individual examples, analytical and experimental approaches to investigate these genes have been hampered by the lack of comprehensive lncRNA annotation. Here, we present and analyze the most complete human lncRNA annotation to date, produced by the GENCODE consortium within the framework of the ENCODE project and comprising 9277 manually annotated genes producing 14,880 transcripts. Our analyses indicate that lncRNAs are generated through pathways similar to that of protein-coding genes, with similar histone-modification profiles, splicing signals, and exon/intron lengths. In contrast to protein-coding genes, however, lncRNAs display a striking bias toward two-exon transcripts, they are predominantly localized in the chromatin and nucleus, and a fraction appear to be preferentially processed into small RNAs. They are under stronger selective pressure than neutrally evolving sequences-particularly in their promoter regions, which display levels of selection comparable to protein-coding genes. Importantly, about one-third seem to have arisen within the primate lineage. Comprehensive analysis of their expression in multiple human organs and brain regions shows that lncRNAs are generally lower expressed than protein-coding genes, and display more tissue-specific expression patterns, with a large fraction of tissue-specific lncRNAs expressed in the brain. Expression correlation analysis indicates that lncRNAs show particularly striking positive correlation with the expression of antisense coding genes. This GENCODE annotation represents a valuable resource for future studies of lncRNAs.

GENCODE: The reference human genome annotation for The ENCODE Project

The GENCODE Consortium aims to identify all gene features in the human genome using a combination of computational analysis, manual annotation, and experimental validation. Since the first public release of this annotation data set, few new protein-coding loci have been added, yet the number of alternative splicing transcripts annotated has steadily increased. The GENCODE 7 release contains 20,687 protein-coding and 9640 long noncoding RNA loci and has 33,977 coding transcripts not represented in UCSC genes and RefSeq. It also has the most comprehensive annotation of long noncoding RNA (lncRNA) loci publicly available with the predominant transcript form consisting of two exons. We have examined the completeness of the transcript annotation and found that 35% of transcriptional start sites are supported by CAGE clusters and 62% of protein-coding genes have annotated polyA sites. Over one-third of GENCODE protein-coding genes are supported by peptide hits derived from mass spectrometry spectra submitted to Peptide Atlas. New models derived from the Illumina Body Map 2.0 RNA-seq data identify 3689 new loci not currently in GENCODE, of which 3127 consist of two exon models indicating that they are possibly unannotated long noncoding loci. GENCODE 7 is publicly available from gencodegenes.org and via the Ensembl and UCSC Genome Browsers.

Current status and new features of the Consensus Coding Sequence database

The Consensus Coding Sequence (CCDS) project (http://www.ncbi.nlm.nih.gov/CCDS/) is a collaborative effort to maintain a dataset of protein-coding regions that are identically annotated on the human and mouse reference genome assemblies by the National Center for Biotechnology Information (NCBI) and Ensembl genome annotation pipelines. Identical annotations that pass quality assurance tests are tracked with a stable identifier (CCDS ID). Members of the collaboration, who are from NCBI, the Wellcome Trust Sanger Institute and the University of California Santa Cruz, provide coordinated and continuous review of the dataset to ensure high-quality CCDS representations. We describe here the current status and recent growth in the CCDS dataset, as well as recent changes to the CCDS web and FTP sites. These changes include more explicit reporting about the NCBI and Ensembl annotation releases being compared, new search and display options, the addition of biologically descriptive information and our approach to representing genes for which support evidence is incomplete. We also present a summary of recent and future curation targets.

Distinct Gene Expression Profiling after Infection of Immature Human Monocyte-Derived Dendritic Cells by the Attenuated Poxvirus Vectors MVA and NYVAC

ABSTRACT Although recombinants based on the attenuated poxvirus vectors MVA and NYVAC are currently in clinical trials, the nature of the genes triggered by these vectors in antigen-presenting cells is poorly characterized. Using microarray technology and various analysis conditions, we compared specific changes in gene expression profiling following MVA and NYVAC infection of immature human monocyte-derived dendritic cells (MDDC). Microarray analysis was performed at 6 h postinfection, since these viruses induced extensive cytopathic effects, rRNA breakdown, and apoptosis at late times postinfection. MVA- and NYVAC-infected MDDC shared upregulation of 195 genes compared to uninfected cells: MVA specifically upregulated 359 genes, and NYVAC upregulated 165 genes. Microarray comparison of NYVAC and MVA infection revealed 544 genes with distinct expression patterns after poxvirus infection and 283 genes specifically upregulated after MVA infection. Both vectors upregulated genes for cytokines, cytokine receptors, chemokines, chemokine receptors, and molecules involved in antigen uptake and processing, including major histocompatibility complex genes. mRNA levels for interleukin 12β (IL-12β), beta interferon, and tumor necrosis factor alpha were higher after MVA infection than after NYVAC infection. The expression profiles of transcription factors such as NF-κB/Rel and STAT were regulated similarly by both viruses; in contrast, OASL, MDA5, and IRIG-I expression increased only during MVA infection. Type I interferon, IL-6, and Toll-like receptor pathways were specifically induced after MVA infection. Following MVA or NYVAC infection in MDDC, we found similarities as well as differences between these virus strains in the expression of cellular genes with immunological function, which should have an impact when these vectors are used as recombinant vaccines.

ERK1/2 MAP kinases promote cell cycle entry by rapid, kinase-independent disruption of retinoblastoma–lamin A complexes

As orchestrators of essential cellular processes like proliferation, ERK1/2 mitogen-activated protein kinase signals impact on cell cycle regulation. A-type lamins are major constituents of the nuclear matrix that also control the cell cycle machinery by largely unknown mechanisms. In this paper, we disclose a functional liaison between ERK1/2 and lamin A whereby cell cycle progression is regulated. We demonstrate that lamin A serves as a mutually exclusive dock for ERK1/2 and the retinoblastoma (Rb) protein. Our results reveal that, immediately after their postactivation entrance in the nucleus, ERK1/2 dislodge Rb from its interaction with lamin A, thereby facilitating its rapid phosphorylation and consequently promoting E2F activation and cell cycle entry. Interestingly, these effects are independent of ERK1/2 kinase activity. We also show that cellular transformation and tumor cell proliferation are dependent on the balance between lamin A and nuclear ERK1/2 levels, which determines Rb accessibility for phosphorylation/inactivation.

Stabilization of a Full-Length Infectious cDNA Clone of Transmissible Gastroenteritis Coronavirus by Insertion of an Intron

ABSTRACT The stable propagation of a full-length transmissible gastroenteritis coronavirus (TGEV) cDNA in Escherichia coli cells as a bacterial artificial chromosome has been considerably improved by the insertion of an intron to disrupt a toxic region identified in the viral genome. The viral RNA was expressed in the cell nucleus under the control of the cytomegalovirus promoter and the intron was efficiently removed during translocation of this RNA to the cytoplasm. The insertion in two different positions allowed stable plasmid amplification for at least 200 generations. Infectious TGEV was efficiently recovered from cells transfected with the modified cDNAs.

Induced magnetic anisotropy in Co–Mn–Si–B amorphous microwires

The effect of both thermal treatments as well as chemical etching treatments on the magnetic behavior has been investigated in nearly-zero magnetostriction glass-coated amorphous (Co1−xMnx)75Si10B15 (x=0.08, 0.09, and 0.10) microwires. Such a small change in x results in rather strong changes in the hysteresis loop parameters including coercitivity, Hc, and initial magnetic permeability, μ15. This effect was ascribed to the change of sign of the magnetostriction constant with a proper combination of the coercivity and relatively high initial permeability for the as-cast x=0.09 sample. Thermal treatment (temperatures 100–200 °C for 0.5–2 h) as well as chemical etching in 20% diluted fluoridric acid with duration from 0.5 up to 50 min modify this magnetic parameters owing to the internal stresses relaxation process. In particular, annealing under applied magnetic field (field annealing) can improve more significantly these magnetic parameters: increasing both Hc and μ15. Such phenomenology can be interprete...

Combining RT-PCR-seq and RNA-seq to catalog all genic elements encoded in the human genome

Within the ENCODE Consortium, GENCODE aimed to accurately annotate all protein-coding genes, pseudogenes, and noncoding transcribed loci in the human genome through manual curation and computational methods. Annotated transcript structures were assessed, and less well-supported loci were systematically, experimentally validated. Predicted exon-exon junctions were evaluated by RT-PCR amplification followed by highly multiplexed sequencing readout, a method we called RT-PCR-seq. Seventy-nine percent of all assessed junctions are confirmed by this evaluation procedure, demonstrating the high quality of the GENCODE gene set. RT-PCR-seq was also efficient to screen gene models predicted using the Human Body Map (HBM) RNA-seq data. We validated 73% of these predictions, thus confirming 1168 novel genes, mostly noncoding, which will further complement the GENCODE annotation. Our novel experimental validation pipeline is extremely sensitive, far more than unbiased transcriptome profiling through RNA sequencing, which is becoming the norm. For example, exon-exon junctions unique to GENCODE annotated transcripts are five times more likely to be corroborated with our targeted approach than with extensive large human transcriptome profiling. Data sets such as the HBM and ENCODE RNA-seq data fail sampling of low-expressed transcripts. Our RT-PCR-seq targeted approach also has the advantage of identifying novel exons of known genes, as we discovered unannotated exons in ~11% of assessed introns. We thus estimate that at least 18% of known loci have yet-unannotated exons. Our work demonstrates that the cataloging of all of the genic elements encoded in the human genome will necessitate a coordinated effort between unbiased and targeted approaches, like RNA-seq and RT-PCR-seq.

Comparison of GENCODE and RefSeq gene annotation and the impact of reference geneset on variant effect prediction

BackgroundA vast amount of DNA variation is being identified by increasingly large-scale exome and genome sequencing projects. To be useful, variants require accurate functional annotation and a wide range of tools are available to this end. McCarthy et al recently demonstrated the large differences in prediction of loss-of-function (LoF) variation when RefSeq and Ensembl transcripts are used for annotation, highlighting the importance of the reference transcripts on which variant functional annotation is based.ResultsWe describe a detailed analysis of the similarities and differences between the gene and transcript annotation in the GENCODE and RefSeq genesets. We demonstrate that the GENCODE Comprehensive set is richer in alternative splicing, novel CDSs, novel exons and has higher genomic coverage than RefSeq, while the GENCODE Basic set is very similar to RefSeq. Using RNAseq data we show that exons and introns unique to one geneset are expressed at a similar level to those common to both. We present evidence that the differences in gene annotation lead to large differences in variant annotation where GENCODE and RefSeq are used as reference transcripts, although this is predominantly confined to non-coding transcripts and UTR sequence, with at most ~30% of LoF variants annotated discordantly. We also describe an investigation of dominant transcript expression, showing that it both supports the utility of the GENCODE Basic set in providing a smaller set of more highly expressed transcripts and provides a useful, biologically-relevant filter for further reducing the complexity of the transcriptome.ConclusionsThe reference transcripts selected for variant functional annotation do have a large effect on the outcome. The GENCODE Comprehensive transcripts contain more exons, have greater genomic coverage and capture many more variants than RefSeq in both genome and exome datasets, while the GENCODE Basic set shows a higher degree of concordance with RefSeq and has fewer unique features. We propose that the GENCODE Comprehensive set has great utility for the discovery of new variants with functional potential, while the GENCODE Basic set is more suitable for applications demanding less complex interpretation of functional variants.

Prophylactic thyroidectomy in MEN 2A syndrome: experience in a single center

BACKGROUND Genetic study of the RET proto-oncogene has modified the management, treatment, and prognosis of medullary thyroid carcinoma (MTC), multiple endocrine neoplasia 2A (MEN 2A), for patients with less advanced tumor stages. Classically, the diagnosis was based on an increase in basal and poststimulus peak calcitonin (bCT and pCT). Prophylactic thyroidectomy, based on results of genetic testing, may reduce recurrences in MTC. STUDY DESIGN Of 82 MTC (MEN 2A) patients genetically diagnosed and surgically treated at our center, 22 received a prophylactic thyroidectomy (RET +, bCT and pCT with normal values and asymptomatic). We analyzed age, gender, phenotype, RET mutation, cervical ultrasound, laboratory tests (bCT, pCT, and CEA), surgery, histologic data, TNM, and followup. RESULTS The 22 patients belonged to 8 families with MTC (MEN 2A). Mean age was 15.2 years (range 5 to 36 years). The RET mutation in 21 patients was Cys-->Tyr and in the remaining patient both in codon 634 in exon 11. The median values of bCT and pCT were 38 pg/mL (range < 15 to 75 pg/mL) and 148.5 pg/mL (range < 15 to 250 pg/mL), respectively. Total thyroidectomy was performed in 8 patients (age < or = 10 years) and associated central neck dissection in 14 patients (age> 10 years). Histologic study showed 7 C-cell hyperplasias and 15 MTCs (8 bilateral); the median size was 0.2 cm (range < 0.1 to 0.7cm); 1 patient had metastatic adenopathies. According to TNM, 7 were stage 0, 14 were stage I, and 1 was stage III. Postsurgery bCT and pCT values were normal in all patients, with a curative rate of 100%. MTC patients compared with C-cell hyperplasia patients were older on average, had higher mean bCT, mean pCT, and mean CEA. CONCLUSIONS Prophylactic thyroidectomy based on genetic testing allows identification and treatment of patients at an early stage of the disease and decreases recurrence rates. pCT values above the upper limit of normal may be markers for the presence of MTC and should be considered in selecting operative procedures for these patients.