Candidate genes for obesity-susceptibility show enriched association within a large genome-wide association study for BMI

Before the advent of genome-wide association studies (GWASs), hundreds of candidate genes for obesity-susceptibility had been identified through a variety of approaches. We examined whether those obesity candidate genes are enriched for associations with body mass index (BMI) compared with non-candidate genes by using data from a large-scale GWAS. A thorough literature search identified 547 candidate genes for obesity-susceptibility based on evidence from animal studies, Mendelian syndromes, linkage studies, genetic association studies and expression studies. Genomic regions were defined to include the genes ±10 kb of flanking sequence around candidate and non-candidate genes. We used summary statistics publicly available from the discovery stage of the genome-wide meta-analysis for BMI performed by the genetic investigation of anthropometric traits consortium in 123 564 individuals. Hypergeometric, rank tail-strength and gene-set enrichment analysis tests were used to test for the enrichment of association in candidate compared with non-candidate genes. The hypergeometric test of enrichment was not significant at the 5% P-value quantile (P = 0.35), but was nominally significant at the 25% quantile (P = 0.015). The rank tail-strength and gene-set enrichment tests were nominally significant for the full set of genes and borderline significant for the subset without SNPs at P < 10(-7). Taken together, the observed evidence for enrichment suggests that the candidate gene approach retains some value. However, the degree of enrichment is small despite the extensive number of candidate genes and the large sample size. Studies that focus on candidate genes have only slightly increased chances of detecting associations, and are likely to miss many true effects in non-candidate genes, at least for obesity-related traits.

Analysis of a human placental lactogen gene promoter

Human placental lactogen (hPL) and human growth hormone (hGH) comprise a multigene family that share $>$90% nucleic acid sequence homology including 500 bp of 5$\sp\prime$ flanking sequence. Despite these similarities, hGH is produced in the anterior pituitary while hPL is expressed in the placenta. For most genes studied to date, regulation of expression occurs by alterations at the level of transcriptional initiation. Nuclear proteins bind specific DNA sequences in the promoter to regulate gene expression. In this study, the hPL$\sb3$ promoter was analyzed for DNA sequences that contribute to its expression. The interaction between the hPL$\sb3$ promoter and nuclear proteins was examined using nuclear extracts from placental and non-placental cells.^ To identify regulatory elements in the promoter of the hPL$\sb3$ gene, 5$\sp\prime$ deletion mutants were constructed by cleaving 1200 bp of upstream sequence with various restriction enzymes. These DNA fragments were ligated 5$\sp\prime$ to a promoterless bacterial gene chloramphenicol acetyltransferase (CAT) and transfected into JEG-3 cells, a human placental choriocarcinoma cell line. The level of CAT activity reflects the ability of the promoter mutants to activate transcription. Deletion of the sequence between $-$142 bp and $-$129 bp, relative to the start of transcription, resulted in an 8-fold decrease in CAT activity. Nuclear proteins from JEG-3, HeLa, and HepG2 (human liver cells), formed specific binding complexes with this region of the hPL$\sb3$ promoter, as shown by gel mobility shift assay. The $-$142 bp to $-$129 bp region contains a sequence similar to that of a variant binding site for the transcription factor Sp1. Sp1-like proteins were identified by DNA binding assay, in the nuclear extracts of the three cell lines. A series of G nucleotides in the hPL$\sb3$ promoter regulatory region were identified by methylation interference assay to interact with the DNA-binding proteins and the pattern obtained is similar to that for other Sp1 binding sites that have been studied. This suggests that hPL$\sb3$ may be transcriptionally regulated by Sp1 or a Sp1-like transacting factor. ^

Regulation of {\it myogenin\/} transcription during mouse embryogenesis

The formation of skeletal muscle during vertebrate development involves the induction of mesoderm and subsequent generation of myoblasts that ultimately differentiate into mature muscles. The recent identification of a group of myogenic regulators that can convert fibroblasts to myoblasts has contributed to our understanding of the molecular events that underlie the establishment of the skeletal muscle phenotype. Members of this group of myogenic regulators share a helix-loop-helix (HLH) motif that mediates DNA binding. The myogenic HLH proteins bind to the consensus sequence CANNTG, referred to as an E-box, and activate muscle-specific transcription. In addition to E-boxes, other motifs, such as the MEF-2 binding site, have been shown to mediate muscle-specific transcription. The myogenic HLH proteins are expressed in the myogenic precursors in somites and limb buds, and in differentiated muscle fibers during embryogenesis, consistent with their roles as regulators for muscle development. The myogenic HLH proteins appear to auto-activate their own and cross-activate one another's expression in cultured cells. Myogenin is one of the myogenic HLH proteins and likely the regulator for terminal muscle differentiation. Myogenin is a common target of diverse regulatory pathways. To search for upstream regulators of myogenin, we studied regulation of myogenin transcription during mouse embryogenesis. We showed that the myogenin promoter contains a binding site for MEF-2, which can mediate indirectly the autoregulation of myogenin transcription. We found that a transgene under the control of a 1.5 kb 5$\sp\prime$ flanking sequence can recapitulate the temporal and spatial expression pattern of the endogenous myogenin gene during mouse embryogenesis. By tracing embryonic cells that activate myogenin-lacZ during embryogenesis, we found no evidence that lacZ was expressed in myogenic precursors migrating from somites to limb buds, suggesting the existence of regulators other than myogenic HLH proteins that can maintain cells in the myogenic lineage. Mutations of an E-box and a MEF-2 site in the myogenin promoter suppressed transcription in subsets of myogenic precursors in mouse embryos. These results suggest that myogenic HLH proteins and MEF-2 participate in separable regulatory pathways controlling myogenin transcription and provide evidence for positional regulation of myogenic regulators in the embryo. ^

Identification of two vitamin D response elements in the rat osteopontin gene

Studies to elucidate the function of vitamin D have demonstrated an important role in regulating bone-related cells, including osteoblasts and osteoclasts. A seemingly paradoxical observation is that 1,25(OH)$\sb2$D$\sb3$, the active metabolite of vitamin D, stimulates bone resorption, yet regulates transcription of genes expressed by osteoblasts. One mechanism that could explain these actions is the upregulation of transcription of osteoblast-specific genes. These gene products could then act as effectors to influence osteoclastic activity. We hypothesized that molecular signals could be deposited directly into the mineralized matrix in the form of noncollagenous proteins, such as osteopontin (OPN). The structure, biosynthesis and localization of OPN suggest that it could function to mediate the molecular "cross talk" between osteoblasts and osteoclasts in response to 1,25(OH)$\sb2$D$\sb3$. To begin to address this hypothesis, elucidation of the molecular mechanisms of action involved in the transactivation of OPN by 1,25(OH)$\sb2$D$\sb3$ is essential.^ In the present study, the rat opn gene was isolated and characterized. Functional analysis by transient transfection of the 5$\sp\prime$ flanking sequences of the rat opn gene fused to the luciferase gene demonstrated that OPN is transcriptionally upregulated by 1,25(OH)$\sb2$D$\sb3$, mediated through two vitamin D response elements (VDRE). Both proximal and distal VDREs are structurally similar (two imperfect direct repeats separated by a 3 nucleotide spacer) and bind protein complexes that include the VDR and retinoid-X receptor (RXR). Isolated VDRE expression constructs produce functional activity of equivalent magnitude of responsiveness to 1,25(OH)$\sb2$D$\sb3$. However, expression constructs containing either VDRE and at least 200 bp of 5$\sp\prime$ and 3$\sp\prime$ flanking sequence demonstrated that the distal VDRE produces an amplitude of response significantly higher than the proximal VDRE. We conclude that the transcriptional upregulation of the opn gene by 1,25(OH)$\sb2$D$\sb3$ involves the transactivation of two VDREs, while maximal responsiveness requires interaction of the VDREs with additional cis-elements contained in the 5$\sp\prime$ sequence. ^

NFkappaB and STAT binding elements participate in regulation of MUC1 expression in normal and transformed mammary epithelial cells

The MUC1 gene encodes a transmembrane mucin glycoprotein that is overexpressed in several cancers of epithelial origin, including those of breast, pancreas, lung, ovary, and colon. Functions of MUC1 include protection of mucosal epithelium, modulation of cellular adhesion, and signal transduction. Aberrantly increased expression of MUC1 in cancer cells promotes tumor progression through adaptation of these functions. Some regulatory elements participating in MUC1 transcription have been described, but the mechanisms responsible for overexpression are largely unknown. A region of MUC1 5′ flanking sequence containing two conserved potential cytokine response elements, an NFκB site at −589/−580 and a STAT binding element (SBE) at −503/−495, has been implicated in high level expression in breast and pancreatic cancer cell lines. Persistent stimulation by proinflammatory cytokines may contribute to increased MUC1 transcription by tumor cells. ^ T47D breast cancer cells and normal human mammary epithelial cells (HMEC) were used to determine the roles of the κB site and SBE in basal and stimulated expression of MUC1. Treatment of T47D cells and HMEC with interferon-γ (IFNγ) alone enhanced MUC1 expression at the level of transcription, and the effect of IFNγ was further stimulated by tumor necrosis factor-α (TNFα). MUC1 responsiveness to these cytokines was modest in T47D cells but clearly evident in HMEC. Transient transfection of T47D cells with mutant MUC1 promoter constructs revealed that the κB site at −589/−580 and the SBE at −503/−495 and were required for cooperative stimulation by TNFα and IFNγ. Electrophoretic mobility shift assays (EMSA) revealed that the synergy was mediated not by cooperative binding of transcription factors but by the independent actions of STAT1α and NFκB p65 on their respective binding sites. Independent mutations in the κB site and SBE abrogated cytokine responsiveness and reduced basal MUC1 promoter activity by 45–50%. However, only the κB site appeared to be constitutively activated in T47D cells, in part by NFκB p65. These findings implicate two cytokine response elements in the 5 ′ flanking region of MUC1, specifically a κB site and a STAT binding element, in overexpression of MUC1 in breast cancer cells. ^

Non B-DNA structures and genomic instability associated with myotonic dystrophy type 2 (CCTG).(CAGG) repeats

There are many diseases associated with the expansion of DNA repeats in humans. Myotonic dystrophy type 2 is one of such diseases, characterized by expansions of a (CCTG)•(CAGG) repeat tract in intron 1 of zinc finger protein 9 (ZNF9) in chromosome 3q21.3. The DM2 repeat tract contains a flanking region 5' to the tract that consists of a polymorphic repetitive sequence (TG)14-25(TCTG)4-11(CCTG) n. The (CCTG)•(CAGG) repeat is typically 11-26 repeats in persons without the disease, but can expand up to 11,000 repeats in affected individuals, which is the largest expansion seen in DNA repeat diseases to date. This DNA tract remains one of the least characterized disease-associated DNA repeats, and mechanisms causing the repeat expansion in humans have yet to be elucidated. Alternative, non B-DNA structures formed by the expanded repeats are typical in DNA repeat expansion diseases. These sequences may promote instability of the repeat tracts. I determined that slipped strand structure formation occurs for (CCTG)•(CAGG) repeats at a length of 42 or more. In addition, Z-DNA structure forms in the flanking human sequence adjacent to the (CCTG)•(CAGG) repeat tract. I have also performed genetic assays in E. coli cells and results indicate that the (CCTG)•(CAGG) repeats are more similar to the highly unstable (CTG)•(CAG) repeat tracts seen in Huntington's disease and myotonic dystrophy type 1, than to those of the more stable (ATTCT)•(AGAAT) repeat tracts of spinocerebellar ataxia type 10. This instability, however, is RecA-independent in the (CCTG)•(CAGG) and (ATTCT)•(AGAAT) repeats, whereas the instability is RecA-dependent in the (CTG)•(CAG) repeats. Structural studies of the (CCTG)•(CAGG) repeat tract and the flanking sequence, as well as genetic selection assays may reveal the mechanisms responsible for the repeat instability in E. coli, and this may lead to a better understanding of the mechanisms contributing to the human disease state. ^

A pair of adjacent glucocorticoid response elements regulate expression of two mouse metallothionein genes

Synthesis of mouse metallothionein (MT)-I and MT-II is transcriptionally induced by the synthetic glucocorticoid, dexamethasone (DEX) or both in vivo as well as in numerous cell lines. However, the location(s) of a glucocorticoid response element (GRE) has not been described. The observation that a marked MT-I gene, as well as heterologous genes, when placed in the context of 17 kb of flanking sequence from the MT locus, are inducible by DEX and lipopolysaccharide in transgenic mice renewed the search for the GRE. Analysis of a series of deletion constructs from this 17-kb region in cultured cells identified a single 455-bp region that conferred DEX induction on a reporter gene. This 455-bp region contains two GREs that bind to the glucocorticoid receptor as assessed by gel mobility shift. Deletion of this fragment from the 17-kb flanking region eliminates the DEX responsiveness of reporter genes. The two GREs, which are located ≈1 kb upstream of the MT-II gene and ≈7 kb upstream of the MT-I gene, are necessary for induction of both genes and can function independently of elements within the proximal promoter region of either gene.

Genomic structure and ecdysone regulation of the prophenoloxidase 1 gene in the malaria vector Anopheles gambiae

Prophenoloxidase, a melanin-synthesizing enzyme, is considered to be an important arthropod immune protein. In mosquitoes, prophenoloxidase has been shown to be involved in refractory mechanisms against malaria parasites. In our study we used Anopheles gambiae, the most important human malaria vector, to characterize the first arthropod prophenoloxidase gene at the genomic level. The complete nucleotide sequence, including the immediate 5′ flanking sequence (−855 bp) of the prophenoloxidase 1 gene, was determined. The gene spans 10 kb and is composed of five exons and four introns coding for a 2.5-kb mRNA. In the 5′ flanking sequence, we found several putative regulatory motifs, two of which were identified as ecdysteroid regulatory elements. Electrophoretic mobility gel-shift assays and supershift assays demonstrated that the Aedes aegypti ecdysone receptor/Ultraspiracle nuclear receptor complex, and, seemingly, the endogenous Anopheles gambiae nuclear receptor complex, was able to bind one of the ecdysteroid response elements. Furthermore, 20-hydroxyecdysone stimulation was shown to up-regulate the transcription of the prophenoloxidase 1 gene in an A. gambiae cell line.

An initiator element mediates autologous downregulation of the human type A γ-aminobutyric acid receptor β1 subunit gene

The regulated expression of type A γ-aminobutyric acid receptor (GABAAR) subunit genes is postulated to play a role in neuronal maturation, synaptogenesis, and predisposition to neurological disease. Increases in GABA levels and changes in GABAAR subunit gene expression, including decreased β1 mRNA levels, have been observed in animal models of epilepsy. Persistent exposure to GABA down-regulates GABAAR number in primary cultures of neocortical neurons, but the regulatory mechanisms remain unknown. Here, we report the identification of a TATA-less minimal promoter of 296 bp for the human GABAAR β1 subunit gene that is neuron specific and autologously down-regulated by GABA. β1 promoter activity, mRNA levels, and subunit protein are decreased by persistent GABAAR activation. The core promoter, 270 bp, contains an initiator element (Inr) at the major transcriptional start site. Three concatenated copies of the 10-bp Inr and its immediate 3′ flanking sequence produce full neural specific activity that is down-regulated by GABA in transiently transfected neocortical neurons. Taking these results together with those of DNase I footprinting, electrophoretic mobility shift analysis, and 2-bp mutagenesis, we conclude that GABA-induced down-regulation of β1 subunit mRNAs involves the differential binding of a sequence-specific basal transcription factor(s) to the Inr. The results support a transcriptional mechanism for the down-regulation of β1 subunit GABAAR gene expression and raises the possibility that altered levels of sequence-specific basal transcription factors may contribute to neurological disorders such as epilepsy.

Aldehyde dehydrogenase class 3 expression: Identification of a cornea-preferred gene promoter in transgenic mice

Aldehyde dehydrogenase class 3 (ALDH3) constitutes 20–40% of the total water-soluble proteins in the mammalian cornea. Here, we show by Northern blot analysis that ALDH3 expression in the mouse is at least 500-fold higher in the cornea than in any other tissue examined, with very low levels of expression detected in the stomach, urinary bladder, ocular lens, and lung. Histochemical localization reveals that this exceptional level of expression in the mouse cornea occurs in the anterior epithelial cells and that little ALDH3 is present in the keratocytes or corneal endothelial cells. A 13-kbp mouse ALDH3 promoter fragment containing >12 kbp of the 5′ flanking sequence, the 40-bp untranslated first exon, and 29 bp of intron 1 directed cat reporter gene expression to tissues that express the endogenous ALDH3 gene, except that transgene promoter activity was higher in the stomach and bladder than in the cornea. By contrast, when driven by a 4.4-kbp mouse ALDH3 promoter fragment [1,050-bp 5′ flanking region, exon 1, intron 1 (3.4 kbp), and 7 bp of exon 2] expression of the cat reporter gene was confined to the corneal epithelial cells, except for very low levels in the liver, effectively reproducing the corneal expression pattern of the endogenous ALDH3 gene. These results indicate that tissue-specific expression of ALDH3 is determined by positive and negative elements in the 5′ flanking region of the gene and suggests putative silencers located in intron 1. We demonstrate regulatory sequences capable of directing cornea-specific gene expression, affording the opportunity for genetic engineering in this transparent tissue.

Retroposon analysis of major cetacean lineages: The monophyly of toothed whales and the paraphyly of river dolphins

SINE (short interspersed element) insertion analysis elucidates contentious aspects in the phylogeny of toothed whales and dolphins (Odontoceti), especially river dolphins. Here, we characterize 25 informative SINEs inserted into unique genomic loci during evolution of odontocetes to construct a cladogram, and determine a total of 2.8 kb per taxon of the flanking sequences of these SINE loci to estimate divergence times among lineages. We demonstrate that: (i) Odontocetes are monophyletic; (ii) Ganges River dolphins, beaked whales, and ocean dolphins diverged (in this order) after sperm whales; (iii) three other river dolphin taxa, namely the Amazon, La Plata, and Yangtze river dolphins, form a monophyletic group with Yangtze River dolphins being the most basal; and (iv) the rapid radiation of extant cetacean lineages occurred some 28–33 million years B.P., in strong accord with the fossil record. The combination of SINE and flanking sequence analysis suggests a topology and set of divergence times for odontocete relationships, offering alternative explanations for several long-standing problems in cetacean evolution.

Efficient expression of protein coding genes from the murine U1 small nuclear RNA promoters.

Few promoters are active at high levels in all cells. Of these, the majority encode structural RNAs transcribed by RNA polymerases I or III and are not accessible for the expression of proteins. An exception are the small nuclear RNAs (snRNAs) transcribed by RNA polymerase II. Although snRNA biosynthesis is unique and thought not to be compatible with synthesis of functional mRNA, we have tested these promoters for their ability to express functional mRNAs. We have used the murine U1a and U1b snRNA gene promoters to express the Escherichia coli lacZ gene and the human alpha-globin gene from either episomal or integrated templates by transfection, or infection into a variety of mammalian cell types. Equivalent expression of beta-galactosidase was obtained from < 250 nucleotides of 5'-flanking sequence containing the complete promoter of either U1 snRNA gene or from the 750-nt cytomegalovirus promoter and enhancer regions. The mRNA was accurately initiated at the U1 start site, efficiently spliced and polyadenylylated, and localized to polyribosomes. Recombinant adenovirus containing the U1b-lacZ chimeric gene transduced and expressed beta-galactosidase efficiently in human 293 cells and airway epithelial cells in culture. Viral vectors containing U1 snRNA promoters may be an attractive alternative to vectors containing viral promoters for persistent high-level expression of therapeutic genes or proteins.

Characterization of promoters and stable transfection by homologous and nonhomologous recombination in Plasmodium falciparum.

Genetic studies of the protozoan parasite Plasmodium falciparum have been severely limited by the inability to introduce or modify genes. In this paper we describe a system of stable transfection of P. falciparum using a Toxoplasma gondii dihydrofolate reductase-thymidylate synthase gene, modified to confer resistance to pyrimethamine, as a selectable marker. This gene was placed under the transcriptional control of the P. falciparum calmodulin gene flanking sequences. Transfected parasites generally maintained plasmids episomally while under selection; however, parasite clones containing integrated forms of the plasmid were obtained. Integration occurred by both homologous and nonhomologous recombination. In addition to the flanking sequence of the P. falciparum calmodulin gene, the 5' sequences of the P. falciparum and P. chabaudi dihydrofolate reductase-thymidylate synthase genes were also shown to be transcriptionally active in P. falciparum. The minimal 5' sequence that possessed significant transcriptional activity was determined for each gene and short sequences containing important transcriptional control elements were identified. These sequences will provide considerable flexibility in the future construction of plasmid vectors to be used for the expression of foreign genes or for the deletion or modification of P. falciparum genes of interest.

The molecular basis of Sanfilippo syndrome type B.

The Sanfilippo syndrome type B is a lysosomal storage disorder caused by deficiency of alpha-N-acetylglucosaminidase; it is characterized by profound mental deterioration in childhood and death in the second decade. For understanding the molecular genetics of the disease and for future development of DNA-based therapy, we have cloned the cDNA and gene encoding alpha-N-acetylglucosaminidase. Cloning started with purification of the bovine enzyme and use of a conserved oligonucleotide sequence to probe a human cDNA library. The cDNA sequence was found to encode a protein of 743 amino acids, with a 20- to 23-aa signal peptide immediately preceding the amino terminus of the tissue enzyme and with six potential N-glycosylation sites. The 8.5-kb gene (NAGLU), interrupted by 5 introns, was localized to the 5'-flanking sequence of a known gene, EDH17B, on chromosome 17q21. Five mutations were identified in cells of patients with Sanfilippo syndrome type B: 503del10, R297X, R626X, R643H, and R674H. The occurrence of a frameshift and a nonsense mutation in homozygous form confirms the identity of the NAGLU gene.

Specific interactions outside the proline-rich core of two classes of Src homology 3 ligands.

Two dodecapeptides belonging to distinct classes of Src homology 3 (SH3) ligands and selected from biased phage display libraries were used to investigate interactions between a specificity pocket in the Src SH3 domain and ligant residues flanking the proline-rich core. The solution structures of c-Src SH3 complexed with these peptides were solved by NMR. In addition to proline-rich, polyproline type II helix-forming core, the class I and II ligands each possesses a flanking sequence that occupies a large pocket between the RT and n-Src loops of the SH3 domain. Structural and mutational analyses illustrate how the two classes of SH3 ligands exploit a specificity pocket on the receptor differently to increase binding affinity and specificity.