Characterization of three splice variants and genomic organization of the mouse BMAL1 gene.

The BMAL1 gene encodes a member of the basic helix-loop-helix/PER-ARNT-SIM (bHLH/PAS) family of transcription factors. It is a key regulator of circadian rhythms. Using sequence information from human BMAL1 (hBMAL1) cDNAs previously reported by our laboratory, we have isolated and characterized cDNAs encoding three splice variants of the mouse BMAL1 (mBMAL1) gene. Of the three splice variants, mBMAL1b extends for 1878 bp in the coding sequence, which is 91% identical to that of hBMAL1b; its deduced amino acid sequence is 626 residues long and is 98% identical to that of hBMAL1b, and sequence identities in the bHLH, PAS-A, and PAS-B regions are 98, 100, and 100%, respectively. mBMAL1b' arises from alternative usage of exon 2, which results in a 7-amino-acid insertion and alternative splice acceptor usage at the intron 9/exon 10 splice junction, which causes an alanine residue deletion. mBMAL1b' encodes 632 amino acids and contains the bHLH/PAS domains. mBMAL1g' is generated by alternative splice acceptor usage at the intron 6/exon 7 splice junction, which results in a 28-bp deletion adjacent to the 5' end of the PAS domain. Since the 28-bp deletion shifts the reading frame, mBMAL1g' is predicted to encode a product of only 222 amino acids that lacks the PAS domain. The tissue distributions of the three splice variants showed some variation. The variations in the tissue distributions and predicted amino acid sequences suggest that the three splice variants may have different functions. Direct sequencing of the genomic mBMAL1 clones indicated that the coding sequence of mBMAL1 spans 32 kb and includes 17 exons. An unusual exon/intron donor sequence was found in intron 14, which begins with GC at the 5' end. Comparison with the bHLH/PAS family genes revealed that the intron/exon splice pattern of mBMAL1 most closely matches that of the mAhr, which suggests that BMAL1 and Ahr belong to the same subclass and may be derived from a common primordial gene.

Yu W ，Ikeda M ，Abe H ，Honma S ，Ebisawa T ，Yamauchi T ，Honma K ，Nomura M ... -  《BIOCHEMICAL AND BIOPHYSICAL RESEARCH COMMUNICATIONS》

cDNA cloning and tissue-specific expression of a novel basic helix-loop-helix/PAS protein (BMAL1) and identification of alternatively spliced variants with alternative translation initiation site usage.

Ikeda M ，Nomura M 《BIOCHEMICAL AND BIOPHYSICAL RESEARCH COMMUNICATIONS》

Molecular cloning and characterization of the human CLOCK gene: expression in the suprachiasmatic nuclei.

The Clock gene is an essential regulator of circadian rhythms. It encodes a member of the basic helix-loop-helix/PER-ARNT-SIM family of transcription factors known to play a central role in the control of diverse cellular events. Previously we described the functional identification and molecular isolation of the Clock gene in the mouse, its interaction with the BMAL1 protein, and the role of this complex as a transcriptional activator in the circadian pacemaker. Here, we report the cloning, exon organization, chromosomal location, and mRNA expression of the human CLOCK gene. The coding sequence of human CLOCK extends for 2538 bp and is 89% identical to its mouse ortholog; its deduced amino acid sequence is 846 residues long and is 96% identical to mouse CLOCK. Radiation hybrid mapping localized human CLOCK to the long arm of human chromosome 4 (4q12). Direct sequencing of a genomic CLOCK clone indicated that the coding sequence of human CLOCK extends over 20 exons and that its intron/exon organization is identical to that of the mouse ortholog. Northern blot analysis indicated widespread expression of two major transcripts of 8 and 10 kb, and in situ hybridization of human brain tissue revealed elevated expression of CLOCK mRNA in the suprachiasmatic nuclei, the locus of circadian control in mammals, and in the cerebellum. Comparison of cDNA clones revealed two single nucleotide polymorphisms in noncoding sequence flanking the CLOCK open reading frame. The central role of Clock in the organization of circadian rhythms suggests that it will be a useful candidate gene for genetic analyses of disorders associated with dysfunction of the circadian system.

Steeves TD ，King DP ，Zhao Y ，Sangoram AM ，Du F ，Bowcock AM ，Moore RY ，Takahashi JS ... -  《GENOMICS》

[Analysis, identification and correction of some errors of model refseqs appeared in NCBI Human Gene Database by in silico cloning and experimental verification of novel human genes].

We found that human genome coding regions annotated by computers have different kinds of many errors in public domain through homologous BLAST of our cloned genes in non-redundant (nr) database, including insertions, deletions or mutations of one base pair or a segment in sequences at the cDNA level, or different permutation and combination of these errors. Basically, we use the three means for validating and identifying some errors of the model genes appeared in NCBI GENOME ANNOTATION PROJECT REFSEQS: (I) Evaluating the support degree of human EST clustering and draft human genome BLAST. (2) Preparation of chromosomal mapping of our verified genes and analysis of genomic organization of the genes. All of the exon/intron boundaries should be consistent with the GT/AG rule, and consensuses surrounding the splice boundaries should be found as well. (3) Experimental verification by RT-PCR of the in silico cloning genes and further by cDNA sequencing. And then we use the three means as reference: (1) Web searching or in silico cloning of the genes of different species, especially mouse and rat homologous genes, and thus judging the gene existence by ontology. (2) By using the released genes in public domain as standard, which should be highly homologous to our verified genes, especially the released human genes appeared in NCBI GENOME ANNOTATION PROJECT REFSEQS, we try to clone each a highly homologous complete gene similar to the released genes in public domain according to the strategy we developed in this paper. If we can not get it, our verified gene may be correct and the released gene in public domain may be wrong. (3) To find more evidence, we verified our cloned genes by RT-PCR or hybrid technique. Here we list some errors we found from NCBI GENOME ANNOTATION PROJECT REFSEQs: (1) Insert a base in the ORF by mistake which causes the frame shift of the coding amino acid. In detail, abase in the ORF of a gene is a redundant insertion, which causes a reading frame shift in the translation of an alternative protein, such as LOC124919 is wrong form of C17 orf32 (with mouse and rat orthologs determined by us). (2) Put together by mistake (with force). This is a wrong assembly of non-relating cDNA segment, such as LOC147007 is wrong form of C17orf32. (3) Mistakenly insert a base or one section of cDNA in the ORF which causes it ending beforehand, only coding cDNA sequence of N-terminal amino acids, incomplete. For example, LOC123722 is wrong form of SPRYD1, and even the human hypothetical gene LOC126250 or PDCD5 is wrong form of our PDCD5 (TFAR19). (4) Incomplete, only coding cDNA sequence of C-terminal amino acids. For example, human LOC149076 and mouse LOC230761 are wrong form of our verified human ZNF362 and mouse Zfp362, respectively. (5) Incomplete, only coding one section of coding protein cDNA sequence of correct gene ORF, lacking N-terminal and C-terminal amino acids sequence, and at the same time, mistakenly anticipates the first non-initiation codon amino acid of the incomplete protein amino acid as the initiation codon, e.g. anticipating L as M. For example, LOC200084 is wrong form of ZNF362. (6) Mistakenly insert a base or one section of cDNA in the ORF, wrongly causing unwanted termination codon before the insertion, so the coding protein lacks the first part of the amino acids. For example, the GenBank Acc. No. AL096883 ( LOCUS No. HS323M22B) is wrong form of an experimentally verified human NM_012263 with mouse ortholog of BC010510 determined. (7) It may regard the polluted genomic sequence as complete gene cDNA sequence and anticipate the so-called single exon gene, even the real one, only a small ORF in the very long single exon mRNA, while there really exists termination code in the same phase of the upper part of the ORF initiation code, no other characters accord with the gene's condition. For example, LOC91126 is wrong form of ZNF362. (8) The anticipated genes only have ORF which has no EST proofs on both terminal sides. Depending on this ORF, a complete gene cDNA with double support of EST and human genome (there are termination codes at the same phase of the upper part of ORF) which indicates the anticipated ORF reference sequence may be incorrect. For example, LOC164395 may be wrong form of novel human gene bankit4590055. (9) A similar but smaller protein-coding gene is anticipated in the range of the human genome sequence that has the support of EST experimental proof, so other new anticipated gene may be incorrect. For example, LOC167563 may be wrong form of CMYA5. However,these errors can be corrected or avoided by using our strategy. Here we give one example in detail: Comparision of the sequence SPRYD1 with human hypothetical gene LOC123722. The TAA bases in the position of 478-480 in LOC123722 cDNA is redundant, which causes a reading frame shift in the translation of an alternative protein. The redundancy of GTAAA of LOC123722 is not supported by our experimental clone,and is almost fully rejected by human EST alignment, and is shown as the next intron sequence by genomic GT/AG organization analysis. The verification of cDNA or genomic DNA sequence of SPRYD1 implies that LOC123722 has a wrong stop codon within its ORF because of the prediction program, thus being not complete cds. To sum up, by combining bioinformatics analyses with experimental verification, we have found that there are many errors of at least nine kinds appeared in NCBI GENOME ANNOTATION PROJECT REFSEQs through BLAST of our cloned genes in non-redundant database, and our strategy is helpful in correcting them, such as LOC14907, LOC200084 and LOC91126 (all of them should be ZNF362, but are three different kinds of wrong forms of ZNF362), three model reference sequences predicted from NCBI contig NT_004511 by automated computational analysis using gene prediction method, or such as LOC124919 and LOC147007 (both should be C17orf32, but are two different kinds of wrong forms of C17orf32), two model reference sequences predicted from NCBI contig NT_010808 by automated computational analysis using gene prediction method. Therefore, the correct identification and annotation of novel human genes may be still a heavy task, which can be finished within a long period of time. So human genome coding regions annotated by computer should be used with caution. The articles published in the past did not clearly point out the existence of mistakes in the NCBI human gene mode reference sequence. At the Seventh International Human Genome Conference held in April 2002, we first published the researching result on this aspect in the communication form of Posterly insert a base or one section of cDNA in the ORF, wrongly causing unwanted termination codon before the insertion, so the coding protein lacks the first part of the amino acids. For example, the GenBank Acc. No. AL096883 ( LOCUS No. HS323M22B) is wrong form of an experimentally verified human NM_012263 with mouse ortholog of BC010510 determined. (7) It may regard the polluted genomic sequence as complete gene cDNA sequence and anticipate the so-called single exon gene, even the real one, only a small ORF in the very long single exon mRNA, while there really exists termination code in the same phase of the upper part of the ORF initiation code, no other characters accord with the gene's condition. For example, LOC91126 is wrong form of ZNF362. (8) The anticipated genes only have ORF which has no EST proofs on both terminal sides. Depending on this ORF, a complete gene cDNA with double support of EST and human genome (there are termination codes at the same phase of the upper part of ORF) which indicates the anticipated ORF reference sequence may be incorrect. For example, LOC164395 may be wrong form of novel human gene bankit4590055. (9) A similar but smaller protein-coding gene is anticipated in the range of the human genome sequence that has the support of EST experimental proof, so other new anticipated gene may be incorrect. For example, LOC167563 may be wrong form of CMYA5. However, these errors can be corrected or avoided by using our strategy. Here we give one example in detail: Comparision of the sequence SPRYD1 with human hypothetical gene LOC123722. The TAA bases in the position of 478-480 in LOC123722 cDNA is redundant, which causes a reading frame shift in the translation of an alternative protein. The redundancy of GTAAA of LOC123722 is not supported by our experimental clone, and is almost fully rejected by human EST alignment, and is shown as the next intron sequence by genomic GT/AG organization analysis. The verification of cDNA or genomic DNA sequence of SPRYD1 implies that LOC123722 has a wrong stop codon within its ORF because of the prediction program, thus being not complete cds. To sum up, by combining bioinformatics analyses with experimental verification, we have found that there are many errors of at least nine kinds appeared in NCBI GENOME ANNOTATION PROJECT REFSEQs through BLAST of our cloned genes in non-redundant database, and our strategy is helpful in correcting them, such as LOC14907, LOC200084 and LOC91126 (all of them should be ZNF362, but are three different kinds of wrong forms of ZNF362), three model reference sequences predicted from NCBI contig NT_004511 by automated computational analysis using gene prediction method, or such as LOC124919 and LOC147007 (both should be C17orf32, but are two different kinds of wrong forms of C17orf32), two model reference sequences predicted from NCBI contig NT_010808 by automated computational analysis using gene prediction method. Therefore, the correct identification and annotation of novel human genes may be still a heavy task, which can be finished within a long period of time. So human genome coding regions annotated by computer should be used with caution. (ABSTRACT TRUNCATED)

Zhang DL ，Ji L ，Li YD 《-》

Genomic organization and tissue-specific expression of splice variants of mouse organic anion transporting polypeptide 2.

cDNAs that code for mouse organic anion transporting polypeptide 2 (oatp2) have been cloned. At least three forms of mouse oatp2 cDNAs containing the same coding sequence were isolated. The common coding sequence is for a protein of 670 amino acids with 12 putative transmembrane domains. The deduced amino acid sequence of the mouse oatp2 shares 89% identity with the reported rat oatp2. Cloning and analysis of mouse oatp2 gene indicates that these isoforms are alternatively spliced products from the same gene. Heterogeneity was observed in the 5'-untranslated region of the cDNAs. Two of the three isoforms lacked the noncoding exon 3 sequence. Northern-blot hybridization analysis using the exon 3-specific probes demonstrated that mouse oatp2 mRNA containing exon 3 sequence is expressed in heart and lung, whereas exon 1-, 2-, and 17-specific probes detected mRNA only in brain and liver. The mouse oatp2 gene consists of 17 exons, including three noncoding exons, and 16 introns. All of the introns are flanked by GT-AG splice sequences except for intron 10 that is flanked by GC-AG splice sequence.

Ogura K ，Choudhuri S ，Klaassen CD 《BIOCHEMICAL AND BIOPHYSICAL RESEARCH COMMUNICATIONS》