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.

Widespread occurrence of non-canonical transcription termination by human RNA polymerase III.

Human RNA polymerase (Pol) III-transcribed genes are thought to share a simple termination signal constituted by four or more consecutive thymidine residues in the coding DNA strand, just downstream of the RNA 3'-end sequence. We found that a large set of human tRNA genes (tDNAs) do not display any T(≥4) stretch within 50 bp of 3'-flanking region. In vitro analysis of tDNAs with a distanced T(≥4) revealed the existence of non-canonical terminators resembling degenerate T(≥5) elements, which ensure significant termination but at the same time allow for the production of Pol III read-through pre-tRNAs with unusually long 3' trailers. A panel of such non-canonical signals was found to direct transcription termination of unusual Pol III-synthesized viral pre-miRNA transcripts in gammaherpesvirus 68-infected cells. Genome-wide location analysis revealed that human Pol III tends to trespass into the 3'-flanking regions of tDNAs, as expected from extensive terminator read-through. The widespread occurrence of partial termination suggests that the Pol III primary transcriptome in mammals is unexpectedly enriched in 3'-trailer sequences with the potential to contribute novel functional ncRNAs.

Ribonucleotide reductase genes of Bacillus prophages: a refuge to introns and intein coding sequences.

The ribonucleotide reductase gene tandem bnrdE/bnrdF in SPbeta-related prophages of different Bacillus spp. isolates presents different configurations of intervening sequences, comprising one to three of six non-homologous splicing elements. Insertion sites of group I introns and intein DNA are clustered in three relatively short segments encoding functionally important domains of the ribonucleotide reductase. Comparison of the bnrdE homologs reveals mutual exclusion of a group I intron and an intein coding sequence flanking the codon that specifies a conserved cysteine. In vivo splicing was demonstrated for all introns. However, for two of them a part of the mRNA precursor molecules remains unspliced. Intergenic bnrdE-bnrdF regions are unexpectedly long, comprising between 238 and 541 nt. The longest encodes a putative polypeptide related to HNH homing endonucleases.

TcruziDB, an Integrated Database, and the WWW Information Server for the Trypanosoma cruzi Genome Project

Data analysis, presentation and distribution is of utmost importance to a genome project. A public domain software, ACeDB, has been chosen as the common basis for parasite genome databases, and a first release of TcruziDB, the Trypanosoma cruzi genome database, is available by ftp from ftp://iris.dbbm.fiocruz.br/pub/genomedb/TcruziDB as well as versions of the software for different operating systems (ftp://iris.dbbm.fiocruz.br/pub/unixsoft/). Moreover, data originated from the project are available from the WWW server at http://www.dbbm.fiocruz.br. It contains biological and parasitological data on CL Brener, its karyotype, all available T. cruzi sequences from Genbank, data on the EST-sequencing project and on available libraries, a T. cruzi codon table and a listing of activities and participating groups in the genome project, as well as meeting reports. T. cruzi discussion lists (tcruzi-l@iris.dbbm.fiocruz.br and tcgenics@iris.dbbm.fiocruz.br) are being maintained for communication and to promote collaboration in the genome project

Molecular genetics of "PRPF31" dominant mutations linked to retinitis pigmentosa

ABSTRACT : The retina is one of the most important human sensory tissues since it detects and transmits all visual information from the outside world to the brain. Retinitis pigmentosa (RP) is the name given to a group of inherited diseases that affect specifically the photoreceptors present in the retina and in many instances lead to blindness. Dominant mutations in PRPF31, a gene that encodes for a pre-mRNA splicing factor, cause retinitis pigmentosa with reduced penetrance. We functionally investigated a novel mutation, identified in a large family with autosomal dominant RP, and 7 other mutations, substitutions and microdeletions, in 12 patients from 7 families with PRPF31-linked RP. Seven mutations lead to PRPF31 mRNA with premature stop codons and one to mRNA lacking the exon containing the initiation codon. Quantification of PRPF31 mRNA and protein levels revealed a significant reduction in cell lines derived from patients, compared to non carriers of mutations in PRPF31. Allelic quantification of PRPF31 mRNA indicated that the level of mutated mRNA is very low compared to wild-type mRNA. No mutant protein was detected and the subnuclear localization of wild-type PRPF31 remains the same in cell lines from patients and controls. Blocking nonsense-mediated mRNA decay in cell lines derived from patients partially restored PRPF31 mutated mRNA but derived proteins were still undetectable, even when protein degradation pathways were inhibited. Our results demonstrated that the vast majority of PRPF31 mutations result in null alleles, since they are subject to surveillance mechanisms that degrade mutated mRNA and possibly block its translation. Altogether, these data indicate that the likely cause of PRPF31-linked RP is haploinsufficiency, rather than a dominant negative effect. Penetrance of PRPF31 mutations has been previously demonstrated to be inversely correlated with the level of PRPF31 mRNA, since high expression of wild-type PRPF31 mRNA protects from the disease. Consequently, we have investigated the genetic modifiers that control the expression of PRPF31 by quantifying PRPF31 mRNA levels in cell lines derived from 200 individuals from 15 families representative of the general population. By linkage analyses we identified a 8.2Mb-region on chromosome 14q21-23 that contains a gene involved in the modulation of PRPF31 expression. We also assessed apreviously-mapped penetrance factor invariably located on the wild-type allele and linked to the PRPF31 locus in asymptomatic patients from different families with RP. We demonstrated that this modifier increases the expression of both PRPF31 alleles already at the pre-mRNA level. Finally, our data suggest that PRPF31 mRNA expression and consequently the penetrance of PRPF31 mutations is modulated by at least 2 diffusible compounds, which act on both PRPF31 alleles during their transcription.

Improving the performance of positive selection inference by filtering unreliable alignment regions.

Errors in the inferred multiple sequence alignment may lead to false prediction of positive selection. Recently, methods for detecting unreliable alignment regions were developed and were shown to accurately identify incorrectly aligned regions. While removing unreliable alignment regions is expected to increase the accuracy of positive selection inference, such filtering may also significantly decrease the power of the test, as positively selected regions are fast evolving, and those same regions are often those that are difficult to align. Here, we used realistic simulations that mimic sequence evolution of HIV-1 genes to test the hypothesis that the performance of positive selection inference using codon models can be improved by removing unreliable alignment regions. Our study shows that the benefit of removing unreliable regions exceeds the loss of power due to the removal of some of the true positively selected sites.

Mutations in Plasmodium falciparum dihydrofolate reductase and dihydropteroate synthase of isolates from the Amazon Region of Brazil

Since the late 1970s pyrimethamine-sulfadoxine (PS; FansidarTM Hoffman-LaRoche, Basel) has been used as first line therapy for uncomplicated malaria in the Amazon basin. Unfortunately, resistance has developed over the last ten years in many regions of the Amazon and PS is no longer recommended for use in Brazil. In vitro resistance to pyrimethamine and cycloguanil (the active metabolite of proguanil) is caused by specific point mutations in Plasmodium falciparum dihydrofolate reductase (DHFR), and in vitro resistance to sulfadoxine has been associated with mutations in dihydropteroate synthase (DHPS). In association with a proguanil-sulfamethoxazole clinical trial in Brazil, we performed a nested mutation-specific polymerase chain reaction to measure the prevalence of DHFR mutations at codons 50, 51, 59, 108 and 164 and DHPS mutations at codons 436, 437, 540, 581 and 613 at three sites in the Brazilian Amazon. Samples from two isolated towns showed a high degree of homogeneity, with the DHFR Arg-50/Ile-51/Asn-108 and DHPS Gly-437/Glu-540/Gly-581 mutant genotype accounting for all infections in Peixoto de Azevedo (n = 15) and 60% of infections in Apiacás (n = 10), State of Mato Grosso. The remaining infections in Apiacás differed from this predominant genotype only by the addition of the Bolivia repeat at codon 30 and the Leu-164 mutation in DHFR. By contrast, 17 samples from Porto Velho, capital city of the State of Rondônia, with much in- and out-migration, showed a wide variety of DHFR and DHPS genotypes.

Protein-binding site at the immunoglobulin mu membrane polyadenylylation signal: possible role in transcription termination.

mRNAs specifying immunoglobulin mu and delta heavy chains are encoded by a single large, complex transcription unit (mu + delta gene). The transcriptional activity of delta gene segments in terminally differentiated, IgM-secreting B lymphocytes is 10-20 times lower than in earlier B-lineage cells expressing delta mRNA. We find that transcription of the mu + delta gene in IgM-secreting murine myeloma cells terminates within a region of 500-1000 nucleotides immediately following the mu membrane (mu m) polyadenylylation site. Transcription decreases only minimally through this region in murine cell lines representative of earlier stages in B-cell development. A DNA fragment containing the mu m polyadenylylation signal gives protein-DNA complexes with different mobilities in gel retardation assays with nuclear extracts from myeloma cells than with nuclear extracts from earlier B-lineage cells. However, using a recently developed "footprinting" procedure in which protein-DNA complexes resolved in gel retardation assays are subjected to nucleolytic cleavage while still in the polyacrylamide gel, we find that the DNA sequences protected by factors from the two cell types are indistinguishable. The factor-binding site on the DNA is located 5' of the mu m polyadenylylation signal AATAAA and includes the 15-nucleotide-long A + T-rich palindrome CTGTAAACAAATGTC. This type of palindromic binding site exhibits orientation-dependent activity consistent with the reported properties of polymerase II termination signals. This binding site is followed by two sets of directly repeated DNA sequences with different helical conformation as revealed by their reactivity with the chemical nuclease 1,10-phenanthroline-copper. The close proximity of these features to the signals for mu m mRNA processing may reflect a linkage of the processes of developmentally regulated mu m polyadenylylation and transcription termination.

An mRNA region of the canine distemper virus fusion protein gene lacking AUG codons can promote protein expression.

Canine distemper virus (CDV) produces a glycosylated type I fusion protein (F) with an internal hydrophobic signal sequence beginning around 115 residues downstream of the first AUG used for translation initiation. Cleavage of the signal sequence yields the F0 molecule, which is cleaved into the F1 and F2 subunits. Surprisingly, when all in-frame AUGs located in the first third of the F gene were mutated a protein of the same molecular size as the F0 molecule was still expressed from both the Onderstepoort (OP) and A75/17-CDV F genes. We designated this protein, which is initiated from a non-AUG codon protein Fx. Site-directed mutagenesis allowed to identify codon 85, a GCC codon coding for alanine, as the most likely position from which translation initiation of Fx occurs in OP-CDV. Deletion analysis demonstrated that at least 60 nucleotides upstream of the GCC codon are required for efficient Fx translation. This sequence is GC-rich, suggesting extensive folding. Secondary structure may therefore be important for translation initiation at codon 85.

Characterization of the murine high Km glucose transporter GLUT2 gene and its transcriptional regulation by glucose in a differentiated insulin-secreting cell line.

In pancreatic beta-cells, the high Km glucose transporter GLUT2 catalyzes the first step in glucose-induced insulin secretion by glucose uptake. Expression of the transporter has been reported to be modulated by glucose either at the protein or mRNA levels. In this study we used the differentiated insulinoma cell line INS-1 which expresses high levels of GLUT2 and show that the expression of GLUT2 is regulated by glucose at the transcriptional level. By run-on transcription assays we showed that glucose induced GLUT2 gene transcription 3-4-fold in INS-1 cells which was paralleled by a 1.7-2.3-fold increase in cytoplasmic GLUT2 mRNA levels. To determine whether glucose regulatory sequences were present in the promoter region of GLUT2, we cloned and characterized a 1.4-kilobase region of mouse genomic DNA located 5' of the translation initiation site. By RNase protection assays and primer extension, we determined that multiple transcription initiation sites were present at positions -55, -64, and -115 from the first coding ATG and which were identified in liver, intestine, kidney, and beta-cells mRNAs. Plasmids were constructed with the mouse promoter region linked to the reporter gene chloramphenicol acetyltransferase (CAT), and transiently and stably transfected in the INS-1 cells. Glucose induced a concentration-dependent increase in CAT activity which reached a maximum of 3.6-fold at 20 mM glucose. Similar CAT constructs made of the human GLUT2 promoter region and the CAT gene displayed the same glucose-dependent increase in transcriptional activity when transfected into INS-1 cells. Comparison of the mouse and human promoter regions revealed sequence identity restricted to a few stretches of sequences which suggests that the glucose responsive element(s) may be conserved in these common sequences.

Discovery of amino acid variants in the human glucose-dependent insulinotropic polypeptide (GIP) receptor: the impact on the pancreatic beta cell responses and functional expression studies in Chinese hamster fibroblast cells.

The two incretins, glucose-dependent insulinotropic polypeptide (GIP) and glucagon-like peptide-1 (GLP-1), are insulinotropic factors released from the small intestine to the blood stream in response to oral glucose ingestion. The insulinotropic effect of GLP-1 is maintained in patients with Type II (non-insulin-dependent) diabetes mellitus, whereas, for unknown reasons, the effect of GIP is diminished or lacking. We defined the exon-intron boundaries of the human GIP receptor, made a mutational analysis of the gene and identified two amino acid substitutions, A207 V and E354Q. In an association study of 227 Caucasian Type II diabetic patients and 224 matched glucose tolerant control subjects, the allelic frequency of the A207 V polymorphism was 1.1% in Type II diabetic patients and 0.7% in control subjects (p = 0.48), whereas the allelic frequency of the codon 354 polymorphism was 24.9% in Type II diabetic patients versus 23.2% in control subjects. Interestingly, the glucose tolerant subjects (6% of the population) who were homozygous for the codon 354 variant had on average a 14% decrease in fasting serum C-peptide concentration (p = 0.01) and an 11% decrease in the same variable 30 min after an oral glucose load (p = 0.03) compared with subjects with the wild-type receptor. Investigation of the function of the two GIP receptor variants in Chinese hamster fibroblasts showed, however, that the GIP-induced cAMP formation and the binding of GIP to cells expressing the variant receptors were not different from the findings in cells expressing the wildtype GIP receptor. In conclusion, amino acid variants in the GIP receptor are not associated with random Type II diabetes in patients of Danish Caucasian origin or with altered GIP binding and GIP-induced cAMP production when stably transfected in Chinese hamster fibroblasts. The finding of an association between homozygosity for the codon 354 variant and reduced fasting and post oral glucose tolerance test (OGTT) serum C-peptide concentrations, however, calls for further investigations and could suggest that GIP even in the fasting state regulates the beta-cell secretory response.