Gene Expression,Intron Density,and Splice Site Strength in <Emphasis Type="Italic">Drosophila</Emphasis> and <Emphasis Type="Italic">Caenorhabditis</Emphasis>

In this paper we investigate the relationships among intron density (number of introns per kilobase of coding sequence), gene expression level, and strength of splicing signals in two species: Drosophila melanogaster and Caenorhabditis elegans. We report a negative correlation between intron density and gene expression levels, opposite to the effect previously observed in human. An increase in splice site strength has been observed in long introns in D. melanogaster. We show this is also true of C. elegans. We also examine the relationship between intron density and splice site strength. There is an increase in splice site strength as the intron structure becomes less dense. This could suggest that introns are not recognized in isolation but could function in a cooperative manner to ensure proper splicing. This effect remains if we control for the effects of alternative splicing on splice site strength. Reviewing Editor: Dr. Nicolas Galtier     

Novel Group I Introns Encoding a Putative Homing Endonuclease in the Mitochondrial <Emphasis Type="Italic">cox1</Emphasis> Gene of Scleractinian Corals

Analyses of mitochondrial sequences revealed the existence of a group I intron in the cytochrome oxidase subunit 1 (cox1) gene in 13 of 41 genera (20 out of 73 species) of corals conventionally assigned to the suborder Faviina. With one exception, phylogenies of the coral cox1 gene and its intron were concordant, suggesting at most two insertions and many subsequent losses. The coral introns were inferred to encode a putative homing endonuclease with a LAGLI-DADG motif as reported for the cox1 group I intron in the sea anemone Metridium senile. However, the coral and sea anemone cox1 group I introns differed in several aspects, such as the intron insertion site and sequence length. The coral cox1 introns most closely resemble the mitochondrial cox1 group I introns of a sponge species, which also has the same insertion site. The coral introns are also more similar to the introns of several fungal species than to that of the sea anemone (although the insertion site differs in the fungi). This suggests either a horizontal transfer between a sponge and a coral or independent transfers from a similar fungal donor (perhaps one with an identical insertion site that has not yet been discovered). The common occurrence of this intron in corals strengthens the evidence for an elevated abundance of group I introns in the mitochondria of anthozoans. [Reviewing Editor: Dr. Niles Lehman]     

Intron and exon length variation in <Emphasis Type="Italic">Arabidopsis</Emphasis>, rice,nematode, and human

As determined by computer sequence analysis, the average exon length in Arabidopsis thaliana, Oryza sativa, Caenorhabditis elegans, and Homo sapiens genes decreases with an increasing number of introns. In A. thaliana and O. sativa, variations in intron and exon lengths with an increasing number of introns are highly correlated. Linear correlation is observed between the total exon length and the number of introns, while the gene length increases in proportion to the number of introns. In human, the average intron and gene lengths depended on the gene density in DNA.     

Correlated Evolution of Synonymous and Nonsynonymous Sites in <Emphasis Type="Italic">Drosophila</Emphasis>

Recent work has shown that Drosophila melanogaster genes with fast-evolving nonsynonymous sites have lower codon usage bias. This pattern has been attributed to interference between positive selection at nonsynonymous sites and weak selection on codon usage. Here we have looked for this correlation in a much larger and less biased dataset, comprising 630 gene pairs from D. melanogaster and D. yakuba. We confirmed that there is a negative correlation between the rate of nonsynonymous substitutions (d_N) and codon bias in D. melanogaster. We then tested the interference hypothesis and other alternative explanations, including one involving gene expression. We found that d_N indeed correlates with the level of gene expression. Given that gene expression is a strong determinant of codon bias, the relationship between d_N and codon bias might be a by-product of gene expression. However, our tests show that none of the hypotheses we consider seem to explain the data fully.This article contains online supplementary material.Reviewing Editor: Dr. John Huelsenbeck     

Variability of Nuclear SSU-rDNA Group Introns Within Septoria Species: Incongruence with Host Sequence Phylogenies

We report structural features and distribution patterns of 26 different group I introns located at three distinct nucleotide positions in nuclear small subunit ribosomal DNA (SSU-rDNA) of 10 Septoria and 4 other anamorphic species related to the teleomorphic genus Mycosphaerella. Secondary structure and sequence characteristics assigned the introns to the common IC1 and IE groups. Intron distribution patterns and phylogenetic relationships strongly suggested that some horizontal transfer events have occurred among the closely related fungal species sampled. To test this hypothesis, we used a comparative approach of intron- and rDNA-based phylogenies through MP- and ML-based topology tests. Our results showed two statistically well-supported major incongruences between the intron and the equivalent internal transcribed spacer (ITS) tree comparisons made. Such absence of a co-evolutive history between group I introns and host sequences is discussed relatively to the intron structures, the mechanisms of intron movement, and the biology of the Mycosphaerella pathogenic fungi. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Reviewing Editor: Debashish Bhattacharya     

Isochore Structures in the Genome of the Plant <Emphasis Type="BoldItalic">Arabidopsis thaliana</Emphasis>

Arabidopsis thaliana is an important model system for the study of plant biology. We have analyzed the complete genome sequences of Arabidopsis by using a newly developed windowless method for the GC content computation, the cumulative GC profile. It is shown that the Arabidopsis genome is organized into a mosaic structure of isochores. All the centromeric regions are located in GC-rich isochores, called centromere-isochores, which are characterized by a high GC content but low gene and T-DNA insertion densities. This characteristic distinguishes centromere-isochores from the other class of GC-rich isochores, called GC-isochores, which have high gene and T-DNA insertion densities. Consequently, 15 isochores have been identified, i.e., 7 AT-isochores, 3 GC-isochores, and 5 centromere-isochores. The genes in centromere-isochores, which have the highest GC content, have much shorter intron lengths and lower intron numbers, compared to those of the other two types. There is also considerable difference in the numbers and lengths of transposable elements (TEs) between AT and GC-isochores, i.e., the TE number (length) of AT-isochores is 6.3 (7.3) times that of GC-isochores. It is generally believed that TEs are accumulated in the regions surrounding the centromeres. However, within these TE-rich regions, there are regions of extremely low TE numbers (TE deserts), which correspond to the positions of centromere-isochores. In addition, a heterochromatic knob is located at the boundary of an AT-isochore. Furthermore, we show that the differences in GC content among isochores are mainly due to the GC content variation of introns, the third codon positions and intergenic regions.[Reviewing Editor: Martin Kreitman]     

Molecular Basis of the Size Polymorphism of the First Intron of the<Emphasis Type="Italic">Adh-1</Emphasis> Gene of the Mediterranean Fruit Fly, <Emphasis Type="Italic">Ceratitis capitata</Emphasis>

The first intron of the gene encoding one of the alcohol dehydrogenase isoenzymes (ADH-1) in Ceratitis capitata is highly polymorphic in size. Five size variants of this intron were isolated from different strains and populations and characterized. Restriction map and sequence analysis showed that the intron size polymorphism is due to the presence or absence of (a) a copy of a defective mariner-like element, postdoc; (b) an 550-bp 3 indel which exhibits no similarity to any known sequence; and (c) a central duplication of 704 bp consisting of part of the 3 end of the postdoc element, the region between postdoc and the 3 indel, and the first 20 bp of the 3 indel. The homologous Adh-1 intron was amplified from the congeneric species, Ceratitis rosa, in order to obtain an outgroup for comparative and phylogenetic analyses. The C. rosa introns were polymorphic in size, ranging from about 1100 to 2000 bp, the major difference between them being the presence or absence of a mariner-like element Crmar2, unrelated to the postdoc element. Phylogenetic analysis suggests that the shorter intron variants in C. capitata may represent the ancestral form of the intron, the longest variants apparently being the most recent.     

Patterns of intron sequence evolution in Drosophila are dependent upon length and GC content

Background

Introns comprise a large fraction of eukaryotic genomes, yet little is known about their functional significance. Regulatory elements have been mapped to some introns, though these are believed to account for only a small fraction of genome wide intronic DNA. No consistent patterns have emerged from studies that have investigated general levels of evolutionary constraint in introns.

Results

We examine the relationship between intron length and levels of evolutionary constraint by analyzing inter-specific divergence at 225 intron fragments in Drosophila melanogaster and Drosophila simulans, sampled from a broad distribution of intron lengths. We document a strongly negative correlation between intron length and divergence. Interestingly, we also find that divergence in introns is negatively correlated with GC content. This relationship does not account for the correlation between intron length and divergence, however, and may simply reflect local variation in mutational rates or biases.

Conclusion

Short introns make up only a small fraction of total intronic DNA in the genome. Our finding that long introns evolve more slowly than average implies that, while the majority of introns in the Drosophila genome may experience little or no selective constraint, most intronic DNA in the genome is likely to be evolving under considerable constraint. Our results suggest that functional elements may be ubiquitous within longer introns and that these introns may have a more general role in regulating gene expression than previously appreciated. Our finding that GC content and divergence are negatively correlated in introns has important implications for the interpretation of the correlation between divergence and levels of codon bias observed in Drosophila.     

Patterns of Base Composition Within the Genes of Drosophila melanogaster

Base composition is not uniform across the genome of Drosophila melanogaster. Earlier analyses have suggested that there is variation in composition in D. melanogaster on both a large scale and a much smaller, within-gene, scale. Here we present analyses on 117 genes which have reliable intron/exon boundaries and no known alternative splicing. We detect significant heterogeneity in G+C content among intron segments from the same gene, as well as a significant positive correlation between the intron and the third codon position G+C content within genes. Both of these observations appear to be due, in part, to an overall decline in intron and third codon position G+C content along Drosophila genes with introns. However, there is also evidence of an increase in third codon position G+C content at the start of genes; this is particularly evident in genes without introns. This is consistent with selection acting against preferred codons at the start of genes. Received: 24 February 1997 / Accepted: 10 November 1997     

Nature of allelic sequence polymorphism at the KIR3DL3 locus

KIR3DL3 is a framework gene of the Leukocyte Receptor Complex, present in all individuals and haplotypes analysed to date. We describe 17 novel KIR3DL3 alleles, including seven single nucleotide polymorphic (SNP) positions within the coding region. Sequence variation within introns included a VNTR within intron 1. As KIR3DL3 mRNA is known to be expressed in decidual NK cells, we investigated the impact of KIR3DL3 allelic variation on pre-eclampsia. No statistical difference in allele frequency or polymorphism was observed between pre-eclampsia patient and control cohorts. Linkage disequilibrium (LD) analysis of exonic SNPs suggested that recombination may be a mechanism of generating sequence diversity within KIR3DL3. A potential recombination hotspot was located within intron 5. A strong LD was detected between polymorphism in exon 6 of KIR3DL3 and the KIR gene −2DL3 or -2DS2 loci, which define the centromeric end of two main haplotypes (A and B) of the KIR cluster. Comparison of primate KIR sequences indicated that the Ig domains of KIR3DL3 are highly conserved between chimpanzee, gorilla and humans. Investigation of KIR3DL3 dN/dS ratios indicated a greater level of synonymous mutations consistent with purifying selection, although positive selection was detected acting on two sites within the stem region.Electronic supplementary material Supplementary material is available for this article at and is accessible for authorized users.     

Gene structure of two-domain arginine kinases from Anthopleura japonicus and Pseudocardium sachalinensis

Unusual two-domain arginine kinases (AKs) arose independently at least two times during molecular evolution of phosphagen kinases: AKs from the primitive sea anemone Anthopleura japonicus and from the clam Pseudocardium sachalinensis. To elucidate its unusual evolution, the structures of Anthopleura and Pseudocardium AK genes have been determined. The Anthopleura gene consisted of 4 exons and 3 introns: two domains are linked by a bridge intron, and each domain contains one intron in different positions. On the other hand, the Pseudocardium gene consisted of 10 exons and 9 introns: two domains are also linked by a bridge intron, and domains 1 and 2 contains 3 and 5 introns, respectively, of which 3 introns are located in exactly same positions. Since the two domains of Pseudocardium AK are estimated to have diverged about 290 million years ago, the 3 introns have been conserved at least for this long. Comparison of intron positions in Anthopleura, Pseudocardium and C. elegans AK genes indicates that there is no intron conserved through the three AK lineages, in sharp contrast to relatively conservative intron positions in creatine kinase (CK) gene family.     

Genome-scale identification of resistance gene analogs and the development of their intron length polymorphism markers in maize

As introns are vulnerable to changes such as insertions and deletions when exposed to various evolutionary forces, they constitute a repository for developing genetic markers based on intron length polymorphisms (ILP). This study developed a set of genetic markers that use the potential intron length polymorphism in resistance gene analogs (RGAs) in Zea mays. By searching the genome of Zea mays B73 for the homologs of 73 R genes which have already been identified in plants, we found 861 RGAs, 632 of which have at least one intron that can serve as putative markers targeting the intron length polymorphism in RGAs (RGA-ILP). We developed 1972 candidate markers via electronic PCR (e-PCR) with primer pairs designed in each pair of exonic regions that flank an intron. Furthermore, the performance of RGA-ILP among four maize inbred lines (Huangzao4, B73, Mo17, and Dan340) was evaluated with 69 pairs of randomly selected primers. Of them, 46.4% showed bands that had discriminating length polymorphism, and between any two of the inbred lines the proportion of polymorphism ranged from 23.2 to 31.9%. To make it convenient to use these markers for those interested in molecular breeding of disease-resistant maize, we provide all related information in a web-based database named MaizeRGA, which is available at .     

Evolution of introns and exons of class II major histocompatibility complex genes of vertebrates

 The phylogenetic relationships and patterns of nucleotide substitution were compared for introns and exons of class II major histocompatibility complex (MHC) genes in three datasets: human DRB1, human DQA1, and cyprinid fish DAB1. In both human DRB1 and cyprinid DAB1, there was strong evidence that recombination events between alleles have occurred in such a way that intron and exon sequences of a given allele do not necessarily share the same evolutionary history. In the case of human DRB1, recombination was found to have homogenized intron 1 and intron 2 sequences relative to exon 2 sequences within lineages of alleles but not between lineages. As a result, mean divergence times of intron sequences are much more recent than those of exonic sequences. Thus, the divergence time of DRB1 introns cannot be used to date that of exons in the same alleles, and the hypothesis that most human DRB1 polymorphism is of very recent origin is not supported. Received: 5 September 1999 / Revised: 30 December 1999     

Complex selection on intron size in Cryptococcus neoformans

We conducted a genome-wide analysis of the roles of mutation and selection in sculpting intron size in the fungal pathogen Cryptococcus neoformans. We find that deletion rate is positively associated with intron length and that insertion rate exhibits a weak negative association with intron length. These patterns suggest that long introns as well as extremely short introns in this unusually intron-rich fungal genome are in mutation-selection disequilibrium and that the proportion of constrained functional sequence in introns does not scale linearly with size. We find that untranslated region introns are longer than coding-region introns and that first introns are substantially longer than subsequent introns, suggesting heterogeneous distribution of constrained functional sequence and/or selective pressures on intron size within genes. In contrast to Drosophila, we find a positive correlation between d(N) and first intron or last intron length and a negative correlation between d(N) and internal intron length. This contrasting pattern may indicate that terminal introns and internal introns are differentially subject to hypothesized selection pressures modulating intron size and provides further evidence of widespread selective constraints on noncoding sequences.     

DETECTION OF A GROUP I (IE) FUNGAL INTRON IN THE GREEN ALGAL GENUS UROSPORA (ULVOPHYCEAE)1

A group IE intron at position 989 (Escherichia coli numbering) in the 18S rRNA gene is described for the first time from the green algae based on primary and secondary structure analysis of the intron. The intron occurred in Urospora neglecta Kornmann, U. wormskioldii (Mertens) Rosenvinge, and in an undescribed taxon closely allied with U. wormskioldii. The Urospora intron was 397 or 399 bp in length, except for one variant, which was 652 bp in length. The latter intron appeared to be a chimera, consisting of a conspecific part, a segment more closely related to another species, and a region of unknown origin. A search of data bases containing introns indicated that the green algal intron is homologous with other 989 introns, all occurring in fungi, which may be the source of the Urospora intron. The introduction of this intron into the Urospora lineage has occurred within the last 14 million years, based on estimates of the earliest divergence of the genus from other members of the family Acrosiphoniaceae, all of which lack 18S introns.     

Evolution of <Emphasis Type="Italic">Pleopsidium</Emphasis> (Lichenized Ascomycota) S943 Group I Introns and the Phylogeography of an Intron-Encoded Putative Homing Endonuclease

The sporadic distribution of nuclear group I introns among different fungal lineages can be explained by vertical inheritance of the introns followed by successive losses, or horizontal transfers from one lineage to another through intron homing or reverse splicing. Homing is mediated by an intron-encoded homing endonuclease (HE) and recent studies suggest that the introns and their associated HE gene (HEG) follow a recurrent cyclical model of invasion, degeneration, loss, and reinvasion. The purpose of this study was to compare this model to the evolution of HEGs found in the group I intron at position S943 of the nuclear ribosomal DNA of the lichen-forming fungus Pleopsidium. Forty-eight S943 introns were found in the 64 Pleopsidium samples from a worldwide screen, 22 of which contained a full-length HEG that encodes a putative 256-amino acid HE, and 2 contained HE pseudogenes. The HEGs are divided into two closely related types (as are the introns that encode them) that differ by 22.6% in their nucleotide sequences. The evolution of the Pleopsidium intron-HEG element shows strong evidence for a cyclical model of evolution. The intron was likely acquired twice in the genus and then transmitted via two or three interspecific horizontal transfers. Close geographical proximity plays an important role in intron-HEG horizontal transfer because most of these mobile elements were found in Europe. Once acquired in a lineage, the intron-HEG element was also vertically transmitted, and occasionally degenerated or was lost. [Reviewing Editor: Dr. Manyuan Long]     

Genetic diversity of dodder (<Emphasis Type="Italic">Cuscuta</Emphasis> spp.) collected from commercial cranberry production as revealed in the<Emphasis Type="Italic">trn</Emphasis>L (UAA) intron

Swamp dodder (Cuscuta gronovii) is a parasitic plant detrimental to cranberries. Observation of emergence of dodder seeds collected from a cultivated cranberry bog in Massachusetts revealed 2 or more peak emergence times during 4 consecutive growing seasons. Molecular methods were used to investigate genetic variation among the emerging dodder seedlings. On emergence, dodder seedlings were collected and analyzed for DNA sequence diversity in thetrnL (UAA) intron, a noncoding region of chloroplast DNA. DNA sequence analysis of 87 seedlings collected during the 1999 and 2000 growing seasons revealed the presence of 2 dodder ecotypes, designated A and B. Comparative DNA sequence analysis indicated that in thetrnL (UAA) intron, the sequence of ecotype A is identical to that ofCuscuta gronovii, whereas the sequence of ecotype B is closest to that ofCuscuta attenuata (99.3% sequence identity; 293 bases considered). ABg/II restriction enzyme cut site was identified that distinguished between thetrnL (UAA) introns of ecotypes A and B. Restriction fragment length polymorphism (RFLP) was used to analyze the sequences of 100 seedlings collected during the growing seasons of 2001 and 2002. Only 10 of the 187 samples were ecotype A, all of which emerged on or before May 7 in the growing seasons. Therefore, the predominant dodder haplotype found in this study may be a close relative ofC. attenuata and notC. gronovii, the common species found in cranberry bogs.     

Low temperature promotes intron retention in two <Emphasis Type="Italic">e-cor</Emphasis> genes of durum wheat

Following the screening of a suppression subtractive library developed from durum wheat plants exposed to low temperature for 6 h, two early cold-regulated (e-cor) genes have been isolated. These genes, coding putatively for a ribokinase (7H8) and a C3H2C3 RING-finger protein (6G2), were characterized by the stress-induced retention of a subset of introns in the mature mRNA. This feature was dependent on cold for 7H8 and on cold and dehydration for 6G2. When other genes, such as the stress-related gene WCOR410c, coding for a dehydrin (one intron), or a gene coding for a putative ATP binding cassette transporter (16 introns) were analyzed, no cold-dependent intron retention was observed. Cold-induced intron retention was not observed in mutants defective in the chloroplast development; nevertheless treatment with cycloheximide in the absence of cold was able to promote intron retention for the 7H8 e-cor gene. These results suggest that the cold-induced intron retention reflects the response of the spliceosoma to specific environmental signals transduced to the splicing protein factors through a chloroplast-dependent pathway. Notably, when the 7H8 Arabidopsis orthologous gene was analyzed, no stress induction in terms of mRNA abundance and no cold-dependent intron retention was detected. Otherwise, 6G2 Arabidopsis homologous sequences sharing the same genomic structure of the durum wheat 6G2 showed a similar intron retention event although not strictly dependent on stress.Electronic Supplementary Material Supplementary material is available for this article at