Modeling sequencing errors by combining Hidden Markov models.

Among the largest resources for biological sequence data is the large amount of expressed sequence tags (ESTs) available in public and proprietary databases. ESTs provide information on transcripts but for technical reasons they often contain sequencing errors. Therefore, when analyzing EST sequences computationally, such errors must be taken into account. Earlier attempts to model error prone coding regions have shown good performance in detecting and predicting these while correcting sequencing errors using codon usage frequencies. In the research presented here, we improve the detection of translation start and stop sites by integrating a more complex mRNA model with codon usage bias based error correction into one hidden Markov model (HMM), thus generalizing this error correction approach to more complex HMMs. We show that our method maintains the performance in detecting coding sequences.

An atlas of human gene expression from massively parallel signature sequencing (MPSS).

We have used massively parallel signature sequencing (MPSS) to sample the transcriptomes of 32 normal human tissues to an unprecedented depth, thus documenting the patterns of expression of almost 20,000 genes with high sensitivity and specificity. The data confirm the widely held belief that differences in gene expression between cell and tissue types are largely determined by transcripts derived from a limited number of tissue-specific genes, rather than by combinations of more promiscuously expressed genes. Expression of a little more than half of all known human genes seems to account for both the common requirements and the specific functions of the tissues sampled. A classification of tissues based on patterns of gene expression largely reproduces classifications based on anatomical and biochemical properties. The unbiased sampling of the human transcriptome achieved by MPSS supports the idea that most human genes have been mapped, if not functionally characterized. This data set should prove useful for the identification of tissue-specific genes, for the study of global changes induced by pathological conditions, and for the definition of a minimal set of genes necessary for basic cell maintenance. The data are available on the Web at http://mpss.licr.org and http://sgb.lynxgen.com.

Selection of patients with germline MLH1 mutated Lynch syndrome by determination of MLH1 methylation and BRAF mutation.

Lynch syndrome is one of the most common hereditary colorectal cancer (CRC) syndrome and is caused by germline mutations of MLH1, MSH2 and more rarely MSH6, PMS2, MLH3 genes. Whereas the absence of MSH2 protein is predictive of Lynch syndrome, it is not the case for the absence of MLH1 protein. The purpose of this study was to develop a sensitive and cost effective algorithm to select Lynch syndrome cases among patients with MLH1 immunohistochemical silencing. Eleven sporadic CRC and 16 Lynch syndrome cases with MLH1 protein abnormalities were selected. The BRAF c.1799T> A mutation (p.Val600Glu) was analyzed by direct sequencing after PCR amplification of exon 15. Methylation of MLH1 promoter was determined by Methylation-Sensitive Single-Strand Conformation Analysis. In patients with Lynch syndrome, there was no BRAF mutation and only one case showed MLH1 methylation (6%). In sporadic CRC, all cases were MLH1 methylated (100%) and 8 out of 11 cases carried the above BRAF mutation (73%) whereas only 3 cases were BRAF wild type (27%). We propose the following algorithm: (1) no further molecular analysis should be performed for CRC exhibiting MLH1 methylation and BRAF mutation, and these cases should be considered as sporadic CRC; (2) CRC with unmethylated MLH1 and negative for BRAF mutation should be considered as Lynch syndrome; and (3) only a small fraction of CRC with MLH1 promoter methylation but negative for BRAF mutation should be true Lynch syndrome patients. These potentially Lynch syndrome patients should be offered genetic counselling before searching for MLH1 gene mutations.

Exome sequencing extends the phenotypic spectrum for ABHD12 mutations: from syndromic to nonsyndromic retinal degeneration.

OBJECTIVE: To identify the genetic causes underlying autosomal recessive retinitis pigmentosa (arRP) and to describe the associated phenotype. DESIGN: Case series. PARTICIPANTS: Three hundred forty-seven unrelated families affected by arRP and 33 unrelated families affected by retinitis pigmentosa (RP) plus noncongenital and progressive hearing loss, ataxia, or both, respectively. METHODS: A whole exome sequencing (WES) analysis was performed in 2 families segregating arRP. A mutational screening was performed in 378 additional unrelated families for the exon-intron boundaries of the ABHD12 gene. To establish a genotype-phenotype correlation, individuals who were homozygous or compound heterozygotes of mutations in ABHD12 underwent exhaustive clinical examinations by ophthalmologists, neurologists, and otologists. MAIN OUTCOME MEASURES: DNA sequence variants, best-corrected visual acuity, visual field assessments, electroretinogram responses, magnetic resonance imaging, and audiography. RESULTS: After a WES analysis, we identified 4 new mutations (p.Arg107Glufs*8, p.Trp159*, p.Arg186Pro, and p.Thr202Ile) in ABHD12 in 2 families (RP-1292 and W08-1833) previously diagnosed with nonsyndromic arRP, which cosegregated with the disease among the family members. Another homozygous mutation (p.His372Gln) was detected in 1 affected individual (RP-1487) from a cohort of 378 unrelated arRP and syndromic RP patients. After exhaustive clinical examinations by neurologists and otologists, the 4 affected members of the RP-1292 had no polyneuropathy or ataxia, and the sensorineural hearing loss and cataract were attributed to age or the normal course of the RP, whereas the affected members of the families W08-1833 and RP-1487 showed clearly symptoms associated with polyneuropathy, hearing loss, cerebellar ataxia, RP, and early-onset cataract (PHARC) syndrome. CONCLUSIONS: Null mutations in the ABHD12 gene lead to PHARC syndrome, a neurodegenerative disease including polyneuropathy, hearing loss, cerebellar ataxia, RP, and early-onset cataract. Our study allowed us to report 5 new mutations in ABHD12. This is the first time missense mutations have been described for this gene. Furthermore, these findings are expanding the spectrum of phenotypes associated with ABHD12 mutations ranging from PHARC syndrome to a nonsyndromic form of retinal degeneration.

Real-time PCR for type-specific identification of herpes simplex in clinical samples: evaluation of type-specific results in the context of CNS diseases.

BACKGROUND: HSV-1 and HSV-2 cause CNS infections of dissimilar clinico-pathological characteristics with prognostic and therapeutic implications. OBJECTIVES: To validate a type-specific real-time PCR that uses MGB/LNA Taqman probes and to review the virologico-clinical data of 25 eligible patients with non-neonatal CNS infections. RESULTS: This real-time PCR was evaluated against conventional PCR (26 CSF and 20 quality controls), and LightCycler assay (51 mucocutaneous, 8 CSF and 32 quality controls) and culture/immunofluorescence (75 mucocutaneous) to assess typing with independent methods. Taqman real-time PCR detected 240 HSV genomes per ml CSF, a level appropriate for the management of patients, and provided unambiguous typing for the 104 positive (62 HSV-1 and 42 HSV-2) out the 160 independent clinical samples tested. HSV type diagnosed by Taqman real-time PCR predicted final diagnosis (meningitis versus encephalitis/meningoencephalitis, p<0.001) in 24/25 patients at time of presentation, in contrast to clinical evaluation. CONCLUSIONS: Our real-time PCR, as a sensitive and specific means for type-specific HSV diagnosis, provided rapid prognostic information for patient management.

Use of Treponema pallidum PCR in Testing of Ulcers for Diagnosis of Primary Syphilis(1.).

Treponema pallidum PCR (Tp-PCR) has been noted as a valid method for diagnosing syphilis. We compared Tp-PCR to a combination of darkfield microscopy (DFM), the reference method, and serologic testing in a cohort of 273 patients from France and Switzerland and found the diagnostic accuracy of Tp-PCR was higher than that for DFM.

Sequence conservation in Plasmodium falciparum alpha-helical coiled coil domains proposed for vaccine development.

BACKGROUND: The availability of the P. falciparum genome has led to novel ways to identify potential vaccine candidates. A new approach for antigen discovery based on the bioinformatic selection of heptad repeat motifs corresponding to alpha-helical coiled coil structures yielded promising results. To elucidate the question about the relationship between the coiled coil motifs and their sequence conservation, we have assessed the extent of polymorphism in putative alpha-helical coiled coil domains in culture strains, in natural populations and in the single nucleotide polymorphism data available at PlasmoDB. METHODOLOGY/PRINCIPAL FINDINGS: 14 alpha-helical coiled coil domains were selected based on preclinical experimental evaluation. They were tested by PCR amplification and sequencing of different P. falciparum culture strains and field isolates. We found that only 3 out of 14 alpha-helical coiled coils showed point mutations and/or length polymorphisms. Based on promising immunological results 5 of these peptides were selected for further analysis. Direct sequencing of field samples from Papua New Guinea and Tanzania showed that 3 out of these 5 peptides were completely conserved. An in silico analysis of polymorphism was performed for all 166 putative alpha-helical coiled coil domains originally identified in the P. falciparum genome. We found that 82% (137/166) of these peptides were conserved, and for one peptide only the detected SNPs decreased substantially the probability score for alpha-helical coiled coil formation. More SNPs were found in arrays of almost perfect tandem repeats. In summary, the coiled coil structure prediction was rarely modified by SNPs. The analysis revealed a number of peptides with strictly conserved alpha-helical coiled coil motifs. CONCLUSION/SIGNIFICANCE: We conclude that the selection of alpha-helical coiled coil structural motifs is a valuable approach to identify potential vaccine targets showing a high degree of conservation.

Molecular detection and identification of zygomycetes species from paraffin-embedded tissues in a murine model of disseminated zygomycosis: a collaborative European Society of Clinical Microbiology and Infectious Diseases (ESCMID) Fungal Infection Study Group (EFISG) evaluation.

The present study was performed to assess the interlaboratory reproducibility of the molecular detection and identification of species of Zygomycetes from formalin-fixed paraffin-embedded kidney and brain tissues obtained from experimentally infected mice. Animals were infected with one of five species (Rhizopus oryzae, Rhizopus microsporus, Lichtheimia corymbifera, Rhizomucor pusillus, and Mucor circinelloides). Samples with 1, 10, or 30 slide cuts of the tissues were prepared from each paraffin block, the sample identities were blinded for analysis, and the samples were mailed to each of seven laboratories for the assessment of sensitivity. A protocol describing the extraction method and the PCR amplification procedure was provided. The internal transcribed spacer 1 (ITS1) region was amplified by PCR with the fungal universal primers ITS1 and ITS2 and sequenced. As negative results were obtained for 93% of the tissue specimens infected by M. circinelloides, the data for this species were excluded from the analysis. Positive PCR results were obtained for 93% (52/56), 89% (50/56), and 27% (15/56) of the samples with 30, 10, and 1 slide cuts, respectively. There were minor differences, depending on the organ tissue, fungal species, and laboratory. Correct species identification was possible for 100% (30 cuts), 98% (10 cuts), and 93% (1 cut) of the cases. With the protocol used in the present study, the interlaboratory reproducibility of ITS sequencing for the identification of major Zygomycetes species from formalin-fixed paraffin-embedded tissues can reach 100%, when enough material is available.

Recent advances and prospective researches on molecular epidemiology of dengue viruses

The determination of amino acid changes in the envelop protein by direct sequencing of either genomic RNA or PCR-amplified cDNA fragments provides useful informations for assessing the genetic variability and the geographic distribution of the actually most widespread dengue-2 serotype. The possible link of variations in the envelope protein-gene and virus virulence is discussed.

Global transcriptome sequencing identifies chlamydospore specific markers in Candida albicans and Candida dubliniensis.

Candida albicans and Candida dubliniensis are pathogenic fungi that are highly related but differ in virulence and in some phenotypic traits. During in vitro growth on certain nutrient-poor media, C. albicans and C. dubliniensis are the only yeast species which are able to produce chlamydospores, large thick-walled cells of unknown function. Interestingly, only C. dubliniensis forms pseudohyphae with abundant chlamydospores when grown on Staib medium, while C. albicans grows exclusively as a budding yeast. In order to further our understanding of chlamydospore development and assembly, we compared the global transcriptional profile of both species during growth in liquid Staib medium by RNA sequencing. We also included a C. albicans mutant in our study which lacks the morphogenetic transcriptional repressor Nrg1. This strain, which is characterized by its constitutive pseudohyphal growth, specifically produces masses of chlamydospores in Staib medium, similar to C. dubliniensis. This comparative approach identified a set of putatively chlamydospore-related genes. Two of the homologous C. albicans and C. dubliniensis genes (CSP1 and CSP2) which were most strongly upregulated during chlamydospore development were analysed in more detail. By use of the green fluorescent protein as a reporter, the encoded putative cell wall related proteins were found to exclusively localize to C. albicans and C. dubliniensis chlamydospores. Our findings uncover the first chlamydospore specific markers in Candida species and provide novel insights in the complex morphogenetic development of these important fungal pathogens.

Use of PCR for the determination of the frequency of the deltaF508 mutation in Brasilian cistic fibrosis patients

The [Delta]F508 mutation in the cystic fibrosis (CF) gene was studied in a population of 18 Brazilian CF patients and their 17 families by use of PCR and differential hybridization with oligonucleotides. In a total of 34 chromosomes considered, 12 (35%) carried the F508 deletion, a frequency much lower than that reported in most other populations. As a consequence, CF in Brazil would be predominantly caused by mutations different from the F508 deletion

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.

Rapid bacterial genome sequencing: methods and applications in clinical microbiology.

The recent advances in sequencing technologies have given all microbiology laboratories access to whole genome sequencing. Providing that tools for the automated analysis of sequence data and databases for associated meta-data are developed, whole genome sequencing will become a routine tool for large clinical microbiology laboratories. Indeed, the continuing reduction in sequencing costs and the shortening of the 'time to result' makes it an attractive strategy in both research and diagnostics. Here, we review how high-throughput sequencing is revolutionizing clinical microbiology and the promise that it still holds. We discuss major applications, which include: (i) identification of target DNA sequences and antigens to rapidly develop diagnostic tools; (ii) precise strain identification for epidemiological typing and pathogen monitoring during outbreaks; and (iii) investigation of strain properties, such as the presence of antibiotic resistance or virulence factors. In addition, recent developments in comparative metagenomics and single-cell sequencing offer the prospect of a better understanding of complex microbial communities at the global and individual levels, providing a new perspective for understanding host-pathogen interactions. Being a high-resolution tool, high-throughput sequencing will increasingly influence diagnostics, epidemiology, risk management, and patient care.