Evolutionary constraint in flanking regions of avian genes

An important comprehension from comparative genomic analysis is that sequence conservation beyond neutral expectations is frequently found outside protein-coding regions, indicating important functional roles of noncoding DNA. Understanding the causes of constraint on noncoding sequence evolution forms an important area of research, not least in light of the importance for understanding the evolution of gene expression. We aligned all orthologous genes of chicken and zebra finch together with 5 kb of their upstream and downstream noncoding sequences, to study the evolution of gene flanking sequences in the avian genome. Using ancestral repeats as a neutral reference, we detected significant evolutionary constraint in the 3' flanking region, highest directly after termination (60%) and then gradually decreasing to about 20% 5 kb downstream. Constraint was higher in annotated 3' untranslated regions (UTRs) than in non-UTRs at the same distance from the stop codon and higher in sequences annotated as microRNA (miRNA)-binding sites than in non-miRNA-binding sites within 3' UTRs. Constraint was also higher when estimated for a smaller data set of genes from more closely related songbird species, indicating turnover of functional elements during avian evolution. On the 5' flanking side constraint was readily seen within the first 125 bp immediately upstream of the start codon (34%) and was about 10% for remaining sequence within 5 kb upstream. Analysis of chicken polymorphism data gave further support for the highest constraint directly before and after the translated region. Finally, we found that genes evolving under the highest constraint measured by d(N)/d(S) also had the highest level of constraint in the 3' flanking region. This study broadens the insights into gene flanking sequence evolution by adding new findings from a vertebrate lineage other than mammals.     

Cloning of the gene for the thyrotropin beta subunit in the Japanese crested ibis,Nipponia nippon

We isolated a putative gene for the thyrotropin beta subunit (TSHbeta) from two types of genomic libraries of the Japanese crested ibis, Nipponia nippon. Exon-intron structure was deduced by comparing the determined sequence with those of TSH beta cDNA of other birds. The deduced amino acid sequence shows extensive similarities to those of the other birds, which assures our assumption that the acquired nucleotide sequence represents the TSHbeta gene. The assembled genomic fragment is 4192 bp in size and consists of 1937 bp of putative 5' flanking region followed by exon-intron structure with three exons and two introns, similar to those observed in rat, human and goldfish counterparts. Locations of introns are also similar to those in mammals and goldfish. Comparison of the 5' flanking region of the ibis TSHbeta gene with those of mammals reveals that several regulatory sequences, such as negative thyroid hormone responsive element (nTRE), Pit-1 responsive element, and AP-1 responsive element, which were characterized in mammalian TSHbeta genes, are also found in the promoter region. This is the first report on the exon-intron structure and 5' flanking region of the TSHbeta gene in an avian species.     

The structure of the human apolipoprotein C-II gene. Electron microscopic analysis of RNA:DNA hybrids, complete nucleotide sequence, and identification of 5' homologous sequences among apolipoprotein genes

Cloned human apo-C-II cDNA was used as a hybridization probe to identify the human apo-C-II gene in a genomic library constructed in our laboratory. The isolated apo-C-II DNA was studied both by electron microscopy and by direct sequence analysis. Ultrastructural morphological analysis of RNA-DNA hybrids revealed that the apo-C-II gene had complex structures because of regions of inverted complementary sequences in and around the gene forming stem-and-loop structures which interfere with the formation of stable RNA:DNA hybrids. Extensive morphological analysis revealed a minimum of 3 intervening sequences (IVS), and their lengths were measured. Direct sequence analysis of the cloned gene confirmed the presence of 3 IVS. There are 4 Alu type sequences in IVS-I. We sequenced 4340 nucleotides which include 545 nucleotides in the 5' flanking region, the entire gene which spans 3320 nucleotides, and 475 nucleotides in the 3' flanking region which also encompasses an additional Alu sequence. The 5' end of the gene was identified by primer extension and sequencing of the primer extended cDNA. Apo-C-II mRNA structure was deduced from the cDNA sequence, the primer extension experiments, and the genomic sequence. It is 494 nucleotides in length. Its sequence differs from previously published sequences in that there are 7 additional nucleotides before the polyadenylate tail. In the 5' flanking region, nucleotides -234 to -213 encompass a GC-rich region which exhibits high homology (greater than 70%) to the 5' flanking regions of the genes of all the apolipoproteins published to date, namely, apo-A-II (-497 to -471), apo-A-I (approximately -196 to -179), apo-E (-409 to -391), and apo-C-III (approximately -116 to -103). This highly conserved region might represent some evolutionarily conserved sequences from these related genes and/or might represent a region with regulatory function.     

Molecular structure and flanking nucleotide sequences of the natural chicken ovomucoid gene

Five independent clones containing the natural chicken ovomucoid gene have been isolated from a chicken gene library. One of these clones, CL21, contains the complete ovomucoid gene and includes more than 3 kb of DNA sequences flanking both termini of the gene. Restriction endonuclease mapping, electron microscopy and direct DNA sequencing analyses of this clone have revealed that the ovomucoid gene is 5.6 kb long and codes for a messenger RNA of 821 nucleotides. The structural gene sequence coding Ifor the mature messenger RNA is split into at least eight segments by a minimum of seven intervening sequences of various sizes. The shortest structural gene segment is only 20 nucleotides long. All seven intervening sequences are located within the peptide coding region of the gene, and the sequences at the 5' and 3' untranslated regions of the mRNA are not interrupted by intervening sequences. The DNA sequences of the regions flanking the 5' and 3' termini of the gene have been determined. Thirty nucleotides before the start of the messenger RNA coding sequence is the heptanucleotide TATATAT, which is also present in a similar location relative to the chicken ovalbumin gene and other unique sequence eucaryotic genes. This sequence resembles that of the Pribnow box in procaryotic genes where a promoter function has been implicated. Seven nucleotides past the 3' end of the gene is the tetranucleotide TTGT, a sequence found to be present at identical locations as either TTTT or TTGT in other eucaryotic genes that have been sequenced. These conserved DNA sequences flanking eucaryotic genes may serve some regulator function in the expression of these genes.     

The human alpha-fetoprotein gene. Sequence organization and the 5' flanking region

The human alpha-fetoprotein (AFP) gene was isolated into three overlapping clones in bacteriophage lambda vectors and its sequence organization analyzed by restriction endonuclease mapping and nucleotide sequencing. The human AFP gene is about 20 kilobase pairs long and contains 15 exons and 14 introns. The overall organization of the human AFP gene is similar to that of the mouse AFP gene, with all but two exons showing identical sizes. Nucleotide sequences at all exon/intron junctions display similarity to the consensus boundary sequence (Breathnach, R., and Chambon, P. (1981) Annu. Rev. Biochem. 50, 349-383), with the GT-AG rule applied to the splicing point. The cap site maps 44 nucleotides upstream from the translation initiation site. The "TATA box" is located 27 nucleotides upstream from the putative cap site and is flanked by sequences with dyad symmetry. The TATA box can thus be placed in the loop portion of a possible stem-loop structure formed by intrastrand base-pairing. Other characteristic nucleotide sequences in the 5' flanking region include a CCAAC pentamer, a 14-base pair (bp) enhancer-like sequence, and a 9-bp sequence homologous to the glucocorticoid responsive element. A long (90 bp) direct repeat and several alternating purine/pyrimidine sequences are also present in the 5' flanking region. A 736-bp sequence of the 5' flanking region adjacent to the cap site of the human AFP gene shows a 61% similarity with the corresponding region of the mouse AFP gene. There are two Alu family sequences and two poly(dT-dG) repeats in the human AFP gene that show different distribution patterns from those in the mouse AFP gene.     

Multiple mRNAs are generated from the chicken lysozyme gene

We have determined the DNA sequence of a 770 by Pst I fragment containing 450 nucleotides of the 5′ flanking region of the chicken lysozyme gene. S1-nuclease mapping was performed to localize the 5′ end of nuclear RNA containing lysozyme-specific sequences and of the mRNA. We present evidence that the 5′ noncoding region of the chicken lysozyme mRNA is heterogeneous in length. The 5′ termini of the different mRNAs map 29, 31 and 53 nucleotides upstream from their common initiation codon. The 5′ ends of lysozyme-specific nuclear RNAs map at positions similar to that of the mRNA. AT-rich regions and sequences similar to the E. coli RNA polymerase recognition sequence are found around 30 and 70 nucleotides upstream from each of these 5′ termini. The AT-rich regions differ, however, from the canonical Goldberg-Hogness box in that they do not contain the extremely conserved TATA sequence motif. Sequence comparison at the 5′ end of the lysozyme, conalbumin and ovalbumin genes reveals only one region of partial homology, 140 nucleotides upstream from the mRNA start sites.     

Complete nucleotide sequence of the chicken alpha A-crystallin gene and its 5' flanking region

The chicken alpha A-crystallin gene and 2.6 kb of its 5' flanking sequence have been isolated and characterized by electron microscopy and sequencing. The structural gene is 4.5 kb long and contains two introns, each approx. 1 kb in length. The first intron divides codons 63 and 64, and the second intron divides codons 104 and 105, as in rodents. There is little indication that the insert exon of rodents (an alternatively spliced sequence) is present in complete form in the chicken alpha A-crystallin gene; small stretches of similarity to this sequence were found throughout the gene. The 5' flanking sequence of the chicken alpha A-crystallin gene shows considerable sequence similarity with other mammalian alpha B-crystallin genes. In addition, one consensus sequence (GCAGCATGCCCTCCTAG) present in the 5' flanking region of the chicken alpha A-crystallin gene was found in the 5' flanking region of most reported crystallin genes.