Assessment of transcript reconstruction methods for RNA-seq.

We evaluated 25 protocol variants of 14 independent computational methods for exon identification, transcript reconstruction and expression-level quantification from RNA-seq data. Our results show that most algorithms are able to identify discrete transcript components with high success rates but that assembly of complete isoform structures poses a major challenge even when all constituent elements are identified. Expression-level estimates also varied widely across methods, even when based on similar transcript models. Consequently, the complexity of higher eukaryotic genomes imposes severe limitations on transcript recall and splice product discrimination that are likely to remain limiting factors for the analysis of current-generation RNA-seq data.

Interactions of Adeno-associated virus Rep78 protein with the host cell

Abstract : Adeno-associated virus (AAV) is a small DNA virus belonging to the familiy of Parvoviridae. Its genome contains two genes : the rep gene encoding four non structural proteins (Rep78, 68, 52 and 40) implicated in transcription, replication and site-specific integration of the viral DNA and the cap gene encoding three capsid proteins. AAV does not cause any disease, but is studied in view of its potential use to treat several diseases. An interesting property of AAV is its antiproliferative effect. Two elements of AAV can inhibit cell growth. Firstly, the single stranded viral DNA is recognized in cells as damaged DNA leading to either a G2 block or cell death depending on p53 status. Secondly, the two larger Rep proteins (Rep78 and 68) also arrest the cell cycle when they are expressed at high levels. Rep78 in particular induces a complete cell cycle arrest in all the phases, including S phase. Such a strong S phase arrest is rarely seen in other conditions. It was thus interesting to determine how Rep78 could induce it. We found that this strong block is the consequence of Rep78's effects on at least two pathways. Rep78 induces a DNA damage response by producing nicks in the cellular chromatin. Furthermore, Rep78 can bind to the cellular phosphatase Cdc25A and prevent its binding to its substrates CDK2 and CDK1, thus inhibiting its activity. A mutational analysis of Rep78 protein determined that its endonuclease activity is responsible for the DNA damage response and its zinc finger domain for Cdc25A inhibition. The combined expression of two mutants each defective for one of these activities, or these two activities obtained independently of Rep78, could restore the complete cell cycle block, indicating that these two effects of Rep78 are likely to explain completely the cell cycle block it induces. Secondly, the lack of pathogenicity of AAV, its broad range of infection and its ability to integrate site-specifically in human chromosome 19 make it an interesting potential vector for gene therapy. However site-specific integration is only possible in the presence of Rep78/68 whose gene is removed in recombinant AAV vectors. In this part of the study, we tried to introduce Rep protein separately from recombinant AAV vectors to promote their site-specific integration. For that purpose, a fusion protein, TAT-Rep, comprising Rep78/68 joined to the human immunodeficiency virus Tat protein was produced. It had the ability to enter cells and remain active there for a short period. Its activity was sufficient to mediate transcription from the p5 promoter, second-strand synthesis of a recombinant AAV and probably site-specific integration. Résumé : Le virus associé à l'adénovirus (AAV) est un petit virus à ADN qui fait partie de la famille des Parvoviridae. Son génome contient deux gènes : le gène rep code pour quatre protéines (Rep78, 68, 52 et 40) qui participent à la transcription, la réplication et l'intégration du virus et le gène cap code pour les trois protéines de capside. AAV ne produit pas de maladie, mais pourrait au contraire être utilisé pour en soigner. Sa bénignité, sa capacité à infecter différents types de cellules et son intégration spécifique en font un vecteur potentiel pour la thérapie génique. Pour qu'il puisse s'intégrer spécifiquement, il a besoin de la protéine Rep78 ou 68, mais ce gène doit être enlevé des vecteurs pour la thérapie génique. Le but de la première partie de cette étude était d'introduire Rep78 ou 68 dans des cellules en même temps qu'un AAV recombinant, mais indépendamment afin de permettre une intégration spécifique. La stratégie utilisée était de produire une protéine de fusion (TAT-Rep) qui peut entrer dans des cellules si elle est présente dans leur milieu. Cette protéine entrait bien dans les cellules et y était active favorisant ainsi l'intégration spécifique. Une deuxième propriété d'AAV, son effet anti-prolifératif, est intéressante dans le cadre de certaines maladies comme le cancer. Deux éléments d'AAV en sont responsables. D'abord, son ADN simple brin active une réponse cellulaire à l'ADN endommagé et arrête les cellules en G2 ou provoque leur mort. De plus, la protéine Rep78 d'AAV peut fortement bloquer le cycle cellulaire à toutes les phases, même en phase S, ce qui est rare. C'est pourquoi nous avons essayé de comprendre cet effet. Nous avons remarqué que Rep78 doit agir sur deux fronts pour obtenir ce fort bloc. D'un côté, Rep78 introduit des coupures simple brin sur l'ADN de la cellule ce qui active une réponse cellulaire à l'ADN endommagé qui passe par ATM. D'un autre côté, Rep78 lie une phosphatase cellulaire, Cdc25A, et l'empêche ainsi de lier ses substrats CDK2 et CDK1 et donc d'être active. Finalement, à l'aide de mutants de Rep78, nous avons déterminé que l'activité endonuclease de Rep78 était nécessaire pour induire une réponse cellulaire via ATM et que le domaine C-terminal appelé «zinc finger » était responsable de la liaison avec Cdc25A. En co-exprimant deux mutants, qui n'ont chacun qu'un des effets de Rep78, ou en obtenant les deux effets de Rep78 indépendamment d'elle, nous avons obtenu un bloc complet du cycle cellulaire similaire à celui obtenu avec Rep78. Il est donc probable que ces deux effets de Rep78 sont suffisants pour expliquer comment elle arrive à arrêter le cycle cellulaire si efficacement.

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 European dimension for the mouse genome mutagenesis program.

The European Mouse Mutagenesis Consortium is the European initiative contributing to the international effort on functional annotation of the mouse genome. Its objectives are to establish and integrate mutagenesis platforms, gene expression resources, phenotyping units, storage and distribution centers and bioinformatics resources. The combined efforts will accelerate our understanding of gene function and of human health and disease.

Analysis of susceptibility to human immunodeficiency virus 1 (HIV-1) by comparative genetics

Summary [résumé français voir ci-dessous] From the beginning of the 20th century the world population has been confronted with the human immune deficiency virus 1 (HIV-1). This virus has the particularity to mutate fast, and could thus evade and adapt to the human host. Our closest evolutionary related organisms, the non-human primates, are less susceptible to HIV-1. In a broader sense, primates are differentially susceptible to various retrovirus. Species specificity may be due to genetic differences among primates. In the present study we applied evolutionary and comparative genetic techniques to characterize the evolutionary pattern of host cellular determinants of HIV-1 pathogenesis. The study of the evolution of genes coding for proteins participating to the restriction or pathogenesis of HIV-1 may help understanding the genetic basis of modern human susceptibility to infection. To perform comparative genetics analysis, we constituted a collection of primate DNA and RNA to allow generation of de novo sequence of gene orthologs. More recently, release to the public domain of two new primate complete genomes (bornean orang-utan and common marmoset) in addition of the three previously available genomes (human, chimpanzee and Rhesus monkey) help scaling up the evolutionary and comparative genome analysis. Sequence analysis used phylogenetic and statistical methods for detecting molecular adaptation. We identified different selective pressures acting on host proteins involved in HIV-1 pathogenesis. Proteins with HIV-1 restriction properties in non-human primates were under strong positive selection, in particular in regions of interaction with viral proteins. These regions carried key residues for the antiviral activity. Proteins of the innate immunity presented an evolutionary pattern of conservation (purifying selection) but with signals of relaxed constrain if we compared them to the average profile of purifying selection of the primate genomes. Large scale analysis resulted in patterns of evolutionary pressures according to molecular function, biological process and cellular distribution. The data generated by various analyses served to guide the ancestral reconstruction of TRIM5a a potent antiviral host factor. The resurrected TRIM5a from the common ancestor of Old world monkeys was effective against HIV-1 and the recent resurrected hominoid variants were more effective against other retrovirus. Thus, as the result of trade-offs in the ability to restrict different retrovirus, human might have been exposed to HIV-1 at a time when TRIM5a lacked the appropriate specific restriction activity. The application of evolutionary and comparative genetic tools should be considered for the systematical assessment of host proteins relevant in viral pathogenesis, and to guide biological and functional studies. Résumé La population mondiale est confrontée depuis le début du vingtième siècle au virus de l'immunodéficience humaine 1 (VIH-1). Ce virus a un taux de mutation particulièrement élevé, il peut donc s'évader et s'adapter très efficacement à son hôte. Les organismes évolutivement le plus proches de l'homme les primates nonhumains sont moins susceptibles au VIH-1. De façon générale, les primates répondent différemment aux rétrovirus. Cette spécificité entre espèces doit résider dans les différences génétiques entre primates. Dans cette étude nous avons appliqué des techniques d'évolution et de génétique comparative pour caractériser le modèle évolutif des déterminants cellulaires impliqués dans la pathogenèse du VIH- 1. L'étude de l'évolution des gènes, codant pour des protéines impliquées dans la restriction ou la pathogenèse du VIH-1, aidera à la compréhension des bases génétiques ayant récemment rendu l'homme susceptible. Pour les analyses de génétique comparative, nous avons constitué une collection d'ADN et d'ARN de primates dans le but d'obtenir des nouvelles séquences de gènes orthologues. Récemment deux nouveaux génomes complets ont été publiés (l'orang-outan du Bornéo et Marmoset commun) en plus des trois génomes déjà disponibles (humain, chimpanzé, macaque rhésus). Ceci a permis d'améliorer considérablement l'étendue de l'analyse. Pour détecter l'adaptation moléculaire nous avons analysé les séquences à l'aide de méthodes phylogénétiques et statistiques. Nous avons identifié différentes pressions de sélection agissant sur les protéines impliquées dans la pathogenèse du VIH-1. Des protéines avec des propriétés de restriction du VIH-1 dans les primates non-humains présentent un taux particulièrement haut de remplacement d'acides aminés (sélection positive). En particulier dans les régions d'interaction avec les protéines virales. Ces régions incluent des acides aminés clé pour l'activité de restriction. Les protéines appartenant à l'immunité inné présentent un modèle d'évolution de conservation (sélection purifiante) mais avec des traces de "relaxation" comparé au profil général de sélection purifiante du génome des primates. Une analyse à grande échelle a permis de classifier les modèles de pression évolutive selon leur fonction moléculaire, processus biologique et distribution cellulaire. Les données générées par les différentes analyses ont permis la reconstruction ancestrale de TRIM5a, un puissant facteur antiretroviral. Le TRIM5a ressuscité, correspondant à l'ancêtre commun entre les grands singes et les groupe des catarrhiniens, est efficace contre le VIH-1 moderne. Les TRIM5a ressuscités plus récents, correspondant aux ancêtres des grands singes, sont plus efficaces contre d'autres rétrovirus. Ainsi, trouver un compromis dans la capacité de restreindre différents rétrovirus, l'homme aurait été exposé au VIH-1 à une période où TRIM5a manquait d'activité de restriction spécifique contre celui-ci. L'application de techniques d'évolution et de génétique comparative devraient être considérées pour l'évaluation systématique de protéines impliquées dans la pathogenèse virale, ainsi que pour guider des études biologiques et fonctionnelles

Clinical, polysomnographic and genome-wide association analyses of narcolepsy with cataplexy: a European Narcolepsy Network study.

The aim of this study was to describe the clinical and PSG characteristics of narcolepsy with cataplexy and their genetic predisposition by using the retrospective patient database of the European Narcolepsy Network (EU-NN). We have analysed retrospective data of 1099 patients with narcolepsy diagnosed according to International Classification of Sleep Disorders-2. Demographic and clinical characteristics, polysomnography and multiple sleep latency test data, hypocretin-1 levels, and genome-wide genotypes were available. We found a significantly lower age at sleepiness onset (men versus women: 23.74 ± 12.43 versus 21.49 ± 11.83, P = 0.003) and longer diagnostic delay in women (men versus women: 13.82 ± 13.79 versus 15.62 ± 14.94, P = 0.044). The mean diagnostic delay was 14.63 ± 14.31 years, and longer delay was associated with higher body mass index. The best predictors of short diagnostic delay were young age at diagnosis, cataplexy as the first symptom and higher frequency of cataplexy attacks. The mean multiple sleep latency negatively correlated with Epworth Sleepiness Scale (ESS) and with the number of sleep-onset rapid eye movement periods (SOREMPs), but none of the polysomnographic variables was associated with subjective or objective measures of sleepiness. Variant rs2859998 in UBXN2B gene showed a strong association (P = 1.28E-07) with the age at onset of excessive daytime sleepiness, and rs12425451 near the transcription factor TEAD4 (P = 1.97E-07) with the age at onset of cataplexy. Altogether, our results indicate that the diagnostic delay remains extremely long, age and gender substantially affect symptoms, and that a genetic predisposition affects the age at onset of symptoms.

Mutation identification by direct comparison of whole-genome sequencing data from mutant and wild-type individuals using k-mers.

Genes underlying mutant phenotypes can be isolated by combining marker discovery, genetic mapping and resequencing, but a more straightforward strategy for mapping mutations would be the direct comparison of mutant and wild-type genomes. Applying such an approach, however, is hampered by the need for reference sequences and by mutational loads that confound the unambiguous identification of causal mutations. Here we introduce NIKS (needle in the k-stack), a reference-free algorithm based on comparing k-mers in whole-genome sequencing data for precise discovery of homozygous mutations. We applied NIKS to eight mutants induced in nonreference rice cultivars and to two mutants of the nonmodel species Arabis alpina. In both species, comparing pooled F2 individuals selected for mutant phenotypes revealed small sets of mutations including the causal changes. Moreover, comparing M3 seedlings of two allelic mutants unambiguously identified the causal gene. Thus, for any species amenable to mutagenesis, NIKS enables forward genetics without requiring segregating populations, genetic maps and reference sequences.

Landscape of transcription in human cells.

Eukaryotic cells make many types of primary and processed RNAs that are found either in specific subcellular compartments or throughout the cells. A complete catalogue of these RNAs is not yet available and their characteristic subcellular localizations are also poorly understood. Because RNA represents the direct output of the genetic information encoded by genomes and a significant proportion of a cell's regulatory capabilities are focused on its synthesis, processing, transport, modification and translation, the generation of such a catalogue is crucial for understanding genome function. Here we report evidence that three-quarters of the human genome is capable of being transcribed, as well as observations about the range and levels of expression, localization, processing fates, regulatory regions and modifications of almost all currently annotated and thousands of previously unannotated RNAs. These observations, taken together, prompt a redefinition of the concept of a gene.

A genome-wide approach accounting for body mass index identifies genetic variants influencing fasting glycemic traits and insulin resistance.

Recent genome-wide association studies have described many loci implicated in type 2 diabetes (T2D) pathophysiology and β-cell dysfunction but have contributed little to the understanding of the genetic basis of insulin resistance. We hypothesized that genes implicated in insulin resistance pathways might be uncovered by accounting for differences in body mass index (BMI) and potential interactions between BMI and genetic variants. We applied a joint meta-analysis approach to test associations with fasting insulin and glucose on a genome-wide scale. We present six previously unknown loci associated with fasting insulin at P < 5 × 10(-8) in combined discovery and follow-up analyses of 52 studies comprising up to 96,496 non-diabetic individuals. Risk variants were associated with higher triglyceride and lower high-density lipoprotein (HDL) cholesterol levels, suggesting a role for these loci in insulin resistance pathways. The discovery of these loci will aid further characterization of the role of insulin resistance in T2D pathophysiology.

Complete loss of vision caused by a giant mucocele of the frontal sinus.

The authors report the case of a 75-year-old man presenting with an exceptionally large giant posttraumatic mucocele of the frontal sinus years after a gunshot blast to the head. The lesion had grown so extensively that the right eye had shrunk and calcified, resulting in total monocular blindness, a complication that has been reported only once. To the best of our knowledge, it is the first time that a giant mucocele of such a large size is reported. We describe how the patient underwent surgical removal of this massive lesion, cranial base reconstruction, and a cosmetic oculoplastic procedure. The etiology, clinical presentation, and possible complications are reviewed, as well as the importance of a regular clinical follow-up and early surgical cure. Although the diagnosis and management of mucoceles are nowadays considered quite standard, the exceptional size of the lesion illustrated here emphasizes the destructive potential of such seemingly indolent lesions. Despite the benign histology of mucoceles, one should never underestimate their morbid potential or be lulled in delaying surgical cure. Large mucoceles should be removed as quickly as possible to prevent such unacceptable complications as permanent visual loss.

Genomics of secondary metabolite production by Pseudomonas fluorescens Pf-5

The complete sequence of the 7.07 Mb genome of the biological control agent Pseudomonas fluorescens Pf-5 is now available, providing a new opportunity to advance knowledge of biological control through genomics. P. fluorescens Pf-5 is a rhizosphere bacterium that suppresses seedling emergence diseases and produces a spectrum of antibiotics toxic to plant-pathogenic fungi and oomycetes. In addition to six known secondary metabolites produced by Pf-5, three novel secondary metabolite biosynthesis gene clusters identified in the genome could also contribute to biological control. The genomic sequence provides numerous clues as to mechanisms used by the bacterium to survive in the spermosphere and rhizosphere. These features include broad catabolic and transport capabilities for utilizing seed and root exudates, an expanded collection of efflux systems for defense against environmental stress and microbial competition, and the presence of 45 outer membrane receptors that should allow for the uptake of iron from a wide array of siderophores produced by soil microorganisms. As expected for a bacterium with a large genome that lives in a rapidly changing environment, Pf-5 has an extensive collection of regulatory genes, only some of which have been characterized for their roles in regulation of secondary metabolite production or biological control. Consistent with its commensal lifestyle, Pf-5 appears to lack a number of virulence and pathogenicity factors found in plant pathogen.

Genome-wide association analyses identify 18 new loci associated with serum urate concentrations.

Elevated serum urate concentrations can cause gout, a prevalent and painful inflammatory arthritis. By combining data from >140,000 individuals of European ancestry within the Global Urate Genetics Consortium (GUGC), we identified and replicated 28 genome-wide significant loci in association with serum urate concentrations (18 new regions in or near TRIM46, INHBB, SFMBT1, TMEM171, VEGFA, BAZ1B, PRKAG2, STC1, HNF4G, A1CF, ATXN2, UBE2Q2, IGF1R, NFAT5, MAF, HLF, ACVR1B-ACVRL1 and B3GNT4). Associations for many of the loci were of similar magnitude in individuals of non-European ancestry. We further characterized these loci for associations with gout, transcript expression and the fractional excretion of urate. Network analyses implicate the inhibins-activins signaling pathways and glucose metabolism in systemic urate control. New candidate genes for serum urate concentration highlight the importance of metabolic control of urate production and excretion, which may have implications for the treatment and prevention of gout.