Rh gene evolution in primates: study of intron sequences

By amplification and sequencing of RH gene intron 4 of various primates we demonstrate that an Alu-Sx-like element has been inserted in the RH gene of the common ancestor of humans, apes, Old World monkeys, and New World monkeys. The study of mouse and lemur intron 4 sequences allowed us to precisely define the insertion point of the Alu-Sx element in intron 4 of the RH gene ancestor common to Anthropoidea. Like humans, chimpanzees and gorillas possess two types of RH intron 4, characterized by the presence (human RHCE and ape RHCE-like genes) or absence (human RHD and ape RHD-like genes) of the Alu-Sx element. This led us to conclude that in the RH common ancestor of humans, chimpanzees, and gorillas, a duplication of the common ancestor gene gave rise to two genes, one differing from the other by a 654-bp deletion encompassing an Alu-Sx element. Moreover, most of chimpanzees and some gorillas posses two types of RHD-like intron 4. The introns 4 of type 1 have a length similar to that of human RHD intron 4, whereas introns 4 of type 2 display an insertion of 12 bp. The latest insertion was not found in the human genome (72 individuals tested). The study of RH intron 3 length polymorphism confirmed that, like humans, chimpanzees and gorillas possess two types of intron 3, with the RHD-type intron 3 being 289 bases shorter than the RHCE intron 3. By amplification and sequencing of regions encompassing introns 3 and 4, we demonstrated that chimpanzee and gorilla RH-like genes displayed associations of introns 3 and 4 distinct to those found in man. Altogether, the results demonstrate that, as in humans, chimpanzee and gorilla RH genes experienced intergenic exchanges.     

Intron loss and gain during evolution of the catalase gene family in angiosperms.

Angiosperms (flowering plants), including both monocots and dicots, contain small catalase gene families. In the dicot, Arabidopsis thaliana, two catalase (CAT) genes, CAT1 and CAT3, are tightly linked on chromosome 1 and a third, CAT2, which is more similar to CAT1 than to CAT3, is unlinked on chromosome 4. Comparison of positions and numbers of introns among 13 angiosperm catalase genomic sequences indicates that intron positions are conserved, and suggests that an ancestral catalase gene common to monocots and dicots contained seven introns. Arabidopsis CAT2 has seven introns; both CAT1 and CAT3 have six introns in positions conserved with CAT2, but each has lost a different intron. We suggest the following sequence of events during the evolution of the Arabidopsis catalase gene family. An initial duplication of an ancestral catalase gene gave rise to CAT3 and CAT1. CAT1 then served as the template for a second duplication, yielding CAT2. Intron losses from CAT1 and CAT3 followed these duplications. One subclade of monocot catalases has lost all but the 5''-most and 3''-most introns, which is consistent with a mechanism of intron loss by replacement of an ancestral intron-containing gene with a reverse-transcribed DNA copy of a fully spliced mRNA. Following this event of concerted intron loss, the Oryza sativa (rice, a monocot) CAT1 lineage acquired an intron in a novel position, consistent with a mechanism of intron gain at proto-splice sites.     

Parental somatic and germ-line mosaicism for a multiexon deletion with unusual endpoints in a type III collagen (COL3A1) allele produces Ehlers-Danlos syndrome type IV in the heterozygous offspring.

Ehlers-Danlos syndrome (EDS) type IV is a dominantly inherited disorder that results from mutations in the type III collagen gene (COL3A1). We studied the structure of the COL3A1 gene of an individual with EDS type IV and that of her phenotypically normal parents. The proband was heterozygous for a 2-kb deletion in COL3A1, while her father was mosaic for the same deletion in somatic and germ cells. In fibroblasts from the father, approximately two-fifths of the COL3A1 alleles carried the deletion, but only 10% of the COL3A1 alleles in white blood cells were of the mutant species. The deletion in the mutant allele extended from intron 7 into intron 11. There was a 12-bp direct repeat in intron 7 and intron 11, the latter about 60 bp 5' to the junction. At the breakpoint there was a duplication of 10 bp from intron 11 separated by an insertion of 4 bp contained within the duplicated sequence. The father was mosaic for the deletion so that the gene rearrangement occurred during his early embryonic development prior to lineage allocation. These findings suggest that at least some of the deletions seen in human genes may occur during replication, rather than as a consequence of meiotic crossing-over, and that they thus have a risk for recurrence when observed de novo.     

Widespread distribution of a group I intron and its three deletion derivatives in the lysin gene of Streptococcus thermophilus bacteriophages

Of 62 Streptococcus thermophilus bacteriophages isolated from various ecological settings, half contain a lysin gene interrupted by a group IA2 intron. Phage mRNA splicing was demonstrated. Five phages possess a variant form of the intron resulting from three distinct deletion events located in the intron-harbored open reading frame (orf 253). The predicted orf 253 gene sequence showed a significantly lower GC content than the surrounding intron and lysin gene sequences, and the predicted protein shared a motif with endonucleases found in phages from both gram-positive and gram-negative bacteria. A comparison of the phage lysin genes revealed a clear division between intron-containing and intron-free alleles, leading to the establishment of a 14-bp consensus sequence associated with intron possession. The conserved intron was not found elsewhere in the phage or S. thermophilus bacterial genomes. Folding of the intron RNA revealed secondary structure elements shared with other phage introns: first, a 38-bp insertion between regions P3 and P4 that can be folded into two stem-loop structures (shared with introns from Bacillus phage SPO1 and relatives); second, a conserved P7.2 region (shared with all phage introns); third, the location of the stop codon from orf 253 in the P8 stem (shared with coliphage T4 and Bacillus phage SPO1 introns); fourth, orf 253, which has sequence similarity with the H-N-H motif of putative endonuclease genes found in introns from Lactococcus, Lactobacillus, and Bacillus phages.     

Comparative analysis of the structural organization of two closely related vitellogenin genes in X. laevis

The structural organization of the two closely related vitellogenin genes A1 and A2 has been determined and compared by electron microscopy. In both genes the mRNA-coding sequence of 6 kb is interrupted 33 times, leading to a total gene length of 21 kb for gene A1 and 16 kb for gene A2. Thus both genes have a mean exon length of 0.175 kb, while the mean intron length is 0.45 kb in gene A1 and 0.31 kb in gene A2. Because the introns interrupt the structural sequence at homologous positions in genes A1 and A2, we suggest that these two genes are the products of a duplication of an ancestral gene which had an intron-exon arrangement similar to that of the extant genes. Since the duplication event, the sequence and length of the analogous introns have changed rapidly, whereas homologous exons have diverged to an extent of only 5% of their sequences. The results suggest different mechanisms of evolution for exons and introns. While the exons evolved primarily by point mutations, such mutations, as well as deletion, insertion and duplication events, were important in the evolution of the introns.     

Intron-size constraint as a mutational mechanism in Rothmund-Thomson syndrome

Rothmund-Thomson syndrome (RTS) is an autosomal recessive disorder caused by deleterious mutations in the RECQL4 gene on chromosome 8. The RECQL4 gene structure is unusual because it contains many small introns <100 bp. We describe a proband with RTS who has a novel 11-bp intronic deletion, and we show that this mutation results in a 66-bp intron too small for proper splicing. Constraint on intron size may represent a general mutational mechanism, since human-genome analysis reveals that approximately 15% of genes have introns <100 bp and are therefore susceptible to size constraint. Thus, monitoring of intron size may allow detection of mutations missed by exon-by-exon approaches.     

Intron evolution in a phylogenetic perspective: Divergent trends in the two copies of the duplicated def gene in Impatiens L. (Balsaminaceae)

The history of MADS box genes is well-known in angiosperms. While duplication events and gene losses occur frequently, gene structure and intron positions are very conserved. We investigated all six introns in a duplicated MADS box gene (deficiens, def) in selected Impatiens taxa, thereby assessing intron features. For the first time, our study provides a comparison of molecular changes in all introns of a gene from a phylogenetic perspective. Interestingly, a uniform pattern of molecular evolution in the introns of each copy was not observed, but intron length increases, decreases, and size retention can be found in each copy. A tendency to accumulate long autapomorphic indels is also present, thus, a longer intron length does not reflect a higher number of parsimony-informative characters. Substitution rates vary between introns of each gene copy. While four of the six introns of def1 exhibit a change in their substitution rate, five of the six def2 introns maintain their rates throughout the genus albeit at different levels. In MADS box genes several regulatory sequences are found residing in introns. Thus, presence of putative regulatory motifs was investigated. Most of them are not conserved in position and usually present in only one of the gene copies. In addition, the potential for phylogenetic reconstruction of introns in both def copies is shortly discussed.     

Structure, chromosome location, and expression of the human smooth muscle (enteric type) gamma-actin gene: evolution of six human actin genes.

Recombinant phages that carry the human smooth muscle (enteric type) gamma-actin gene were isolated from human genomic DNA libraries. The amino acid sequence deduced from the nucleotide sequence matches those of cDNAs but differs from the protein sequence previously reported at one amino acid position, codon 359. The gene containing one 5' untranslated exon and eight coding exons extends for 27 kb on human chromosome 2. The intron between codons 84 and 85 (site 3) is unique to the two smooth muscle actin genes. In the 5' flanking region, there are several CArG boxes and E boxes, which are regulatory elements in some muscle-specific genes. Hybridization with the 3' untranslated region, which is specific for the human smooth muscle gamma-actin gene, suggests the single gene in the human genome and specific expressions in enteric and aortic tissues. From characterized molecular structures of the six human actin isoform genes, we propose a hypothesis of evolutionary pathway of the actin gene family. A presumed ancestral actin gene had introns at least sites 1, 2, and 4 through 8. Cytoplasmic actin genes may have directly evolved from it through loss of introns at sites 5 and 6. However, through duplication of the ancestral actin gene with substitutions of many amino acids, a prototype of muscle actin genes had been created. Subsequently, striated muscle actin and smooth muscle actin genes may have evolved from this prototype by loss of an intron at site 4 and acquisition of a new intron at site 3, respectively.     

Inactivation of human keratin genes: the spectrum of mutations in the sequence of an acidic keratin pseudogene

Keratins are cytoskeletal proteins encoded by a multigene family. We have identified the first human keratin pseudogene and determined its complete nucleotide sequence. Sequence comparisons indicate that the pseudogene arose from a very recent duplication of the 50-kd keratin (K14) gene. The coding and the intron sequences of the two genes are 95% and 93% identical, respectively. Although the sequence of the regulatory region in the pseudogene is virtually identical to that in the 50-kd functional gene, several deleterious mutations have been identified in the pseudogene. There are three frameshifts in the coding regions, one of which is a perfect 8-bp duplication. A single-base-pair deletion in the first exon and a single-base-pair insertion in the penultimate exon also result in frameshifts. The three remaining deleterious mutations interfere with the mRNA processing signals: two alter the intron/exon boundaries, and the third disrupts the polyadenylation signal. These mutations clearly identify the sequence as a human keratin pseudogene.     

Elimination of a group II intron from a plastid gene causes a mutant phenotype

Group II introns are found in bacteria and cell organelles (plastids, mitochondria) and are thought to represent the evolutionary ancestors of spliceosomal introns. It is generally believed that group II introns are selfish genetic elements that do not have any function. Here, we have scrutinized this assumption by analyzing two group II introns that interrupt a plastid gene (ycf3) involved in photosystem assembly. Using stable transformation of the plastid genome, we have generated mutant plants that lack either intron 1 or intron 2 or both. Interestingly, the deletion of intron 1 caused a strong mutant phenotype. We show that the mutants are deficient in photosystem I and that this deficiency is directly related to impaired ycf3 function. We further show that, upon deletion of intron 1, the splicing of intron 2 is strongly inhibited. Our data demonstrate that (i) the loss of a group II intron is not necessarily phenotypically neutral and (ii) the splicing of one intron can depend on the presence of another.     

Chromosomal arrangement of the duck alpha-globin genes and primary structure of the embryonic alpha-globin gene pi

A recombinant lambda Charon 4A bacteriophage, D alpha G-1, carrying the genes coding for the duck embryonic (pi') and adult (alpha A, alpha D) alpha-like globins was isolated from a previously constructed duck DNA recombinant library. The three globin genes are transcribed from the same DNA strand and are arranged in the order of their expression during development: 5'-pi'-alpha D-alpha A-3'. We have determined the complete nucleotide sequence of the duck pi'-globin gene, including the flanking regions. Due to the unusual length of intron 1 (963 bp) and intron 2 (568 bp) the 2167-bp duck pi'-globin gene is by far the largest among all known mammalian or avian alpha- and beta-globin genes. For instance, the duck pi'-globin gene introns are almost twice as long as those of the chicken pi'-globin genes. A surprisingly high degree of nucleotide sequence homology (88%) has been found for the 5' flanking region (positions -1 to -223) of the duck and chicken pi'-globin gene.     

Prevalence of intron gain over intron loss in the evolution of paralogous gene families

The mechanisms and evolutionary dynamics of intron insertion and loss in eukaryotic genes remain poorly understood. Reconstruction of parsimonious scenarios of gene structure evolution in paralogous gene families in animals and plants revealed numerous gains and losses of introns. In all analyzed lineages, the number of acquired new introns was substantially greater than the number of lost ancestral introns. This trend held even for lineages in which vertical evolution of genes involved more intron losses than gains, suggesting that gene duplication boosts intron insertion. However, dating gene duplications and the associated intron gains and losses based on the molecular clock assumption showed that very few, if any, introns were gained during the last ~100 million years of animal and plant evolution, in agreement with previous conclusions reached through analysis of orthologous gene sets. These results are generally compatible with the emerging notion of intensive insertion and loss of introns during transitional epochs in contrast to the relative quiet of the intervening evolutionary spans.     

Structural organization of the 5' region of the thyroglobulin gene. Evidence for intron loss and "exonization" during evolution

More than one third of thyroglobulin (1190 residues out of 2750) is made of one peptide motif repeated ten times in tandem. Segments unrelated to the motif interrupt this structure at various places. The corresponding gene region, which extends over 40 x 10(3) bases, was studied in detail. All exon borders and exon/intron junctions were localized precisely and sequenced, and their positions were correlated with the repetitive organization of the protein. When intron positions were compiled on a consensus sequence of all repeats, three categories of introns were observed. Except between repeats numbers 5 and 6, an intron was invariably found within the Cys codon making the limit of each motif. This category of intron most probably reflects the serial duplication events responsible for the evolution of this region of the gene. All other introns, except no. 2, are found at positions were the repetitive structure is disrupted by "inserted" peptides. We present the hypothesis that this second category of introns was already present in the original unit before the first duplication. Thereafter, they would have experienced either complete loss (some units do not contain any intron) or partial or total exonization, resulting in the slipping of intronic material into coding sequence. Intron no. 2, finally, separates motif no. 1 at a position on the boundary between two segments presenting sequence homology. This last type of intron probably reflects an initial duplication event at the origin of a primordial thyroglobulin gene motif. With all these characteristics, the thyroglobulin gene is presented as a paradigm for the analysis of the fate of introns in gene evolution.     

High rate of recent intron gain and loss in simultaneously duplicated Arabidopsis genes

We examined the gene structure of a set of 2563 Arabidopsis thaliana paralogous pairs that were duplicated simultaneously 20-60 MYA by tetraploidy. Out of a total of 23,164 introns in these genes, we found that 10,004 pairs have been conserved and 578 introns have been inserted or deleted in the time since the duplication event. This intron insertion/deletion rate of 2.7 x 10(-3) to 9.1 x 10(-4) per site per million years is high in comparison to previous studies. At least 56 introns were gained and 39 lost based on parsimony analysis of the phylogenetic distribution of these introns. We found weak evidence that genes undergoing intron gain and loss are biased with respect to gene ontology terms. Gene pairs that experienced at least 2 intron insertions or deletions show evidence of enrichment for membrane location and transport and transporter activity function. We do not find any relationship of intron flux to expression level or G + C content of the gene. Detection of a bias in the location of intron gains and losses within a gene depends on the method of measurement: an intragene method indicates that events (specifically intron losses) are biased toward the 3' end of the gene. Despite the relatively recent acquisition of these introns, we found only one case where we could identify the mechanism of intron origin--the TOUCH3 gene has experienced 2 tandem, partial, internal gene duplications that duplicated a preexisting intron and also created a novel, alternatively spliced intron that makes use of a duplicated pair of cryptic splice sites.     

Conservation of the positions of metazoan introns from sponges to humans

Sponges (phylum Porifera) are the phylogenetic oldest Metazoa still extant. They can be considered as reference animals (Urmetazoa) for the understanding of the evolutionary processes resulting in the creation of Metazoa in general and also for the metazoan gene organization in particular. In the marine sponge Suberites domuncula, genes encoding p38 and JNK kinases contain nine and twelve introns, respectively. Eight introns in both genes share the same positions and the identical phases. One p38 intron slipped for six bases and the JNK gene has three more introns. However, the sequences of the introns are not conserved and the introns in JNK gene are generally much longer. Introns interrupt most of the conserved kinase subdomains I-XI and are found in all three phases (0, 1 and 2). We analyzed in details p38 and JNK genes from human, Caenorhabditis elegans and Drosophila melanogaster and found in most genes introns at the positions identical to those in sponge genes. The exceptions are two p38 genes from D. melanogaster that have lost all introns in the coding sequence. The positions of 11 introns in each of four human p38 genes are fully conserved and ten introns occupy identical positions as the introns in sponge p38 or JNK genes. The same is true for nine, out of ten introns in the human JNK-1 gene. The introns in human p38 and JNK genes are on average more than ten times longer than corresponding introns in sponges. It was proposed that yeast HOG1-like kinases (from i.e. Saccharomyces cerevisiae and Emericella nidulans) and metazoan p38 and JNK kinases are orthologues. p38 and JNK genes were created after the split from fungi by the duplication and diversification of the HOG1-like progenitor gene. Our results further support the common origin of p38 and JNK genes and speak in favor of a very early time of duplication. The ancestral gene contained at least ten introns, which are still present at the very conserved positions in p38 and JNK genes of extant animals. Four of these introns are present at the same positions in the HOG-like gene in the fungus E. nidulans. The others probably entered the ancestral gene after the split of fungi, but before the duplication of the gene and before the creation of the common, urmetazoan progenitor of all multicellular animals. A second gene coding for an immune molecule is described, the allograft inflammatory factor, which likewise showed a highly conserved exon/intron structure in S. domuncula and in human. These data show that the intron/exon borders are highly conserved in genes from sponges to human.     

Positional cloning of rice semidwarfing gene, sd-1: rice "green revolution gene" encodes a mutant enzyme involved in gibberellin synthesis.

A rice semidwarfing gene, sd-1, known as the "green revolution gene," was isolated by positional cloning and revealed to encode gibberellin 20-oxidase, the key enzyme in the gibberellin biosynthesis pathway. Analysis of 3477 segregants using several PCR-based marker technologies, including cleaved amplified polymorphic sequence, derived-CAPS, and single nucleotide polymorphisms revealed 1 ORF in a 6-kb candidate interval. Normal-type rice cultivars have an identical sequence in this region, consisting of 3 exons (558, 318, and 291 bp) and 2 introns (105 and 1471 bp). Dee-Geo-Woo-Gen-type sd-1 mutants have a 383-bp deletion from the genome (278-bp deletion from the expressed sequence), from the middle of exon 1 to upstream of exon 2, including a 105-bp intron, resulting in a frame-shift that produces a termination codon after the deletion site. The radiation-induced sd-1 mutant Calrose 76 has a 1-bp substitution in exon 2, causing an amino acid substitution (Leu [CTC] to Phe [TTC]). Expression analysis suggests the existence of at least one more locus of gibberellin 20-oxidase which may prevent severe dwarfism from developing in sd-1 mutants.