Identification of a third distinct estrogen receptor and reclassification of estrogen receptors in teleosts

This paper describes three distinct estrogen receptor (ER) subtypes: ERα, ERβ, and a unique type, ERγ, cloned from a teleost fish, the Atlantic croaker Micropogonias undulatus; the first identification of a third type of classical ER in vertebrate species. Phylogenetic analysis shows that ERγ arose through gene duplication from ERβ early in the teleost lineage and indicates that ERγ is present in other teleosts, although it has not been recognized as such. The Atlantic croaker ERγ shows amino acid differences in regions important for ligand binding and receptor activation that are conserved in all other ERγs. The three ER subtypes are genetically distinct and have different distribution patterns in Atlantic croaker tissues. In addition, ERβ and ERγ fusion proteins can each bind estradiol-17β with high affinity. The presence of three functional ERs in one species expands the role of ER multiplicity in estrogen signaling systems and provides a unique opportunity to investigate the dynamics and mechanisms of ER evolution.

The essential protein encoded by the UL31 gene of herpes simplex virus 1 depends for its stability on the presence of UL34 protein

To pursue an earlier observation that the protein encoded by the UL34 gene binds to intermediate chain of dynein, we constructed a series of mutants from which sequences encoding the entire protein (ΔUL34) or amino-terminal [UL34Δ(3–119)] or carboxyl-terminal [UL34Δ(245–275)] domains were deleted. The mutant lacking the sequence encoding the carboxyl-terminal domain grew in all cell lines tested. The two other mutants replicated only in cell type-dependent manner and poorly. Rescue of ΔUL34 mutant with a fragment that does not encompass the UL31 ORF restored wild-type phenotype. UL34 protein interacts physically with UL31, and the UL31 deletion mutant appears to have a phenotype similar to that of UL34 deletion mutant. Experiments designed to determine whether the phenotypes of the deletion mutants have a common base revealed that cells infected with the ΔUL34 mutant accumulate UL31 RNA but not the corresponding protein. The UL31 protein accumulated, however, to near wild-type virus-infected cell levels in cells infected with ΔUL34 mutant and treated with the MG132 proteosomal inhibitor at 6 h after infection. This is evidence that the stability of an essential viral protein requires the presence of another protein. The observation raises the bar for identification of gene function on the basis of analyses of the phenotype of mutants in which the gene has been deleted or rendered inoperative.

Regional and strain-specific gene expression mapping in the adult mouse brain

To determine the genetic causes and molecular mechanisms responsible for neurobehavioral differences in mice, we used highly parallel gene expression profiling to detect genes that are differentially expressed between the 129SvEv and C57BL/6 mouse strains at baseline and in response to seizure. In addition, we identified genes that are differentially expressed in specific brain regions. We found that approximately 1% of expressed genes are differentially expressed between strains in at least one region of the brain and that the gene expression response to seizure is significantly different between the two inbred strains. The results lead to the identification of differences in gene expression that may account for distinct phenotypes in inbred strains and the unique functions of specific brain regions.

Identification of a renal-specific oxido-reductase in newborn diabetic mice

Aldose reductase (ALR2), a NADPH-dependent aldo-keto reductase (AKR), is widely distributed in mammalian tissues and has been implicated in complications of diabetes, including diabetic nephropathy. To identify a renal-specific reductase belonging to the AKR family, representational difference analyses of cDNA from diabetic mouse kidney were performed. A full-length cDNA with an ORF of 855 nt and yielding a ≈1.5-kb mRNA transcript was isolated from a mouse kidney library. Human and rat homologues also were isolated, and they had ≈91% and ≈97% amino acid identity with mouse protein. In vitro translation of the cDNA yielded a protein product of ≈33 kDa. Northern and Western blot analyses, using the cDNA and antirecombinant protein antibody, revealed its expression exclusively confined to the kidney. Like ALR2, the expression was up-regulated in diabetic kidneys. Its mRNA and protein expression was restricted to renal proximal tubules. The gene neither codistributed with Tamm–Horsfall protein nor aquaporin-2. The deduced protein sequence revealed an AKR-3 motif located near the N terminus, unlike the other AKR family members where it is confined to the C terminus. Fluorescence quenching and reactive blue agarose chromatography studies revealed that it binds to NADPH with high affinity (KdNADPH = 66.9 ± 2.3 nM). This binding domain is a tetrapeptide (Met-Ala-Lys-Ser) located within the AKR-3 motif that is similar to the other AKR members. The identified protein is designated as RSOR because it is renal-specific with properties of an oxido-reductase, and like ALR2 it may be relevant in the renal complications of diabetes mellitus.

Building a dictionary for genomes: Identification of presumptive regulatory sites by statistical analysis

The availability of complete genome sequences and mRNA expression data for all genes creates new opportunities and challenges for identifying DNA sequence motifs that control gene expression. An algorithm, “MobyDick,” is presented that decomposes a set of DNA sequences into the most probable dictionary of motifs or words. This method is applicable to any set of DNA sequences: for example, all upstream regions in a genome or all genes expressed under certain conditions. Identification of words is based on a probabilistic segmentation model in which the significance of longer words is deduced from the frequency of shorter ones of various lengths, eliminating the need for a separate set of reference data to define probabilities. We have built a dictionary with 1,200 words for the 6,000 upstream regulatory regions in the yeast genome; the 500 most significant words (some with as few as 10 copies in all of the upstream regions) match 114 of 443 experimentally determined sites (a significance level of 18 standard deviations). When analyzing all of the genes up-regulated during sporulation as a group, we find many motifs in addition to the few previously identified by analyzing the subclusters individually to the expression subclusters. Applying MobyDick to the genes derepressed when the general repressor Tup1 is deleted, we find known as well as putative binding sites for its regulatory partners.

Human sulfite oxidase R160Q: Identification of the mutation in a sulfite oxidase-deficient patient and expression and characterization of the mutant enzyme

Sulfite oxidase catalyzes the terminal reaction in the degradation of sulfur amino acids. Genetic deficiency of sulfite oxidase results in neurological abnormalities and often leads to death at an early age. The mutation in the sulfite oxidase gene responsible for sulfite oxidase deficiency in a 5-year-old girl was identified by sequence analysis of cDNA obtained from fibroblast mRNA to be a guanine to adenine transition at nucleotide 479 resulting in the amino acid substitution of Arg-160 to Gln. Recombinant protein containing the R160Q mutation was expressed in Escherichia coli, purified, and characterized. The mutant protein contained its full complement of molybdenum and heme, but exhibited 2% of native activity under standard assay conditions. Absorption spectroscopy of the isolated molybdenum domains of native sulfite oxidase and of the R160Q mutant showed significant differences in the 480- and 350-nm absorption bands, suggestive of altered geometry at the molybdenum center. Kinetic analysis of the R160Q protein showed an increase in Km for sulfite combined with a decrease in kcat resulting in a decrease of nearly 1,000-fold in the apparent second-order rate constant kcat/Km. Kinetic parameters for the in vitro generated R160K mutant were found to be intermediate in value between those of the native protein and the R160Q mutant. Native sulfite oxidase was rapidly inactivated by phenylglyoxal, yielding a modified protein with kinetic parameters mimicking those of the R160Q mutant. It is proposed that Arg-160 attracts the anionic substrate sulfite to the binding site near the molybdenum.

The 3′ untranslated region of a rice α-amylase gene functions as a sugar-dependent mRNA stability determinant

In plants, sugar feedback regulation provides a mechanism for control of carbohydrate allocation and utilization among tissues and organs. The sugar repression of α-amylase gene expression in rice provides an ideal model for studying the mechanism of sugar feedback regulation. We have shown previously that sugar repression of α-amylase gene expression in rice suspension cells involves control of both transcription rate and mRNA stability. The α-amylase mRNA is significantly more stable in sucrose-starved cells than in sucrose-provided cells. To elucidate the mechanism of sugar-dependent mRNA turnover, we have examined the effect of αAmy3 3′ untranslated region (UTR) on mRNA stability by functional analyses in transformed rice suspension cells. We found that the entire αAmy3 3′ UTR and two of its subdomains can independently mediate sugar-dependent repression of reporter mRNA accumulation. Analysis of reporter mRNA half-lives demonstrated that the entire αAmy3 3′ UTR and the two subdomains each functioned as a sugar-dependent destabilizing determinant in the turnover of mRNA. Nuclear run-on transcription analysis further confirmed that the αAmy3 3′ UTR and the two subdomains did not affect the transcription rate of promoter. The identification of sequence elements in the α-amylase mRNA that dictate the differential stability has very important implications for the study of sugar-dependent mRNA decay mechanisms.

Identification of genes differentially regulated by interferon α, β, or γ using oligonucleotide arrays

The pleiotropic activities of interferons (IFNs) are mediated primarily through the transcriptional regulation of many downstream effector genes. The mRNA profiles from IFN-α, -β, or -γ treatments of the human fibrosarcoma cell line, HT1080, were determined by using oligonucleotide arrays with probe sets corresponding to more than 6,800 human genes. Among these were transcripts for known IFN-stimulated genes (ISGs), the expression of which were consistent with previous studies in which the particular ISG was characterized as responsive to either Type I (α, β) or Type II (γ) IFNs, or both. Importantly, many novel IFN-stimulated genes were identified that were diverse in their known biological functions. For instance, several novel ISGs were identified that are implicated in apoptosis (including RAP46/Bag-1, phospholipid scramblase, and hypoxia inducible factor-1α). Furthermore, several IFN-repressed genes also were identified. These results demonstrate the usefulness of oligonucleotide arrays in monitoring mammalian gene expression on a broad and unprecedented scale. In particular, these findings provide insights into the basic mechanisms of IFN actions and ultimately may contribute to better therapeutic uses for IFNs.

Repression of transcription mediated by dual elements in the CCAAT/enhancer binding protein α gene

During adipocyte differentiation, the expression of C/EBPα is activated, which in turn serves to transcriptionally activate numerous adipocyte genes. A previous search for cis elements that regulate transcription of the C/EBPα gene led to the identification of a potential repressive element within the proximal 5′ flanking region of the gene. Nuclear extracts from 3T3-L1 preadipocytes, but not adipocytes, were found to contain a factor, CUP (C/EBPα undifferentiated protein), that binds to this site (the CUP-1 site). In the present investigation, we show that C/EBPα promoter-luciferase constructs containing both the proximal 5′ flanking and the entire 5′ untranslated regions of the gene exhibit an expression pattern during adipocyte differentiation comparable to that of the endogenous C/EBPα gene. Mutation of the CUP-1 site in these constructs had little effect on reporter gene expression; however, when this mutation was combined with deletion of the 5′ untranslated region, reporter gene expression by preadipocytes was dramatically up-regulated. Consistent with this finding, a second CUP binding site (the CUP-2 site) was identified in the 5′ untranslated region. Although mutation of either CUP element in constructs containing both the 5′ flanking and 5′ untranslated region had little effect on reporter gene transcription, mutation of both CUP elements markedly activated transcription. Thus, it appears that dual CUP regulatory elements repress transcription of the C/EBPα gene prior to induction of the adipocyte differentiation program.

Plastid-localized acetyl-CoA carboxylase of bread wheat is encoded by a single gene on each of the three ancestral chromosome sets

5′-End fragments of two genes encoding plastid-localized acetyl-CoA carboxylase (ACCase; EC 6.4.1.2) of wheat (Triticum aestivum) were cloned and sequenced. The sequences of the two genes, Acc-1,1 and Acc-1,2, are 89% identical. Their exon sequences are 98% identical. The amino acid sequence of the biotin carboxylase domain encoded by Acc-1,1 and Acc-1,2 is 93% identical with the maize plastid ACCase but only 80–84% identical with the cytosolic ACCases from other plants and from wheat. Four overlapping fragments of cDNA covering the entire coding region were cloned by PCR and sequenced. The wheat plastid ACCase ORF contains 2,311 amino acids with a predicted molecular mass of 255 kDa. A putative transit peptide is present at the N terminus. Comparison of the genomic and cDNA sequences revealed introns at conserved sites found in the genes of other plant multifunctional ACCases, including two introns absent from the wheat cytosolic ACCase genes. Transcription start sites of the plastid ACCase genes were estimated from the longest cDNA clones obtained by 5′-RACE (rapid amplification of cDNA ends). The untranslated leader sequence encoded by the Acc-1 genes is at least 130–170 nucleotides long and is interrupted by an intron. Southern analysis indicates the presence of only one copy of the gene in each ancestral chromosome set. The gene maps near the telomere on the short arm of chromosomes 2A, 2B, and 2D. Identification of three different cDNAs, two corresponding to genes Acc-1,1 and Acc-1,2, indicates that all three genes are transcriptionally active.

Identification of two novel transmembrane γ-carboxyglutamic acid proteins expressed broadly in fetal and adult tissues

The proline-rich γ-carboxyglutamic acid (Gla) proteins (PRGPs) 1 and 2 are the founding members of a family of vitamin K-dependent single-pass integral membrane proteins characterized by an extracellular amino terminal domain of approximately 45 amino acids that is rich in Gla. The intracellular carboxyl terminal region of these two proteins contains one or two copies of the sequence PPXY, a motif present in a variety of proteins involved in such diverse cellular functions as signal transduction, cell cycle progression, and protein turnover. In this report, we describe the cloning of the cDNAs for two additional human transmembrane Gla proteins (TMG) of 20–24 kDa named TMG3 and TMG4. These two proteins possess extracellular Gla domains with 13 or 9 potential Gla residues, respectively, followed by membrane-spanning hydrophobic regions and cytoplasmic carboxyl terminal regions that contain PPXY motifs. This emerging family of integral membrane Gla proteins includes proline-rich Gla protein (PRGP) 1, PRGP2, TMG3, and TMG4, all of which are characterized by broad and variable distribution in both fetal and adult tissues. Members of this family can be grouped into two subclasses on the basis of their gene organization and amino acid sequence. These observations suggest novel physiological functions for vitamin K beyond its known role in the biosynthesis of proteins involved in blood coagulation and bone development. The identification and characterization of these proteins may allow a more complete understanding of the teratogenic consequences of exposure in utero to vitamin K antagonists, such as warfarin-based anticoagulants.

Identification of two short internal ribosome entry sites selected from libraries of random oligonucleotides

Sequences that control translation of mRNA may play critical roles in regulating protein levels. One such element is the internal ribosome entry site (IRES). We previously showed that a 9-nt segment in the 5′ leader sequence of the mRNA encoding Gtx homeodomain protein could function as an IRES. To identify other short sequences with similar properties, we designed a selection procedure that uses a retroviral vector to express dicistronic mRNAs encoding enhanced green and cyan fluorescent proteins as the first and second cistrons, respectively. Expression of the second cistron was dependent upon the intercistronic sequences and was indicative of IRES activity. B104 cells were infected with two retroviral libraries that contained random sequences of 9 or 18 nt in the intercistronic region. Cells expressing both cistrons were sorted, and sequences recovered from selected cells were reassayed for IRES activity in a dual luciferase dicistronic mRNA. Two novel IRESes were identified by this procedure, and both contained segments with complementarity to 18S rRNA. When multiple copies of either segment were linked together, IRES activities were dramatically enhanced. Moreover, these synthetic IRESes were differentially active in various cell types. These properties are similar to those of the previously identified 9-nt IRES module from Gtx mRNA. These results provide further evidence that short nucleotide sequences can function as IRESes and support the idea that some cellular IRESes may be composed of shorter functional modules. The ability to identify IRES modules with specific expression properties may be useful in the design of vectors for biotechnology and gene therapy.

Identification of sample-specific sequences in mammalian cDNA and genomic DNA by the novel ligation-mediated subtraction (Limes)

The representational difference analysis (RDA) and other subtraction techniques are used to enrich sample-specific sequences by elimination of ubiquitous sequences existing in both the sample of interest (tester) and the subtraction partner (driver). While applying the RDA to genomic DNA of cutaneous lymphoma cells in order to identify tumor relevant alterations, we predominantly isolated repetitive sequences and artificial repeat-mediated fusion products of otherwise independent PCR fragments (PCR hybrids). Since these products severely interfered with the isolation of tester-specific fragments, we developed a considerably more robust and efficient approach, termed ligation-mediated subtraction (Limes). In first applications of Limes, genomic sequences and/or transcripts of genes involved in the regulation of transcription, such as transforming growth factor β stimulated clone 22 related gene (TSC-22R), cell death and cytokine production (caspase-1) or antigen presentation (HLA class II sequences), were found to be completely absent in a cutaneous lymphoma line. On the assumption that mutations in tumor-relevant genes can affect their transcription pattern, a protocol was developed and successfully applied that allows the identification of such sequences. Due to these results, Limes may substitute/supplement other subtraction/comparison techniques such as RDA or DNA microarray techniques in a variety of different research fields.

Identification of human DNA helicase V with the far upstream element-binding protein

The properties of human DNA helicase V (HDH V) were studied in greater detail following an improved purification procedure. From 450 g of cultured cells, <0.1 mg of pure protein was isolated. HDH V unwinds DNA unidirectionally by moving in the 3′ to 5′ direction along the bound strand in an ATP- and Mg2+-dependent fashion. The enzyme is not processive and can also unwind partial RNA–RNA duplexes such as HDH IV and HDH VIII. The Mr determined by SDS–PAGE (66 kDa) corresponds to that measured under native conditions, suggesting that HDH V exists as a monomer in the nucleus. Microsequencing of the purified HDH V shows that this enzyme is identical to the far upstream element-binding protein (FBP), a protein that stimulates the activity of the c-myc gene by binding specifically to the ‘FUSE’ DNA region localized upstream of its promoter. The sequence of HDH V/FBP contains RGG motifs like HDH IV/nucleolin, HDH VIII/G3BP as well as other human RNA and DNA helicases identified by other laboratories.

The TIGR Gene Indices: analysis of gene transcript sequences in highly sampled eukaryotic species

While genome sequencing projects are advancing rapidly, EST sequencing and analysis remains a primary research tool for the identification and categorization of gene sequences in a wide variety of species and an important resource for annotation of genomic sequence. The TIGR Gene Indices (http://www.tigr.org/tdb/tgi.shtml) are a collection of species-specific databases that use a highly refined protocol to analyze EST sequences in an attempt to identify the genes represented by that data and to provide additional information regarding those genes. Gene Indices are constructed by first clustering, then assembling EST and annotated gene sequences from GenBank for the targeted species. This process produces a set of unique, high-fidelity virtual transcripts, or Tentative Consensus (TC) sequences. The TC sequences can be used to provide putative genes with functional annotation, to link the transcripts to mapping and genomic sequence data, to provide links between orthologous and paralogous genes and as a resource for comparative sequence analysis.