Heterogeneous yet similar introns reside in identical positions of the rRNA genes in natural isolates of the archaeon Aeropyrum pernix

Some archaeal ribosomal DNA (rDNA) introns carry homing endonuclease-like genes and are therefore assumed to propagate by "intron homing". A previous study demonstrated that three introns are located within the rRNA operon (arnSL) of Aeropyrum pernix strain K1, two of which, Ialpha and Igamma, harbor open reading frames (ORFs) encoding putative LAGLIDADG-type endonucleases. In an effort to understand further the rDNA intron distribution in natural A. pernix populations, 11 A. pernix strains were isolated from marine hydrothermal biotopes, and comparative nucleotide sequence analysis of the arnSL alleles was performed. Of the 11 isolates, eight contained multiple introns, and three patterns of intron insertion were found. Three novel introns, Idelta (62 bp in length), Ivarepsilon (122 bp) and Izeta (57 bp) were identified. They were all ORF-less, but their predicted RNA secondary structure at the exon-intron junctions was consistent with the bulge-helix-bulge motif. The insertion positions and the terminal inverted repeat sequences of Idelta and Izeta were in agreement with those of Ialpha and Igamma, respectively. This suggests that these intron variants were generated by large indels (insertions/deletions) during their evolution.     

Equal G and C contents in histone genes indicate selection pressures on mRNA secondary structure

Summary Protein-specific versus taxon-specific patterns of nucleotide frequencies were studied in histone genes. The third positions of codons have a (well-known) taxon-specific G+C level and a histone type-specific G/C ratio. This ratio counterbalances the G/C ratio in the first and second positions so that the overall G and C levels in the coding region become approximately equal. The compensation of the G/C ratio indicates a selection pressure at the mRNA level rather than a selection pressure or mutation bias at the DNA level or a selection pressure on codon usage. The structure of histone mRNAs is compatible with the hypothesis that the G/C compensation is due to selection pressures on mRNA secondary structure. Nevertheless, no specific motifs seem to have been selected, and the free energy of the secondary structures is only slightly lower than that expected on the basis of nucleotide frequencies.Offprint requests to: M. A. Huynen     

Complementary addressed modification and cleavage of a single stranded DNA fragment with alkylating oligonucleotide derivatives.

A single stranded DNA fragment was modified with alkylating derivatives of oligonucleotides complementary to a certain nucleotide sequences in the fragment. The derivatives carried aromatic 2-chloroethylamino groups at their 3'- or 5'-terminal nucleotide residues. Some of the derivatives carried both alkylating group and intercalating phenazine group which stabilized complementary complexes. It was found that these oligonucleotide derivatives modify the DNA fragment in a specific way near the target complementary nucleotide sequences, and the DNA fragment can be cleaved at the alkylated nucleotides positions. Alkylating derivatives carrying phenazine groups were found to be the most efficient in reaction with the DNA fragment.     

Nuclear and mitochondrial genes contain similar phylogenetic signal for pigeons and doves (Aves: Columbiformes)

Molecular systematic studies generally assume that gene trees are reasonable estimates of species trees. We tested the validity of this assumption in the pigeons and doves (Aves: Columbiformes) by comparing phylogenies derived from nuclear (beta-fibrinogen intron 7) and mitochondrial (cytochrome b) genes. Trees derived from the two genes when analyzed separately contained many nodes in common. A partition homogeneity test revealed no significant incongruence between trees derived from the two genes; so, we combined nuclear and mitochondrial data in subsequent phylogenetic analyses. The resulting tree, which was highly resolved and generally well supported, contained a strong biogeographic component. The rate of nucleotide substitution for the nuclear intron was approximately six times slower than that of cytochrome b. This resulted in a much higher consistency index for trees derived from the intron because of the low level of multiple substitution. However, the degree of resolution and support for trees reconstructed from the two genes was similar. We also examined the transition and transversion substitution rates for the genes. Third position transversions for cytochrome b accumulated linearly with intron divergence, suggesting low levels of multiple substitution for third position transversions.     

Structure of a hemoglobin gene cluster and nucleotide sequence of three hemoglobin genes from the midge Chironomus thummi piger (Diptera, Insecta)

The aquatic larvae of the genus Chironomus (Diptera, Insecta) contain at least 12 different hemoglobin (Hb) variants in their hemolymph. In the present study we have analysed the structure and part of the nucleotide sequence of a Hb gene cluster cloned from the genomic DNA of Chironomus thummi piger. The cluster contains probably 6 different genes, separated by intergenic regions of various lengths. The nucleotide sequence of three putative Hb genes including the intergenic regions is presented. The inferred amino-acid sequences show clearly that two of these putative genes code for subvariants of the Hb variant VIIB. The third gene codes for a so far unknown Hb protein. As known already for other chironomid Hb genes, there are no intron sequences present in the coding regions.     

Nucleotide sequence of the cro-cII-oop region of bacteriophage 434 DNA.

The nucleotide sequence of a 869 bp segment of phage 434 DNA including the regulatory genes cro and cII is presented and compared with the corresponding part of the phage lambda DNA sequence. The 434 cro protein as deduced from the DNA sequence is a highly basic protein of 71 amino acid residues with a calculated molecular weight of 8089. While the cro gene sequences of phage 434 and lambda DNA are very different, the nuleotide sequences to the right of the lambda imm434 boundary show differences only at 11 out of 512 positions. Nucleotide substitutions in the cII gene occur with one exception in the third positions of the respective codons and only one out of several DNA regulatory signals located in this region of the phage genomes is affected by these nucleotide substitutions.     

Structure of the human prealbumin gene

Using cloned human prealbumin cDNA as a probe, Southern blot hybridization of human genomic DNA revealed that the prealbumin gene consists of an unique, single-copy DNA. The nucleotide sequences of the entire human prealbumin gene, including both 581 base pairs of the 5'- and 95 base pairs of the 3'-flanking sequences, were determined. The gene spans about 7.0 kilobase pairs and consists of four exons and three introns. As in most eukaryotic genes, the consensus TATA and CAAT sequences are found 30 and 101 nucleotides, respectively, upstream from the putative cap site, and a polyadenylation signal sequence AA-TAAA is found in the 3'-untranslated region. Unexpectedly, two independent open reading frames provided with respective regulatory sequences were found within the gene: one in the first intron and the other in the third intron.     

Characterization of novel and uncharacterized p53 SNPs in the Chinese population--intron 2 SNP co-segregates with the common codon 72 polymorphism

Multiple single nucleotide polymorphisms (SNPs) have been identified in the tumor suppressor gene p53, though the relevance of many of them is unclear. Some of them are also differentially distributed in various ethnic populations, suggesting selective functionality. We have therefore sequenced all exons and flanking regions of p53 from the Singaporean Chinese population and report here the characterization of some novel and uncharacterized SNPs - four in intron 1 (nucleotide positions 8759/10361/10506/11130), three in intron 3 (11968/11969/11974) and two in the 3'UTR (19168/19514). Allelic frequencies were determined for all these and some known SNPs, and were compared in a limited scale to leukemia and lung cancer patient samples. Intron 2 (11827) and 7 (14181/14201) SNPs were found to have a high minor allele frequency of between 26-47%, in contrast to the lower frequencies found in the US population, but similar in trend to the codon 72 polymorphism (SNP12139) that shows a distribution pattern correlative with latitude. Several of the SNPs were linked, such as those in introns 1, 3 and 7. Most interestingly, we noticed the co-segregation of the intron 2 and the codon 72 SNPs, the latter which has been shown to be expressed in an allele-specific manner, suggesting possible regulatory cross-talk. Association analysis indicated that the T/G alleles in both the co-segregating intron 7 SNPs and a 4tagSNP haplotype was strongly associated increased susceptibility to lung cancer in non-smoker females [OR: 1.97 (1.32, 3.394)]. These data together demonstrate high SNP diversity in p53 gene between different populations, highlighting ethnicity-based differences, and their association with cancer risk.     

Intron Distribution in Ancient Paralogs Supports Random Insertion and Not Random Loss

The intron positions of ten different protein families were examined to determine (the statistical likelihood of) whether spliceosomal introns are the result of random insertion events into previously intronless genes, on the one hand, or the result of random loss from common ancestral introns, on the other. The number of expected matches for the alternative scenarios was calculated for a binomial distribution by considering currently observed introns relative to all possible locations for insertion or loss. Introns occurring at approximately the same location (hereafter called a ``match') were tallied for each of the paired proteins. Matches were identified by their positions in the multiple alignment and were defined as any two introns occurring within a window of 11 possible nucleotide positions, thereby allowing for possible alignment errors and ``intron sliding.' Matches were tallied from the raw data and compared with the expected number of matches for the two different scenarios. The results suggest that the distribution of introns in genes encoding proteins is due to random insertion and not random loss. Received: 8 September 1996 / Accepted: 24 January 1997     

Characterization of mutations in the factor VIII gene by direct sequencing of amplified genomic DNA

In order to search for mutations resulting in hemophilia A that are not detectable by restriction analysis, three regions of the factor VIII gene were chosen for direct sequence analysis. Short segments of genomic DNA of 127 unrelated patients with hemophilia A were amplified by polymerase chain reaction. A total of 136,017 nucleotides were sequenced, and four mutations leading to the disease were found: a frameshift at codon 360 due to deletion of two nucleotides (GA), a nonsense codon 1705 due to a C----T transition, and two missense codons at positions 1699 and 1708. The first missense mutation (A----T) results in a Tyr----Phe substitution at a putative von Willebrand factor binding site. The second results in an Arg----Cys substitution at a thrombin cleavage site. In addition, we identified three rare sequence variants: a silent C----T transition at codon 34 which does not result in an amino acid change, a G----C change at codon 345 (Val----Leu), and an A----G change at the third nucleotide of intron 14. Direct sequence analysis of amplified DNA is a powerful but labor-intensive method of identifying mutations in large genes such as the human factor VIII gene.     

New polymorphic sites in the structure of the human phenylalanine hydroxylase gene

A previously unknown sequence of the human phenylalanine hydroxylase (PAH) gene intron 7 (GeneBank AN AF204239) has been reported. Screening of the group of phenylketonuria patients from Nobosibirsk region for polymorphic sites within intron 7 revealed single nucleotide substitutions at intron positions 332, 451, 574 and 791. Polymorphic site at intron position 791 corresponds to one of the eight restriction sites (MspI) utilized for haplotype construction. Analysis of the MspI allele frequencies in 29 phenylketonuria patients showed that the frequency of the MspI+ allele in this group was 79.4%. Polymorphic sites at nucleotide position +97 from the beginning of intron 10, and at nucleotide position -54 from the end of intron 5, were also described. The polymorphic sites revealed can be used as markers for identification of the PAH alleles in population genetic studies, and also serve for diagnostics of phenylketonuria (PKU). The presence of numerous nucleotide substitutions within the intronic sequences confirms highly polymorphic structure of the PAH gene.     

Genetic analysis of structure and function in phage T4 tRNASer

We have determined the nucleotide sequences of 55 spontaneous mutations that inactivate a suppressor gene of phage T4 tRNASer. Most of the mutations caused substitutions or deletions of single nucleotides at 18 different positions in the tRNA. Two of three mutations that allowed the synthesis of mature tRNA had nucleotide substitutions at the junction of the dihydrouridine and anticodon stems, suggesting that this region of tRNASer is important for aminoacylation. The third mutation that synthesized tRNA had a nucleotide deletion in the anticodon loop, which presumably affected the translational capacity of the tRNA. We also sequenced 58 spontaneous reversion mutations derived from strains with the inactive suppressor genes. Some of these regenerated the initial tRNA sequence, while other generated a second-site mutation in the tRNA. These second-site mutations restored helical base-pairings to the tRNA that had been eliminated by the initial mutations. The new base-pairings involved G.C and A.U, and the A.C wobble pair at certain positions in the tRNA. This finding establishes the existence of A.C wobble pair in tRNA helices.     

Analysis of sequences from the extremely A + T-rich genome of Plasmodium falciparum

The genome of the human malaria parasite Plasmodium falciparum has an A + T content of about 82%, higher than any other organism whose DNA has been characterized. Computer analysis of 36 kb of available nucleotide sequences from this species showed that the coding regions, with an A + T content of 69.0%, are flanked by more A + T-rich regions of 86.0% A + T. Within the coding sequences, the A/T ratio was 1.68 in the mRNA sense strand, and overall A + T content in the three codon positions increased in the order 1st-2nd-3rd position. Codons with T or especially A in the third position were strongly preferred. Codon usage among individual parasite genes was very similar compared to genes from other species. Dinucleotide frequencies for the parasite DNA were close to those expected for a random sequence with the known base composition, except that the CpG frequency in the coding sequences was low.     

Differential Mismatch Repair Can Explain the Disproportionalities between Physical Distances and Recombination Frequencies of Cyc1 Mutations in Yeast

Recombination rates have been examined in two-point crosses of various defined cyc1 mutations that cause the loss or nonfunction of iso-1-cytochrome c in the yeast Saccharomyces cerevisiae. Recombinants arising by three different means were investigated, including X-ray induced mitotic recombination, spontaneous mitotic recombination, and meiotic recombination. Heteroallelic diploid strains were derived by crossing cyc1 mutants containing a series of alterations at or near the same site to cyc1 mutants containing alterations at various distances. Marked disproportionalities between physical distances and recombination frequencies were observed with certain cyc1 mutations, indicating that certain mismatched bases can significantly affect recombination. The marker effects were more pronounced when the two mutational sites of the heteroalleles were within about 20 base pairs, but separated by at least 4 base pairs. Two alleles, cyc1-163 and cyc1-166, which arose by G.C----C.G transversions at nucleotide positions 3 and 194, respectively, gave rise to especially high rates of recombination. Other mutations having different substitutions at the same nucleotide positions were not associated with abnormally high recombination frequencies. We suggest that these marker effects are due to the lack of repair of either G/G or C/C mismatched base pairs, while the other mismatched base pair of the heteroallele undergoes substantial repair. Furthermore, we suggest that diminished recombination frequencies are due to the concomitant repair of both mismatches within the same DNA tract.     

On the frequencies of nucleotides and nucleotide substitutions in conservative regulatory DNA sequences

We present data on the frequencies of nucleotides and nucleotide substitutions in conservative DNA regions involved in the regulation of gene expression. Data on prokaryotes and eukaryotes are considered separately. In both cases DNA strands complementary to those which serve as templates for RNA-polymerase have low frequencies of cytosine. The most conservative positions also have an increased frequency of adenine. Various substitutions in the series of homologous regulatory DNA sequences, as compared to their consensuses, have different frequencies. In prokaryotes guanine in a consensus sequence is substituted for at the lowest and adenine at the highest frequency, whereas in eukaryotes cytosine is substituted for at the lowest and guanine at the highest frequency. In both cases the nucleotides substituted for are most frequently replaced with cytosine. Deviations from consensus sequences tend to cluster in adjacent positions. The more pronounced the consequences of a nucleotide substitution are the higher is the frequency of substitutions in adjacent positions. Possible explanations for these phenomena are discussed.     

Genomic DNA sequence of the cystic fibrosis transmembrane conductance regulator (CFTR) gene.

The gene responsible for cystic fibrosis, the most common severe autosomal recessive disorder, is located on the long arm of human chromosome 7, region q31-q32. The gene has recently been identified and shown to be approximately 250 kb in size. To understand the structure and to provide the basis for a systematic analysis of the disease-causing mutations in the gene, genomic DNA clones spanning different regions of the previously reported cDNA were isolated and used to determine the coding regions and sequences of intron/exon boundaries. A total of 22,708 bp of sequence, accounting for approximately 10% of the entire gene, was obtained. Alignment of the genomic DNA sequence with the cDNA sequence showed perfect colinearity between the two and a total of 27 exons, each flanked by consensus splice signals. A number of repetitive elements, including the Alu and Kpn families and simple repeats, such as (GT)17, (GATT)7, and (TA)14, were detected in close vicinity of some of the intron/exon boundaries. At least three of the simple repeats were found to be polymorphic in the population. Although an internal amino acid sequence homology could be detected between the two halves of the predicted polypeptide, especially in the regions of the two putative nucleotide-binding folds (NBF1 and NBF2), the lack of alignment of the nucleotide sequence as well as the different positions of the exon/intron boundaries does not seem to support the hypothesis of a recent gene duplication event. To facilitate detection of mutations by direct sequence analysis of genomic DNA, 28 sets of oligonucleotide primers were designed and tested for their ability to amplify individual exons and the immediately flanking sequences in the introns.     

Molecular evolution of eukaryotic genomes: hemiascomycetous yeast spliceosomal introns

As part of the exploratory sequencing program Génolevures, visual scrutinisation and bioinformatic tools were used to detect spliceosomal introns in seven hemiascomycetous yeast species. A total of 153 putative novel introns were identified. Introns are rare in yeast nuclear genes (<5% have an intron), mainly located at the 5′ end of ORFs, and not highly conserved in sequence. They all share a clear non-random vocabulary: conserved splice sites and conserved nucleotide contexts around splice sites. Homologues of metazoan snRNAs and putative homologues of SR splicing factors were identified, confirming that the spliceosomal machinery is highly conserved in eukaryotes. Several introns’ features were tested as possible markers for phylogenetic analysis. We found that intron sizes vary widely within each genome, and according to the phylogenetic position of the yeast species. The evolutionary origin of spliceosomal introns was examined by analysing the degree of conservation of intron positions in homologous yeast genes. Most introns appeared to exist in the last common ancestor of present day yeast species, and then to have been differentially lost during speciation. However, in some cases, it is difficult to exclude a possible sliding event affecting a pre-existing intron or a gain of a novel intron. Taken together, our results indicate that the origin of spliceosomal introns is complex within a given genome, and that present day introns may have resulted from a dynamic flux between intron conservation, intron loss and intron gain during the evolution of hemiascomycetous yeasts.