Mutations at a single codon in Mad homology 2 domain of SMAD4 cause Myhre syndrome.

Myhre syndrome (MIM 139210) is a developmental disorder characterized by short stature, short hands and feet, facial dysmorphism, muscular hypertrophy, deafness and cognitive delay. Using exome sequencing of individuals with Myhre syndrome, we identified SMAD4 as a candidate gene that contributes to this syndrome on the basis of its pivotal role in the bone morphogenetic pathway (BMP) and transforming growth factor (TGF)-β signaling. We identified three distinct heterozygous missense SMAD4 mutations affecting the codon for Ile500 in 11 individuals with Myhre syndrome. All three mutations are located in the region of SMAD4 encoding the Mad homology 2 (MH2) domain near the site of monoubiquitination at Lys519, and we found a defect in SMAD4 ubiquitination in fibroblasts from affected individuals. We also observed decreased expression of downstream TGF-β target genes, supporting the idea of impaired TGF-β-mediated transcriptional control in individuals with Myhre syndrome.

Improved efficiency of a Salmonella-based vaccine against human papillomavirus type 16 virus-like particles achieved by using a codon-optimized version of L1.

Cervical cancer results from cervical infection by human papillomaviruses (HPVs), especially HPV16. An effective vaccine against these HPVs is expected to have a dramatic impact on the incidence of this cancer and its precursor lesions. The leading candidate, a subunit prophylactic HPV virus-like particle (VLP) vaccine, can protect women from HPV infection. An alternative improved vaccine that avoids parenteral injection, that is efficient with a single dose, and that induces mucosal immunity might greatly facilitate vaccine implementation in different settings. In this study, we have constructed a new generation of recombinant Salmonella organisms that assemble HPV16 VLPs and induce high titers of neutralizing antibodies in mice after a single nasal or oral immunization with live bacteria. This was achieved through the expression of a HPV16 L1 capsid gene whose codon usage was optimized to fit with the most frequently used codons in Salmonella. Interestingly, the high immunogenicity of the new recombinant bacteria did not correlate with an increased expression of L1 VLPs but with a greater stability of the L1-expressing plasmid in vitro and in vivo in absence of antibiotic selection. Anti-HPV16 humoral and neutralizing responses were also observed with different Salmonella enterica serovar Typhimurium strains whose attenuating deletions have already been shown to be safe after oral vaccination of humans. Thus, our findings are a promising improvement toward a vaccine strain that could be tested in human volunteers.

A generalized mechanistic codon model.

Models of codon evolution have attracted particular interest because of their unique capabilities to detect selection forces and their high fit when applied to sequence evolution. We described here a novel approach for modeling codon evolution, which is based on Kronecker product of matrices. The 61 × 61 codon substitution rate matrix is created using Kronecker product of three 4 × 4 nucleotide substitution matrices, the equilibrium frequency of codons, and the selection rate parameter. The entities of the nucleotide substitution matrices and selection rate are considered as parameters of the model, which are optimized by maximum likelihood. Our fully mechanistic model allows the instantaneous substitution matrix between codons to be fully estimated with only 19 parameters instead of 3,721, by using the biological interdependence existing between positions within codons. We illustrate the properties of our models using computer simulations and assessed its relevance by comparing the AICc measures of our model and other models of codon evolution on simulations and a large range of empirical data sets. We show that our model fits most biological data better compared with the current codon models. Furthermore, the parameters in our model can be interpreted in a similar way as the exchangeability rates found in empirical codon models.

KCM, a general framework for codon-based model of evolution

With the advancement of high-throughput sequencing and dramatic increase of available genetic data, statistical modeling has become an essential part in the field of molecular evolution. Statistical modeling results in many interesting discoveries in the field, from detection of highly conserved or diverse regions in a genome to phylogenetic inference of species evolutionary history Among different types of genome sequences, protein coding regions are particularly interesting due to their impact on proteins. The building blocks of proteins, i.e. amino acids, are coded by triples of nucleotides, known as codons. Accordingly, studying the evolution of codons leads to fundamental understanding of how proteins function and evolve. The current codon models can be classified into three principal groups: mechanistic codon models, empirical codon models and hybrid ones. The mechanistic models grasp particular attention due to clarity of their underlying biological assumptions and parameters. However, they suffer from simplified assumptions that are required to overcome the burden of computational complexity. The main assumptions applied to the current mechanistic codon models are (a) double and triple substitutions of nucleotides within codons are negligible, (b) there is no mutation variation among nucleotides of a single codon and (c) assuming HKY nucleotide model is sufficient to capture essence of transition- transversion rates at nucleotide level. In this thesis, I develop a framework of mechanistic codon models, named KCM-based model family framework, based on holding or relaxing the mentioned assumptions. Accordingly, eight different models are proposed from eight combinations of holding or relaxing the assumptions from the simplest one that holds all the assumptions to the most general one that relaxes all of them. The models derived from the proposed framework allow me to investigate the biological plausibility of the three simplified assumptions on real data sets as well as finding the best model that is aligned with the underlying characteristics of the data sets. -- Avec l'avancement de séquençage à haut débit et l'augmentation dramatique des données géné¬tiques disponibles, la modélisation statistique est devenue un élément essentiel dans le domaine dé l'évolution moléculaire. Les résultats de la modélisation statistique dans de nombreuses découvertes intéressantes dans le domaine de la détection, de régions hautement conservées ou diverses dans un génome de l'inférence phylogénétique des espèces histoire évolutive. Parmi les différents types de séquences du génome, les régions codantes de protéines sont particulièrement intéressants en raison de leur impact sur les protéines. Les blocs de construction des protéines, à savoir les acides aminés, sont codés par des triplets de nucléotides, appelés codons. Par conséquent, l'étude de l'évolution des codons mène à la compréhension fondamentale de la façon dont les protéines fonctionnent et évoluent. Les modèles de codons actuels peuvent être classés en trois groupes principaux : les modèles de codons mécanistes, les modèles de codons empiriques et les hybrides. Les modèles mécanistes saisir une attention particulière en raison de la clarté de leurs hypothèses et les paramètres biologiques sous-jacents. Cependant, ils souffrent d'hypothèses simplificatrices qui permettent de surmonter le fardeau de la complexité des calculs. Les principales hypothèses retenues pour les modèles actuels de codons mécanistes sont : a) substitutions doubles et triples de nucleotides dans les codons sont négligeables, b) il n'y a pas de variation de la mutation chez les nucléotides d'un codon unique, et c) en supposant modèle nucléotidique HKY est suffisant pour capturer l'essence de taux de transition transversion au niveau nucléotidique. Dans cette thèse, je poursuis deux objectifs principaux. Le premier objectif est de développer un cadre de modèles de codons mécanistes, nommé cadre KCM-based model family, sur la base de la détention ou de l'assouplissement des hypothèses mentionnées. En conséquence, huit modèles différents sont proposés à partir de huit combinaisons de la détention ou l'assouplissement des hypothèses de la plus simple qui détient toutes les hypothèses à la plus générale qui détend tous. Les modèles dérivés du cadre proposé nous permettent d'enquêter sur la plausibilité biologique des trois hypothèses simplificatrices sur des données réelles ainsi que de trouver le meilleur modèle qui est aligné avec les caractéristiques sous-jacentes des jeux de données. Nos expériences montrent que, dans aucun des jeux de données réelles, tenant les trois hypothèses mentionnées est réaliste. Cela signifie en utilisant des modèles simples qui détiennent ces hypothèses peuvent être trompeuses et les résultats de l'estimation inexacte des paramètres. Le deuxième objectif est de développer un modèle mécaniste de codon généralisée qui détend les trois hypothèses simplificatrices, tandis que d'informatique efficace, en utilisant une opération de matrice appelée produit de Kronecker. Nos expériences montrent que sur un jeux de données choisis au hasard, le modèle proposé de codon mécaniste généralisée surpasse autre modèle de codon par rapport à AICc métrique dans environ la moitié des ensembles de données. En outre, je montre à travers plusieurs expériences que le modèle général proposé est biologiquement plausible.

Functional studies of a germ-line polymorphism at codon 47 within the p53 gene.

A rare germ-line polymorphism in codon 47 of the p53 gene replaces the wild-type proline (CCG) with a serine (TCG). Restriction analysis of 101 human samples revealed the frequency of the rare allele to be 0% (n = 69) in Caucasians and 4.7% (3/64, n = 32) among African-Americans. To investigate the consequence of this amino acid substitution, a cDNA construct (p53 mut47ser) containing the mutation was introduced into a lung adenocarcinoma cell line (Calu-6) that does not express p53. A growth suppression similar to that obtained after introduction of a wild-type p53 cDNA construct was observed, in contrast to the result obtained by introduction of p53 mut143ala. Furthermore, expression of neither p53 mut47ser nor wild-type p53 was tolerated by growing cells. In transient expression assays, both mut47ser and wild-type p53 activated the expression of a reporter gene linked to a p53 binding sequence (PG13-CAT) and inhibited the expression of the luciferase gene under the control of the Rous sarcoma virus promoter (RSVluc). In the same assay, mut143ala did not activate the expression of PG13-CAT and produced only a slight inhibitory effect on RSVluc. These findings indicate that the p53 variant with a serine at codon 47 should be considered as a rare germ-line polymorphism that does not alter the growth-suppression activity of p53.

Transition mutation in codon 248 of the p53 tumor suppressor gene induced by reactive oxygen species and a nitric oxide-releasing compound.

Exposing the human bronchial epithelial cell line BEAS-2B to the nitric oxide (NO) donor sodium 1-(N,N-diethylamino)diazen-1-ium-1, 2-diolate (DEA/NO) at an initial concentration of 0.6 mM while generating superoxide ion at the rate of 1 microM/min with the hypoxanthine/xanthine oxidase (HX/XO) system induced C:G-->T:A transition mutations in codon 248 of the p53 gene. This pattern of mutagenicity was not seen by 'fish-restriction fragment length polymorphism/polymerase chain reaction' (fish-RFLP/PCR) on exposure to DEA/NO alone, however, exposure to HX/XO led to various mutations, suggesting that co-generation of NO and superoxide was responsible for inducing the observed point mutation. DEA/NO potentiated the ability of HX/XO to induce lipid peroxidation as well as DNA single- and double-strand breaks under these conditions, while 0.6 mM DEA/NO in the absence of HX/XO had no significant effect on these parameters. The results show that a point mutation seen at high frequency in certain common human tumors can be induced by simultaneous exposure to reactive oxygen species and a NO source.

Flippase (FLP) recombinase-mediated marker recycling in the dermatophyte Arthroderma vanbreuseghemii.

Biological processes can be elucidated by investigating complex networks of relevant factors and genes. However, this is not possible in species for which dominant selectable markers for genetic studies are unavailable. To overcome the limitation in selectable markers for the dermatophyte Arthroderma vanbreuseghemii (anamorph: Trichophyton mentagrophytes), we adapted the flippase (FLP) recombinase-recombination target (FRT) site-specific recombination system from the yeast Saccharomyces cerevisiae as a selectable marker recycling system for this fungus. Taking into account practical applicability, we designed FLP/FRT modules carrying two FRT sequences as well as the flp gene adapted to the pathogenic yeast Candida albicans (caflp) or a synthetic codon-optimized flp (avflp) gene with neomycin resistance (nptII) cassette for one-step marker excision. Both flp genes were under control of the Trichophyton rubrum copper-repressible promoter (PCTR4). Molecular analyses of resultant transformants showed that only the avflp-harbouring module was functional in A. vanbreuseghemii. Applying this system, we successfully produced the Ku80 recessive mutant strain devoid of any selectable markers. This strain was subsequently used as the recipient for sequential multiple disruptions of secreted metalloprotease (fungalysin) (MEP) or serine protease (SUB) genes, producing mutant strains with double MEP or triple SUB gene deletions. These results confirmed the feasibility of this system for broad-scale genetic manipulation of dermatophytes, advancing our understanding of functions and networks of individual genes in these fungi.

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.

Evolutionary switch and genetic convergence on rbcL following the evolution of C4 photosynthesis.

Rubisco is responsible for the fixation of CO2 into organic compounds through photosynthesis and thus has a great agronomic importance. It is well established that this enzyme suffers from a slow catalysis, and its low specificity results into photorespiration, which is considered as an energy waste for the plant. However, natural variations exist, and some Rubisco lineages, such as in C4 plants, exhibit higher catalytic efficiencies coupled to lower specificities. These C4 kinetics could have evolved as an adaptation to the higher CO2 concentration present in C4 photosynthetic cells. In this study, using phylogenetic analyses on a large data set of C3 and C4 monocots, we showed that the rbcL gene, which encodes the large subunit of Rubisco, evolved under positive selection in independent C4 lineages. This confirms that selective pressures on Rubisco have been switched in C4 plants by the high CO2 environment prevailing in their photosynthetic cells. Eight rbcL codons evolving under positive selection in C4 clades were involved in parallel changes among the 23 independent monocot C4 lineages included in this study. These amino acids are potentially responsible for the C4 kinetics, and their identification opens new roads for human-directed Rubisco engineering. The introgression of C4-like high-efficiency Rubisco would strongly enhance C3 crop yields in the future CO2-enriched atmosphere.

Clinical and genetic investigation of a large Tunisian family with complete achromatopsia: identification of a new nonsense mutation in GNAT2 gene.

Complete achromatopsia is a rare autosomal recessive disease associated with CNGA3, CNGB3, GNAT2 and PDE6C mutations. This retinal disorder is characterized by complete loss of color discrimination due to the absence or alteration of the cones function. The purpose of the present study was the clinical and the genetic characterization of achromatopsia in a large consanguineous Tunisian family. Ophthalmic evaluation included a full clinical examination, color vision testing and electroretinography. Linkage analysis using microsatellite markers flanking CNGA3, CNGB3, GNAT2 and PDE6C genes was performed. Mutations were screened by direct sequencing. A total of 12 individuals were diagnosed with congenital complete achromatopsia. They are members of six nuclear consanguineous families belonging to the same large consanguineous family. Linkage analysis revealed linkage to GNAT2. Mutational screening of GNAT2 revealed three intronic variations c.119-69G>C, c.161+66A>T and c.875-31G>C that co-segregated with a novel mutation p.R313X. An identical GNAT2 haplotype segregating with this mutation was identified, indicating a founder mutation. All patients were homozygous for the p.R313X mutation. This is the first report of the clinical and genetic investigation of complete achromatopsia in North Africa and the largest family with recessive achromatopsia involving GNAT2; thus, providing a unique opportunity for genotype-phenotype correlation for this extremely rare condition.

Pervasive positive selection on duplicated and nonduplicated vertebrate protein coding genes.

A stringent branch-site codon model was used to detect positive selection in vertebrate evolution. We show that the test is robust to the large evolutionary distances involved. Positive selection was detected in 77% of 884 genes studied. Most positive selection concerns a few sites on a single branch of the phylogenetic tree: Between 0.9% and 4.7% of sites are affected by positive selection depending on the branches. No functional category was overrepresented among genes under positive selection. Surprisingly, whole genome duplication had no effect on the prevalence of positive selection, whether the fish-specific genome duplication or the two rounds at the origin of vertebrates. Thus positive selection has not been limited to a few gene classes, or to specific evolutionary events such as duplication, but has been pervasive during vertebrate evolution.

Complete hypogonadotropic hypogonadism associated with a novel inactivating mutation of the gonadotropin-releasing hormone receptor.

In this study, we describe a patient with a phenotype of complete hypogonadotropic hypogonadism who presented primary failure of pulsatile GnRH therapy, but responded to exogenous gonadotropin administration. This patient bore a novel point mutation (T for A) at codon 168 of the gene encoding the GnRH receptor (GnRH-R), resulting in a serine to arginine change in the fourth transmembrane domain of the receptor. This novel mutation was present in the homozygous state in the patient, whereas it was in the heterozygous state in both phenotypically normal parents. When introduced into the complementary DNA coding for the GnRH-R, this mutation resulted in the complete loss of the receptor-mediated signaling response to GnRH. In conclusion, we report the first mutation of the GnRH-R gene that can induce a total loss of function of this receptor and is associated with a phenotype of complete hypogonadotropic hypogonadism.

Chondrodysplasia and Abnormal Joint Development Associated with Mutations in IMPAD1, Encoding the Golgi-Resident Nucleotide Phosphatase, gPAPP.

We used whole-exome sequencing to study three individuals with a distinct condition characterized by short stature, chondrodysplasia with brachydactyly, congenital joint dislocations, cleft palate, and facial dysmorphism. Affected individuals carried homozygous missense mutations in IMPAD1, the gene coding for gPAPP, a Golgi-resident nucleotide phosphatase that hydrolyzes phosphoadenosine phosphate (PAP), the byproduct of sulfotransferase reactions, to AMP. The mutations affected residues in or adjacent to the phosphatase active site and are predicted to impair enzyme activity. A fourth unrelated patient was subsequently found to be homozygous for a premature termination codon in IMPAD1. Impad1 inactivation in mice has previously been shown to produce chondrodysplasia with abnormal joint formation and impaired proteoglycan sulfation. The human chondrodysplasia associated with gPAPP deficiency joins a growing number of skeletoarticular conditions associated with defective synthesis of sulfated proteoglycans, highlighting the importance of proteoglycans in the development of skeletal elements and joints.

Novel testis germ cell-specific transcript of the CREB gene contains an alternatively spliced exon with multiple in-frame stop codons.

The gene encoding the cAMP-responsive transcription factor CREB consists of multiple small exons some of which undergo alternative RNA splicing. We describe the finding of a novel transcript of the CREB gene expressed at high levels in the germ cells of the rat testis. The transcript contains an alternatively spliced exon inserted within the sequence encoding the transcriptional transactivation domain of CREB and this exon contains multiple in-frame stop codons. Furthermore, the exon is conserved in both rat and human genes (75% nucleotide identity). Although the function(s) of this RNA or the truncated CREB protein predicted to result from the translation of this unusual transcript is unknown, the high level of expression in the testicular germ cells and remarkable conservation of sequences in rat and human suggests that it may have a unique biological function in these cells.

Periodic accumulation of cdc15 mRNA is not necessary for septation in Schizosaccharomyces pombe.

Analysis of Schizosaccharomyces pombe mutants that are defective in septum formation and cytokinesis has identified the product of the cdc15 gene as a key element in formation of a division septum. S. pombe cells lacking cdc15p function cannot assemble a functional medial ring, and do not make a division septum. cdc15 mRNA accumulates periodically during the cell cycle, peaking after entry into mitosis, and increased expression of the gene in G2-arrested cells can promote F-actin ring formation. Here, we have investigated the effects of mutations that block cell division upon the expression of cdc15 in synchronised cell populations, and analysed the expression of cdc15 when septum formation is induced by ectopic activation of the septation signalling network. We concluded the following: (i) the septation signalling network genes are not required for periodic accumulation of cdc15 mRNA; (ii) induction of septum formation in G2-arrested cells by activation of the septation signalling network does not result in accumulation of cdc15 mRNA, which is therefore not a prerequisite for septum formation; (iii) failure to turn off septum formation at the end of mitosis results in continued expression of cdc15; and (iv) periodic accumulation of cdc15 mRNA is mediated by a 97 bp region 5' to the mRNA start site.