The CYSMA web server: An example of integrative tool for in silico analysis of missense variants identified in Mendelian disorders.

Exome sequencing used for molecular diagnosis of Mendelian disorders considerably increases the number of missense variants of unclear significance, whose pathogenicity can be assessed by a variety of prediction tools. As the performance of algorithms may vary according to the datasets, complementary specific resources are needed to improve variant interpretation. As a model, we were interested in the cystic fibrosis transmembrane conductance regulator gene (CFTR) causing cystic fibrosis, in which at least 40% of missense variants are reported. Cystic fibrosis missense analysis (CYSMA) is a new web server designed for online estimation of the pathological relevance of CFTR missense variants. CYSMA generates a set of computationally derived data, ranging from evolutionary conservation to functional observations from three-dimensional structures, provides all available allelic frequencies, clinical observations, and references for functional studies. Compared to software classically used in analysis pipelines on a dataset of 141 well-characterized missense variants, CYSMA was the most efficient tool to discriminate benign missense variants, with a specificity of 85%, and very good sensitivity of 89%. These results suggest that such integrative tools could be adapted to numbers of genes involved in Mendelian disorders to improve the interpretation of missense variants identified in the context of diagnosis.

NF-Y controls fidelity of transcription initiation at gene promoters through maintenance of the nucleosome-depleted region.

Faithful transcription initiation is critical for accurate gene expression, yet the mechanisms underlying specific transcription start site (TSS) selection in mammals remain unclear. Here, we show that the histone-fold domain protein NF-Y, a ubiquitously expressed transcription factor, controls the fidelity of transcription initiation at gene promoters in mouse embryonic stem cells. We report that NF-Y maintains the region upstream of TSSs in a nucleosome-depleted state while simultaneously protecting this accessible region against aberrant and/or ectopic transcription initiation. We find that loss of NF-Y binding in mammalian cells disrupts the promoter chromatin landscape, leading to nucleosomal encroachment over the canonical TSS. Importantly, this chromatin rearrangement is accompanied by upstream relocation of the transcription pre-initiation complex and ectopic transcription initiation. Further, this phenomenon generates aberrant extended transcripts that undergo translation, disrupting gene expression profiles. These results suggest NF-Y is a central player in TSS selection in metazoans and highlight the deleterious consequences of inaccurate transcription initiation.

Ribo-seq enlightens codon usage bias.

Codon usage is biased between lowly and highly expressed genes in a genome-specific manner. This universal bias has been well assessed in some unicellular species, but remains problematic to assess in more complex species. We propose a new method to compute codon usage bias based on genome wide translational data. A new technique based on sequencing of ribosome protected mRNA fragments (Ribo-seq) allowed us to rank genes and compute codon usage bias with high precision for a great variety of species, including mammals. Genes ranking using Ribo-Seq data confirms the influence of the tRNA pool on codon usage bias and shows a decreasing bias in multicellular species. Ribo-Seq analysis also makes possible to detect preferred codons without information on genes function.

Intron retention-dependent gene regulation in Cryptococcus neoformans.

The biological impact of alternative splicing is poorly understood in fungi, although recent studies have shown that these microorganisms are usually intron-rich. In this study, we re-annotated the genome of C. neoformans var. neoformans using RNA-Seq data. Comparison with C. neoformans var. grubii revealed that more than 99% of ORF-introns are in the same exact position in the two varieties whereas UTR-introns are much less evolutionary conserved. We also confirmed that alternative splicing is very common in C. neoformans, affecting nearly all expressed genes. We also observed specific regulation of alternative splicing by environmental cues in this yeast. However, alternative splicing does not appear to be an efficient method to diversify the C. neoformans proteome. Instead, our data suggest the existence of an intron retention-dependent mechanism of gene expression regulation that is not dependent on NMD. This regulatory process represents an additional layer of gene expression regulation in fungi and provides a mechanism to tune gene expression levels in response to any environmental modification.

Analysis of the genome and transcriptome of Cryptococcus neoformans var. grubii reveals complex RNA expression and microevolution leading to virulence attenuation.

Cryptococcus neoformans is a pathogenic basidiomycetous yeast responsible for more than 600,000 deaths each year. It occurs as two serotypes (A and D) representing two varieties (i.e. grubii and neoformans, respectively). Here, we sequenced the genome and performed an RNA-Seq-based analysis of the C. neoformans var. grubii transcriptome structure. We determined the chromosomal locations, analyzed the sequence/structural features of the centromeres, and identified origins of replication. The genome was annotated based on automated and manual curation. More than 40,000 introns populating more than 99% of the expressed genes were identified. Although most of these introns are located in the coding DNA sequences (CDS), over 2,000 introns in the untranslated regions (UTRs) were also identified. Poly(A)-containing reads were employed to locate the polyadenylation sites of more than 80% of the genes. Examination of the sequences around these sites revealed a new poly(A)-site-associated motif (AUGHAH). In addition, 1,197 miscRNAs were identified. These miscRNAs can be spliced and/or polyadenylated, but do not appear to have obvious coding capacities. Finally, this genome sequence enabled a comparative analysis of strain H99 variants obtained after laboratory passage. The spectrum of mutations identified provides insights into the genetics underlying the micro-evolution of a laboratory strain, and identifies mutations involved in stress responses, mating efficiency, and virulence.

Hydrophobic pulses predict transmembrane helix irregularities and channel transmembrane units.

BACKGROUND: Few high-resolution structures of integral membranes proteins are available, as crystallization of such proteins needs yet to overcome too many technical limitations. Nevertheless, prediction of their transmembrane (TM) structure by bioinformatics tools provides interesting insights on the topology of these proteins. METHODS: We describe here how to extract new information from the analysis of hydrophobicity variations or hydrophobic pulses (HPulses) in the sequence of integral membrane proteins using the Hydrophobic Pulse Predictor, a new tool we developed for this purpose. To analyze the primary sequence of 70 integral membrane proteins we defined two levels of analysis: G1-HPulses for sliding windows of n = 2 to 6 and G2-HPulses for sliding windows of n = 12 to 16. RESULTS: The G2-HPulse analysis of 541 transmembrane helices allowed the definition of the new concept of transmembrane unit (TMU) that groups together transmembrane helices and segments with potential adjacent structures. In addition, the G1-HPulse analysis identified helix irregularities that corresponded to kinks, partial helices or unannotated structural events. These irregularities could represent key dynamic elements that are alternatively activated depending on the channel status as illustrated by the crystal structures of the lactose permease in different conformations. CONCLUSIONS: Our results open a new way in the understanding of transmembrane secondary structures: hydrophobicity through hydrophobic pulses strongly impacts on such embedded structures and is not confined to define the transmembrane status of amino acids.

p.Ser1235Arg should no longer be considered as a cystic fibrosis mutation: results from a large collaborative study.

Among the 1700 mutations reported in the cystic fibrosis transmembrane conductance regulator (CFTR) gene, a missense mutation, p.Ser1235Arg, is a relatively frequent finding. To clarify its clinical significance, we collected data from 104 subjects heterozygous for the mutation p.Ser1235Arg from the French CF network, addressed for various indications including classical CF, atypical phenotypes or carrier screening in subjects with or without a family history. Among them, 26 patients (5 having CF, 10 CBAVD (congenital bilateral absence of the vas deferens) and 11 with CF-like symptoms) and 14 healthy subjects were compound heterozygous for a second CFTR mutation. An exhaustive CFTR gene analysis identified a second mutation in cis of p.Ser1235Arg in all CF patients and in 81.8% CBAVD patients. Moreover, epidemiological data from >2100 individuals found a higher frequency of p.Ser1235Arg in the general population than in CF or CBAVD patients. These data, added to the fact that in silico analysis and functional assays suggest a benign nature of this substitution, give several lines of evidence against an association of p.Ser1235Arg with CF or CBAVD.

UMD-CFTR: a database dedicated to CF and CFTR-related disorders.

With the increasing knowledge of cystic fibrosis (CF) and CFTR-related diseases (CFTR-RD), the number of sequence variations in the CFTR gene is constantly raising. CF and particularly CFTR-RD provide a particular challenge because of many unclassified variants and identical genotypes associated with different phenotypes. Using the Universal Mutation Database (UMD) software we have constructed UMD-CFTR (freely available at the URL: http://www.umd.be/CFTR/), the first comprehensive relational CFTR database that allows an in-depth analysis and annotation of all variations identified in individuals whose CFTR genes have been analyzed extensively. The system has been tested on the molecular data from 757 patients (540 CF and 217 CBAVD) including disease-causing, unclassified, and nonpathogenic alterations (301 different sequence variations) representing 3,973 entries. Tools are provided to assess the pathogenicity of mutations. UMD-CFTR also offers a number of query tools and graphical views providing instant access to the list of mutations, their frequencies, positions and predicted consequences, or correlations between genotypes, haplotypes, and phenotypes. UMD-CFTR offers a way to compile not only disease-causing genotypes but also haplotypes. It will help the CFTR scientific and medical communities to improve sequence variation interpretation, evaluate the putative influence of haplotypes on mutations, and correlate molecular data with phenotypes.