Pancreatic-specific expression of the glucose transporter type 2 gene: identification of cis-elements and islet-specific trans-acting factors.

A defect in glucose sensing of the pancreatic beta-cells has been observed in several animal models of type II diabetes and has been correlated with a reduced gene expression of the glucose transporter type 2 (Glut2). In a transgenic mouse model, expression of Glut2 antisense RNA in pancreatic beta-cells has recently been shown to be associated with an impaired glucose-induced insulin secretion and the development of diabetes. To identify factors that may be involved in the specific decrease of Glut2 in the beta-cells of the diabetic animal, an attempt was made to localize the cis-elements and trans-acting factors involved in the control of Glut2 expression in the endocrine pancreas. It was demonstrated by transient transfection studies that only 338 base pairs (bp) of the murine Glut2 proximal promoter are needed for reporter gene expression in pancreatic islet-derived cell lines, whereas no activity was detected in nonpancreatic cells. Three cis-elements, GTI, GTII, and GTIII, have been identified by DNAse I footprinting and gel retardation experiments within these 338 bp. GTI and GTIII bind distinct but ubiquitously expressed trans-acting factors. On the other hand, nuclear proteins specifically expressed in pancreatic cell lines interact with GTII, and their relative abundance correlates with endogenous Glut2 expression. These GTII-binding factors correspond to nuclear proteins of 180 and 90 kilodaltons as defined by Southwestern analysis. The 180-kilodalton factor is present in pancreatic beta-cell lines but not in an alpha-cell line. Mutation of the GTI or GTIII cis-elements decreases transcriptional activity directed by the 338-bp promoter, whereas mutation of GTII increases gene transcription. Thus negative and positive regulatory sequences are identified within the proximal 338 bp of the GLUT2 promoter and may participate in the islet-specific expression of the gene by binding beta-cell specific trans-acting factors.

Mutation screening of the urocortin gene: identification of new single nucleotide polymorphisms and association studies with obesity in French Caucasians.

A linkage between obesity-related phenotypes and the 2p21-23 locus has been reported previously. The urocortin (UCN) gene resides at this interval, and its protein decreases appetite behavior, suggesting that UCN may be a candidate gene for susceptibility to obesity. We localized the UCN gene by radiation hybrid mapping, and the surrounding markers were genotyped in a collection of French families. Evidence for linkage was shown between the marker D2S165 and leptin levels (LOD score, 1.34; P = 0.006) and between D2S2247 and the z-score of body mass index (LOD score, 1.829; P = 0.0019). The gene was screened for SNPs in 96 obese patients. Four new variants were established. Two single nucleotide polymorphisms were located in the promoter (-535 A-->G, -286 G-->A), one in intron 1 (+31 C-->G), and one in the 3'-untranslated region (+34 C-->T). Association studies in cohorts of 722 unrelated obese and 381 control subjects and transmission disequilibrium tests, performed for the two frequent promoter polymorphisms, in 120 families (894 individuals) showed that no association was present between these variants and obesity, obesity-related phenotypes, and diabetes. Thus, our analyses of the genetic variations of the UCN gene suggest that, at least in French Caucasians, they do not represent a major cause of obesity.

Hepatocyte nuclear factor-4 alpha regulates the human apolipoprotein AV gene: Identification of a novel response element and involvement in the control by peroxisome proliferator-activated receptor-gamma coactivator-1 alpha, AMP-activated protein kinase, and mitogen-activated protein kinase pathway

The recently discovered apolipoprotein AV (apoAV) gene has been reported to be a key player in modulating plasma triglyceride levels. Here we identify the hepatocyte nuclear factor-4 (HNF-4 ) as a novel regulator of human apoAV gene. Inhibition of HNF-4 expression by small interfering RNA resulted in down-regulation of apoAV. Deletion, mutagenesis, and binding assays revealed that HNF-4 directly regulates human apoAV promoter through DR1 [a direct repeat separated by one nucleotide (nt)], and via a novel element for HNF-4 consisting of an inverted repeat separated by 8 nt (IR8). In addition, we show that the coactivator peroxisome proliferator-activated receptor- coactivator-1 was capable of stimulating the HNF-4 -dependent transactivation of apoAV promoter. Furthermore, analyses in human hepatic cells demonstrated that AMP-activated protein kinase (AMPK) and the MAPK signaling pathway regulate human apoAV expression and suggested that this regulation may be mediated, at least in part, by changes in HNF-4 . Intriguingly, EMSAs and mice with a liver-specific disruption of the HNF-4 gene revealed a species-distinct regulation of apoAV by HNF-4 , which resembles that of a subset of HNF-4 target genes. Taken together, our data provide new insights into the binding properties and the modulation of HNF-4 and underscore the role of HNF-4 in regulating triglyceride metabolism.

Candidate gene identification following linkage: Search for hypertension susceptibility genes

Following up genetic linkage studies to identify the underlying susceptibility gene(s) for complex disease traits is an arduous yet biologically and clinically important task. Complex traits, such as hypertension, are considered polygenic with many genes influencing risk, each with small effects. Chromosome 2 has been consistently identified as a genomic region with genetic linkage evidence suggesting that one or more loci contribute to blood pressure levels and hypertension status. Using combined positional candidate gene methods, the Family Blood Pressure Program has concentrated efforts in investigating this region of chromosome 2 in an effort to identify underlying candidate hypertension susceptibility gene(s). Initial informatics efforts identified the boundaries of the region and the known genes within it. A total of 82 polymorphic sites in eight positional candidate genes were genotyped in a large hypothesis-generating sample consisting of 1640 African Americans, 1339 whites, and 1616 Mexican Americans. To adjust for multiple comparisons, resampling-based false discovery adjustment was applied, extending traditional resampling methods to sibship samples. Following this adjustment for multiple comparisons, SLC4A5, a sodium bicarbonate transporter, was identified as a primary candidate gene for hypertension. Polymorphisms in SLC4A5 were subsequently genotyped and analyzed for validation in two populations of African Americans (N = 461; N = 778) and two of whites (N = 550; N = 967). Again, SNPs within SLC4A5 were significantly associated with blood pressure levels and hypertension status. While not identifying a single causal DNA sequence variation that is significantly associated with blood pressure levels and hypertension status across all samples, the results further implicate SLC4A5 as a candidate hypertension susceptibility gene, validating previous evidence for one or more genes on chromosome 2 that influence hypertension related phenotypes in the population-at-large. The methodology and results reported provide a case study of one approach for following up the results of genetic linkage analyses to identify genes influencing complex traits. ^

Generation of longer cDNA fragments from serial analysis of gene expression tags for gene identification

We have developed a technique called the generation of longer cDNA fragments from serial analysis of gene expression (SAGE) tags for gene identification (GLGI), to convert SAGE tags of 10 bases into their corresponding 3′ cDNA fragments covering hundred bases. A primer containing the 10-base SAGE tag is used as the sense primer, and a single base anchored oligo(dT) primer is used as an antisense primer in PCR, together with Pfu DNA polymerase. By using this approach, a cDNA fragment extending from the SAGE tag toward the 3′ end of the corresponding sequence can be generated. Application of the GLGI technique can solve two critical issues in applying the SAGE technique: one is that a longer fragment corresponding to a SAGE tag, which has no match in databases, can be generated for further studies; the other is that the specific fragment corresponding to a SAGE tag can be identified from multiple sequences that match the same SAGE tag. The development of the GLGI method provides several potential applications. First, it provides a strategy for even wider application of the SAGE technique for quantitative analysis of global gene expression. Second, a combined application of SAGE/GLGI can be used to complete the catalogue of the expressed genes in human and in other eukaryotic species. Third, it can be used to identify the 3′ cDNA sequence from any exon within a gene. It can also be used to confirm the reality of exons predicted by bioinformatic tools in genomic sequences. Fourth, a combined application of SAGE/GLGI can be applied to define the 3′ boundary of expressed genes in the genomic sequences in human and in other eukaryotic genomes.

Mutagenesis and gene identification in Dictyostelium by shotgun antisense.

We have developed a mutagenesis technique that uses antisense cDNA to identify genes required for development in Dictyostelium discoideum. We transformed Dictyostelium cells with a cDNA library made from the mRNA of vegetative and developing cells. The cDNA was cloned in an antisense orientation immediately downstream of a vegetative promoter, so that in transformed cells the promoter will drive the synthesis of an antisense RNA transcript. We find that individual transformants typically contain one or occasionally two antisense cDNAs. Using this mutagenesis technique, we have generated mutants that fail to aggregate, aggregate but fail to form fruiting bodies, or aggregate but form abnormal fruiting bodies. The individual cDNA molecules from the mutants were identified and cloned using PCR. Initial sequence analysis of the PCR products from 35 mutants has identified six novel Dictyostelium genes, each from a transformant with one antisense cDNA. When the PCR-isolated antisense cDNAs were ligated into the antisense vector and the resulting constructs transformed into cells, the phenotypes of the transformed cells matched those of the original mutants from which each cDNA was obtained. We made homologous recombinant gene disruption transformants for three of the novel genes, in each case generating mutants with phenotypes indistinguishable from those of the original antisense transformants. Shotgun antisense thus is a rapid way to identify genes in Dictyostelium and possibly other organisms.

Sequencing and identification of expressed Schistosoma mansoni genes by random selection of cDNA clones from a directional library

We have initiated a gene discovery program in Schistosoma mansoni based on the technique of Expressed Sequence Tags (ESTs), i.e. partial sequences of cDNAs obtained from single passes in automatic DNA sequencers. ESTs can be used to identify genese onf the basis of their homology whith sequences from other species deposited in DNA or protein databases. Trasncripts with sequences without matches in teh databases may represent novel parasite-specific genes. This approach has shown to be very efficient and in less than two years a broad range of novel genes has already been ascertained, more than doubling the number of known S. mansoni genes.

An assessment of gene prediction accuracy in large DNA sequences

One of the first useful products from the human genome will be a set of predicted genes. Besides its intrinsic scientific interest, the accuracy and completeness of this data set is of considerable importance for human health and medicine. Though progress has been made on computational gene identification in terms of both methods and accuracy evaluation measures, most of the sequence sets in which the programs are tested are short genomic sequences, and there is concern that these accuracy measures may not extrapolate well to larger, more challenging data sets. Given the absence of experimentally verified large genomic data sets, we constructed a semiartificial test set comprising a number of short single-gene genomic sequences with randomly generated intergenic regions. This test set, which should still present an easier problem than real human genomic sequence, mimics the approximately 200kb long BACs being sequenced. In our experiments with these longer genomic sequences, the accuracy of GENSCAN, one of the most accurate ab initio gene prediction programs, dropped significantly, although its sensitivity remained high. Conversely, the accuracy of similarity-based programs, such as GENEWISE, PROCRUSTES, and BLASTX was not affected significantly by the presence of random intergenic sequence, but depended on the strength of the similarity to the protein homolog. As expected, the accuracy dropped if the models were built using more distant homologs, and we were able to quantitatively estimate this decline. However, the specificities of these techniques are still rather good even when the similarity is weak, which is a desirable characteristic for driving expensive follow-up experiments. Our experiments suggest that though gene prediction will improve with every new protein that is discovered and through improvements in the current set of tools, we still have a long way to go before we can decipher the precise exonic structure of every gene in the human genome using purely computational methodology.

Identification of canine papillomavirus type 1 (CPV1) DNA in dogs with cutaneous papillomatosis

Canine oral papillomavirus (COPV), also known as Canine Papillomavirus type 1 (CPV1), induces papillomas at the mucous membranes of the oral cavity and at the haired skin of dogs. The classification of Papillomavirus (PV) types is based on the L1 capsid protein and nucleotide sequence; so far, 14 CPV types have been described in several countries, but the molecular characterization of CPV in Brazil is lacking. This study investigated the presence of the PV in seven papillomas from four mixed breed dogs from Londrina/PR, Southern Brazil, by partial sequencing of the L1 gene. Seven exophytic cutaneous lesions were surgically removed and processed for histopathological and molecular characterization. Histopathology confirmed the lesions as viral papillomas due to typical histological features. Polymerase Chain Reaction (PCR) assay using the FAP59 and FAP64 primers targeted the L1 gene followed by sequence analysis of the amplicons identified CPV1 in all evaluated papilloma samples. This study represents the first description of CPV1 DNA associated with canine papillomatosis in Brazil.

Skeletal muscle protein content of lipolytic inhibitor G(0)/G(1) switch gene-2 protein: the effect of endurance training

The first and rate-limiting step of lipolysis is the removal of the first fatty acid from a triglyceride molecule; it is catalyzed by adipose triglyceride lipase (ATGL). ATGL is co-activated by comparative gene identification-58 (CGI-58) and inhibited by the G(0)/G(1) switch gene-2 protein (G0S2). G0S2 has also recently been identified as a positive regulator of oxidative phosphorylation within the mitochondria. Previous research has demonstrated in cell culture, a dose dependent mechanism for inhibition by G0S2 on ATGL. However our data is not consistent with this hypothesis. There was no change in G0S2 protein content during an acute lipolytic inducing set of contractions in both whole muscle, and isolated mitochondria yet both ATGL and G0S2 increase following endurance training, in spite of the fact that there should be increased reliance on intramuscular lipolysis. Therefore, inhibition of ATGL by G0S2 appears to be regulated through more complicated intracellular or post-translation regulation.

Regulation of gene expression in the dinoflagellate Lingulodinium polyedrum

Les dinoflagellés sont des eucaryotes unicellulaires que l’on retrouve autant en eau douce qu’en milieu marin. Ils sont particulièrement connus pour causer des fleurs d’algues toxiques nommées ‘marée-rouge’, ainsi que pour leur symbiose avec les coraux et pour leur importante contribution à la fixation du carbone dans les océans. Au point de vue moléculaire, ils sont aussi connus pour leur caractéristiques nucléaires uniques, car on retrouve généralement une quantité immense d’ADN dans leurs chromosomes et ceux-ci sont empaquetés et condensés sous une forme cristalline liquide au lieu de nucléosomes. Les gènes encodés par le noyau sont souvent présents en multiples copies et arrangés en tandem et aucun élément de régulation transcriptionnelle, y compris la boite TATA, n’a encore été observé. L’organisation unique de la chromatine des dinoflagellés suggère que différentes stratégies sont nécessaires pour contrôler l’expression des gènes de ces organismes. Dans cette étude, j’ai abordé ce problème en utilisant le dinoflagellé photosynthétique Lingulodinium polyedrum comme modèle. L. polyedrum est d’un intérêt particulier, car il a plusieurs rythmes circadiens (journalier). À ce jour, toutes les études sur l’expression des gènes lors des changements circadiens ont démontrées une régulation à un niveau traductionnel. Pour mes recherches, j’ai utilisé les approches transcriptomique, protéomique et phosphoprotéomique ainsi que des études biochimiques pour donner un aperçu de la mécanique de la régulation des gènes des dinoflagellés, ceci en mettant l’accent sur l’importance de la phosphorylation du système circadien de L. polyedrum. L’absence des protéines histones et des nucléosomes est une particularité des dinoflagellés. En utilisant la technologie RNA-Seq, j’ai trouvé des séquences complètes encodant des histones et des enzymes modifiant les histones. L polyedrum exprime donc des séquences conservées codantes pour les histones, mais le niveau d’expression protéique est plus faible que les limites de détection par immunodétection de type Western. Les données de séquençage RNA-Seq ont également été utilisées pour générer un transcriptome, qui est une liste des gènes exprimés par L. polyedrum. Une recherche par homologie de séquences a d’abord été effectuée pour classifier les transcrits en diverses catégories (Gene Ontology; GO). Cette analyse a révélé une faible abondance des facteurs de transcription et une surprenante prédominance, parmi ceux-ci, des séquences à domaine Cold Shock. Chez L. polyedrum, plusieurs gènes sont répétés en tandem. Un alignement des séquences obtenues par RNA-Seq avec les copies génomiques de gènes organisés en tandem a été réalisé pour examiner la présence de transcrits polycistroniques, une hypothèse formulée pour expliquer le manque d’élément promoteur dans la région intergénique de la séquence de ces gènes. Cette analyse a également démontré une très haute conservation des séquences codantes des gènes organisés en tandem. Le transcriptome a également été utilisé pour aider à l’identification de protéines après leur séquençage par spectrométrie de masse, et une fraction enrichie en phosphoprotéines a été déterminée comme particulièrement bien adapté aux approches d’analyse à haut débit. La comparaison des phosphoprotéomes provenant de deux périodes différentes de la journée a révélée qu’une grande partie des protéines pour lesquelles l’état de phosphorylation varie avec le temps est reliées aux catégories de liaison à l’ARN et de la traduction. Le transcriptome a aussi été utilisé pour définir le spectre des kinases présentes chez L. polyedrum, qui a ensuite été utilisé pour classifier les différents peptides phosphorylés qui sont potentiellement les cibles de ces kinases. Plusieurs peptides identifiés comme étant phosphorylés par la Casein Kinase 2 (CK2), une kinase connue pour être impliquée dans l’horloge circadienne des eucaryotes, proviennent de diverses protéines de liaison à l’ARN. Pour évaluer la possibilité que quelques-unes des multiples protéines à domaine Cold Shock identifiées dans le transcriptome puissent moduler l’expression des gènes de L. polyedrum, tel qu’observé chez plusieurs autres systèmes procaryotiques et eucaryotiques, la réponse des cellules à des températures froides a été examinée. Les températures froides ont permis d’induire rapidement un enkystement, condition dans laquelle ces cellules deviennent métaboliquement inactives afin de résister aux conditions environnementales défavorables. Les changements dans le profil des phosphoprotéines seraient le facteur majeur causant la formation de kystes. Les phosphosites prédits pour être phosphorylés par la CK2 sont la classe la plus fortement réduite dans les kystes, une découverte intéressante, car le rythme de la bioluminescence confirme que l’horloge a été arrêtée dans le kyste.

Role of the Protein Tyrosine Kinase 7 gene in human neural tube defects

Les anomalies du tube neural (ATN) sont des anomalies développementales où le tube neural reste ouvert (1-2/1000 naissances). Afin de prévenir cette maladie, une connaissance accrue des processus moléculaires est nécessaire. L’étiologie des ATN est complexe et implique des facteurs génétiques et environnementaux. La supplémentation en acide folique est reconnue pour diminuer les risques de développer une ATN de 50-70% et cette diminution varie en fonction du début de la supplémentation et de l’origine démographique. Les gènes impliqués dans les ATN sont largement inconnus. Les études génétiques sur les ATN chez l’humain se sont concentrées sur les gènes de la voie métabolique des folates du à leur rôle protecteur dans les ATN et les gènes candidats inférés des souris modèles. Ces derniers ont montré une forte association entre la voie non-canonique Wnt/polarité cellulaire planaire (PCP) et les ATN. Le gène Protein Tyrosine Kinase 7 est un membre de cette voie qui cause l’ATN sévère de la craniorachischisis chez les souris mutantes. Ptk7 interagit génétiquement avec Vangl2 (un autre gène de la voie PCP), où les doubles hétérozygotes montrent une spina bifida. Ces données font de PTK7 comme un excellent candidat pour les ATN chez l’humain. Nous avons re-séquencé la région codante et les jonctions intron-exon de ce gène dans une cohorte de 473 patients atteints de plusieurs types d’ATN. Nous avons identifié 6 mutations rares (fréquence allélique <1%) faux-sens présentes chez 1.1% de notre cohorte, dont 3 sont absentes dans les bases de données publiques. Une variante, p.Gly348Ser, a agi comme un allèle hypermorphique lorsqu'elle est surexprimée dans le modèle de poisson zèbre. Nos résultats impliquent la mutation de PTK7 comme un facteur de risque pour les ATN et supporte l'idée d'un rôle pathogène de la signalisation PCP dans ces malformations.

Identity by descent and candidate gene mapping of Richieri-Costa and Pereira syndrome

The Richieri-Costa-Pereira syndrome is a rare autosomal recessive disorder characterized by short stature, Robin sequence, cleft mandible, pre/postaxial anomalies and clubfoot. of 15 families reported with this disorder 14 are from Brazil suggesting a founder effect. We studied 15 families using identity-by-descent as a hypothesis to attempt gene localization We have examined through linkage analysis 497 polymorphicmarkers and also performed direct sequencing of exons for 10 candidate genes selected on the basis of their expression in the developing mandible and limb. No evidence for allele sharing at any locus tested or mutations in candidate genes was found. Additional higher resolution mapping, new families and other candidate genes might improve future chances of gene identification. (C) 2003 Wiley-Liss, Inc.

Comparative analysis of noncoding sequences of orthologous bovine and human gene pairs

Genomic sequence comparison across species has enabled the elucidation of important coding and regulatory sequences encoded within DNA. Of particular interest are the noncoding regulatory sequences, which influence gene transcriptional and posttranscriptional processes. A phylogenetic footprinting strategy was employed to identify noncoding conservation patterns of 39 human and bovine orthologous genes. Seventy-three conserved noncoding sequences were identified that shared greater than 70% identity over at least 100 bp. Thirteen of these conserved sequences were also identified in the mouse genome. Evolutionary conservation of noncoding sequences across diverse species may have functional significance, and these conserved sequences may be good candidates for regulatory elements.

Identification and complete sequencing of novel human transcripts through the use of mouse orthologs and testis cDNA sequences

The correct identification of all human genes, and their derived transcripts, has not yet been achieved, and it remains one of the major aims of the worldwide genomics community. Computational programs suggest the existence of 30,000 to 40,000 human genes. However, definitive gene identification can only be achieved by experimental approaches. We used two distinct methodologies, one based on the alignment of mouse orthologous sequences to the human genome, and another based on the construction of a high-quality human testis cDNA library, in an attempt to identify new human transcripts within the human genome sequence. We generated 47 complete human transcript sequences, comprising 27 unannotated and 20 annotated sequences. Eight of these transcripts are variants of previously known genes. These transcripts were characterized according to size, number of exons, and chromosomal localization, and a search for protein domains was undertaken based on their putative open reading frames. In silico expression analysis suggests that some of these transcripts are expressed at low levels and in a restricted set of tissues.