Reconstruction of the phylogenetic position of larch (Larix sukaczewii Dylis) by sequencing data for the trnK intron of chloroplast DNA

To reconstruct the systematic relationships of larch Larix sukaczewii, we used the chloroplast trnK intron sequences of L. decidua, L. sukaczewii, L. sibirica, L. czekanovskii, and L. gmelinii. Analysis of phylogenetic trees constructed using the maximum parsimony and maximum likelihood methods showed a clear divergence of the trnK intron sequences between L. sukaczewii and L. sibirica. This divergence reaches intraspecific level, which supports a previously published hypothesis on the taxonomic isolation of L. sukaczewii.     

Caryophyllales: Evaluating phylogenetic signal in trnK intron versus matK

Abstract We assess the phylogenetic information in trnK intron at the ordinal level using the Caryophyllales and compare it with that derived from matK. The trnK gene is split into two exons by an intron that includes the matK gene. The plastid trnK is a tRNA gene encoding Lysine^(UUU), whereas the matK gene is a putative group II intron maturase. The two regions are usually coamplified, and trnK intron is partially sequenced but its sequences are often excluded from phylogenetic reconstruction at deep historic levels. This study shows that the two regions are comparable in proportion of variable sites, possess a comparable pattern of substitution rates per site, and display similar phylogenetic informativeness profiles and per‐site informativeness. Phylogenetic analyses show strong congruence between phylogenetic trees based on matK and trnK intron partitioned datasets from 45 genera representing 30 of the 34 recognized Caryophyllales families. The trnK intron alone provides a relatively well‐resolved topology for the order. Combining the trnK intron with matK sequence data resulted in six most parsimonious trees, differing only in the placement of Claytonia (Portulacaceae) within the noncore group. A well‐supported major basal split in the order into core and noncore Caryophyllales with Rhabdodendraceae, Simmondsiaceae, and Asteropeiaceae as sister to remaining core lineages is evident in partitioned and combined analyses. The placement of these three families has been disputable, impacting the overall backbone topology of the Caryophyllales. This study demonstrates the cost effectiveness of using the trnK intron along with matK (both substitutions and insertions/deletions) at deeper phylogenetic level.     

Geographical origin of Leucobryum boninense Sull. & Lesq. (Leucobryaceae,Musci) endemic to the Bonin Islands,Japan

Leucobryum boninense is endemic to the Bonin Islands, Japan, and its related species are widely distributed in Asia and the Pacific. We aimed to clarify the phylogenetic relationships among Leucobryum species and infer the origin of L. boninense. We also describe the utility of the chloroplast trnK intron including matK for resolving the phylogenetic relationships among Leucobryum species, as phylogenetic analyses using trnK intron and/or matK have not been performed well in bryophytes to date. Fifty samples containing 15 species of Leucobryum from Asia and the Pacific were examined for six chloroplast DNA regions including rbcL, rps4, partial 5′ trnK intron, matK, partial 3′ trnK intron, and trnL‐F intergenic spacer plus one nuclear DNA region including ITS. A molecular phylogenetic tree showed that L. boninense made a clade with L. scabrum from Japan, Taiwan and, Hong Kong; L. javense which is widely distributed in East and Southeast Asia, and L. pachyphyllum and L. seemannii restricted to the Hawaii Islands, as well as with L. scaberulum from the Ryukyus, Japan, Taiwan, and southeastern China. Leucobryum boninense from various islands of the Bonin Islands made a monophylic group that was closely related to L. scabrum and L. javense from Japan. Therefore, L. boninense may have evolved from L. scabrum from Japan, Taiwan, or Hong Kong, or L. javense from Japan. We also described the utility of trnK intron including matK. A percentage of the parsimony‐informative characters in trnK intron sequence data (5.8%) was significantly higher than that from other chloroplast regions, rbcL (2.4%) and rps4 (3.2%) sequence data. Nucleotide sequence data of the trnK intron including matK are more informative than other chloroplast DNA regions for identifying the phylogenetic relationships among Leucobryum species.     

Phylogenetic Use of Noncoding Regions in the GenusGentianaL.: ChloroplasttrnL (UAA) Intron versus Nuclear Ribosomal Internal Transcribed Spacer Sequences

Sequence divergence was estimated within noncoding sequences of both chloroplast DNA (cpDNA)trnL (UAA) intron and nuclear ribosomal DNA (nrDNA) internal transcribed spacer sequences (ITS1 and ITS2) for 10 species of the genusGentianaL. (Gentianaceae). Comparisons of evolutionary rates among these sequences (cpDNA versus nrDNA, ITS1 versus ITS2) were performed. It appears that sequence divergence is on average two to three times higher in ITSs than in thetrnL intron sequences and higher in ITS1 than in ITS2. Both the cpDNA intron and ITSs of nrDNA give concordant phylogenetic trees. However, the ITS-based phylogeny displays higher bootstrap values. At the intrageneric level, at least inGentiana,ITSs (especially ITS2) sequences seem to be more appropriate in the assessment of plant phylogenies. Nevertheless, the cpDNAtrnL intron seems to be preferable at the intergeneric level.     

Molecular Phylogeny of Magnolia (Magnoliaceae) Inferred from cpDNA Sequences and Evolutionary Divergence of the Floral Scents

Magnolia (disjuncts), however, have similar chemical profiles. A molecular phylogeny of Magnoliaceae was constructed to reveal phylogenetic relationships of taxa by sequencing the trnK intron (including the matK coding region), psbA-trnH, and atpB-rbcL intergenic spacer regions of chloroplast DNA from 25 Magnolia, two Michelia, and two Liriodendron taxa. The psbA-trnH spacer region showed twice the sequence divergence (0.0157) of the trnK intron (0.0073) or the matK coding region (0.0077). The strict consensus tree constructed from the combined data set (ca. 3,700 bp) indicated the genus Magnolia was polyphyletic containing Michelia species as ingroup. The clade of Magnolia liliifera var. obovata, M. coco, and M. delavayi formed the first branch. Among the remaining species, two additional large clades were recognized, i.e., one comprised of American evergreen Magnolia species and another of subgenus Yulania. The relationship among sect. Rytidospermum taxa was not clearly resolved. Parsimonious mapping of the floral scent chemical characters was performed onto the molecular phylogenetic tree to discuss evolutionary trends of the floral scent chemistries. Received 7 December 1998/ Accepted in revised form 18 May 1999     

Phylogeography of Larix sukaczewii Dyl. and Larix sibirica L. inferred from nucleotide variation of nuclear genes

We investigated phylogeography of Larix sukaczewii and Larix sibirica using nucleotide variation at three following nuclear gene regions: 5.8 S rDNA including two internal transcribed spacers (ITS), cinnamyl alcohol dehydrogenase (CAD), and phytochrome-O (PHYO). We also included sequences of the 4-coumarate: coenzyme A ligase (4CL) gene region obtained in our recent study. CAD and PHYO showed very low nucleotide variation, but ITS and 4CL had levels of variation similar to those reported for other conifers. Pleistocene refugia have been hypothesized to exist in the Southern Urals and South Central Siberia, where four out of nine of the investigated populations occur. We found moderate to high levels of population differentiation (F _ST = 0.115 – 0.531) in some pairwise comparisons suggesting limited gene flow and independent evolution of some refugial populations. In L. sukaczewii, low levels of differentiation were found among populations from areas glaciated during the Pleistocene, indicating their recent origin. Our results also suggest these populations were created by migrants from multiple, genetically distinct refugia. Furthermore, some haplotypes observed in populations from previously glaciated areas were not found in putative refugial populations, suggesting these populations might have contributed little to the extant populations created after the Last Glacial Maximum. Some authors regard L. sukaczewii and L. sibirica as a single species, while others consider them as separate species. The observed conspicuous differences in haplotype composition and distribution between L. sukaczewii and L. sibirica, together with high values of F _ST between populations of the two species, appear to support the latter classification. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. Ismael A. Khatab and Kariyawasam K.G.U. Hemamali contributed equally to this work.     

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     

DNA barcoding as a tool for early warning and monitoring alien duckweeds (Lemna sp.pl.): the case of Central Italy

Aquatic habitats are vulnerable to the invasion of alien species, so early warning protocols are necessary for eradication. The presence in Italy of two alien duckweeds in freshwaters has been documented: Lemna minuta, that showed high invasivity, and L. valdiviana, still confined to south Lazio. These two species may be mistaken for each other and for the domestic L. minor and L. gibba due to morphological variation. Here, we assess the applicability of DNA barcoding as a complement to morphological analysis for monitoring the spread of alien Lemna. We chose two chloroplast genome sequences for their ability to discriminate all Lemna species: the 5’ intron of the trnK gene and the matK gene. Among 48 samples of Lemna collected at 20 sites in Central Italy, 20 were identified as L. minor, 19 as L. minuta, five as L. trisulca and four as L. gibba. L. minuta was present at most sampling sites; in particular, at six locations of Lake Trasimeno, eight L. minuta samples were found. We demonstrate that DNA sequence analyses with cost-effective barcoding techniques can effectively support expert efforts in species determination for an early alert system of invasive Lemna species.     

A multi-locus plastid phylogenetic analysis of the pantropical genus Diospyros (Ebenaceae), with an emphasis on the radiation and biogeographic origins of the New Caledonian endemic species

We aimed to clarify phylogenetic relationships within the pantropical genus Diospyros (Ebenaceae sensu lato), and ascertain biogeographical patterns in the New Caledonian endemic species. We used DNA sequences from eight plastid regions (rbcL, atpB, matK, ndhF, trnK intron, trnL intron, trnL-trnF spacer, and trnS-trnG spacer) and included 149 accessions representing 119 Diospyros species in our analysis. Results from this study confirmed the monophyly of Diospyros with good support and provided a clearer picture of the relationships within the genus than in previous studies. Evidence from phylogenetic analyses suggests that Diospyros colonized New Caledonia multiple times. The four lineages of Diospyros in New Caledonia also differ in their degree of diversification. The molecular data indicate that one lineage is paleoendemic and derived from an ancient Australian species. The other three lineages are more closely related to several Southeast Asian species; two of them are neoendemics, and one has radiated rapidly and recently.     

Chloroplast DNA of black pine retains a residual inverted repeat lacking rRNA genes: nucleotide sequences of trnQ, trnK, psbA, trnI and trnH and the absence of rps16

Summary A physical map of black pine (Pinus thunbergii) chloroplast DNA (120 kb) was constructed and two separate portions of its nucleotide sequence were determined. One portion contains trnQ-UUG, ORF510, ORF83, trnK-UUU (ORF515 in the trnK intron), ORF22, psbA, trnI-CAU (on the opposing strand) and trnH-GUG, in that order. Sequence analysis of another portion revealed the presence of a 495 by inverted repeat containing trnI-CAU and the 3 end of psbA but lacking rRNA genes. The position of trnI-CAU is unique because most chloroplast DNAs have no gene between psbA and trnH (trnI-CAU is usually located further downstream). Black pine chloroplast DNA lacks rps16, which has been found between trnQ and trnK in angiosperm chloroplast DNAs, but possesses ORF510 instead. This ORF is highly homologous to ORF513 found in the corresponding region of liverwort chloroplast DNA and ORF563 located downstream from trnT in Chlamydomonas moewusii chloroplast DNA. A possible pathway for the evolution of black pine chloroplast DNA is discussed.     

Different evolutionary pathway of B*570101 and B*5801 (B17 group) alleles based in intron sequences

Two theories about MHC allele generation have been put forward: (1) point mutation diversification and/or (2) gene conversion events. A model supporting the existence of both of these mechanisms is shown in this paper; the possible evolution of the HLA-B*570101 and HLA-B*5801 alleles (which belong to the HLA-B17 serology group) is studied. The hypothesis favoured is that gene conversion events have originated these alleles, because intron sequences are also analysed. Evolution by point mutation should only be accepted if flanking introns have also been sequenced.The nucleotide sequence data (exons and introns) reported in this paper have been sequenced in our laboratory. They are in the GenBank nucleotide sequence database and have been assigned the accession numbers: B*150101—(a) exon 1, L79939; (b) exon 2 and exon 3, L48400; (c) intron 1, L76249; (d) intron 2, L42468; B*1515—(a) exon 1, exon 2 and exon 3, L49343; (b) intron 1 and intron 2, L76254; B*1539—(a) exon 2, AF033501; (b) exon 3, AF033502; (c) intron 1, AF034961; (d) intron 2 AF034962; B*350101—(a) exon 1, exon 2 and exon 3, L63544; (b) intron 1, L79921; (c) intron 2, L57505; B*510101—(a) exon 1, L77204; (b) exon 2 and exon 3, L47985; (c) intron 1, L76245; (d) intron 2, L42469; B*520102—(a) exon 1, L77205; (b) exon 2 and exon 3, L47984; (c) intron 1, L76244; (d) intron 2, L76251; B*5301—(a) exon 1, intron 1, exon 2, intron 2 and exon 3, U90566; B*1302—(a) intron 1, exon 2, intron 2, exon 3, AF196182; B*400101/02—(a) exon 2 and exon 3, L79937; (b) intron 1, L79919; (c) intron 2, L76629; B*4101—(a) intron 1, exon 2, intron 2 and exon 3, U90560; B*4102 (a) intron 1, exon 2, intron 2 and exon 3, AF 126199; B*4501—(a) intron 1, exon 2, intron 2 and exon 3, U90562; B*570101—(a) intron 1, exon 2, intron 2 and exon 3, AF196183; B*5801—(a) intron 1, exon 2, intron 2 and exon 3, AF196184All exon sequences were officially assigned as confirmatory by the WHO Nomenclature Committee in December 2003: B*1302, B*150101, B*350101, B*400101/02, B*4101, B*510101, B*570101, B*5801, B*5301, B*4501, B*520102, B*1515, B*4102 and B*1539. This follows the agreed policy that, subject to the conditions stated in the Nomenclature Report [Marsh et al. (2002) Tissue Antigens 60:407–464], names will be assigned to new sequences as they are identified. Lists of such new names will be published in the following WHO Nomenclature Report     

Patterns of DNA-Sequence Divergence Between <Emphasis Type="Italic">Drosophila miranda</Emphasis> and <Emphasis Type="Italic">D. pseudoobscura</Emphasis>

Contrary to the classical view, a large amount of non-coding DNA seems to be selectively constrained in Drosophila and other species. Here, using Drosophila miranda BAC sequences and the Drosophila pseudoobscura genome sequence, we aligned coding and non-coding sequences between D. pseudoobscura and D. miranda, and investigated their patterns of evolution. We found two patterns that have previously been observed in comparisons between Drosophila melanogaster and its relatives. First, there is a negative correlation between intron divergence and intron length, suggesting that longer non-coding sequences may contain more regulatory elements than shorter sequences. Our other main finding is a negative correlation between the rate of non-synonymous substitutions (d _N) and codon usage bias (F _op), showing that fast-evolving genes have a lower codon usage bias, consistent with strong positive selection interfering with weak selection for codon usage.     

Ancestral Population Sizes and Species Divergence Times in the Primate Lineage on the Basis of Intron and BAC End Sequences

The effective sizes of ancestral populations and species divergence times of six primate species (humans, chimpanzees, gorillas, orangutans, and representatives of Old World monkeys and New World monkeys) are estimated by applying the two-species maximum likelihood (ML) method to intron sequences of 20 different loci. Examination of rate heterogeneity of nucleotide substitutions and intragenic recombination identifies five outrageous loci (ODC1, GHR, HBE, INS, and HBG). The estimated ancestral polymorphism ranges from 0.21 to 0.96% at major divergences in primate evolution. One exceptionally low polymorphism occurs when African and Asian apes diverged. However, taking into consideration the possible short generation times in primate ancestors, it is concluded that the ancestral population size in the primate lineage was no smaller than that of extant humans. Furthermore, under the assumption of 6 million years (myr) divergence between humans and chimpanzees, the divergence time of humans from gorillas, orangutans, Old World monkeys, and New World monkeys is estimated as 7.2, 18, 34, and 65 myr ago, respectively, which are generally older than traditional estimates. Beside the intron sequences, three other data sets of orthologous sequences are used between the human and the chimpanzee comparison. The ML application to these data sets including 58,156 random BAC end sequences (BES) shows that the nucleotide substitution rate is as low as 0.6–0.8 × 10^–9 per site per year and the extent of ancestral polymorphism is 0.33–0.51%. With such a low substitution rate and short generation time, the relatively high extent of polymorphism suggests a fairly large effective population size in the ancestral lineage common to humans and chimpanzees.[Reviewing Editor: Dr. Magnus Nordborg]     

Evolutionary relationship of plant catalase genes inferred from exon-intron structures: isozyme divergence after the separation of monocots and dicots

In order to understand the molecular evolution of catalase genes in higher plants, we compared the exon-intron structures of 12 genomic sequences from six plant species. It was assumed that the putative single primordial catalase gene had seven introns, because only those catalase genes having this structure are found in the monocotyledonae and dicotyledonae classes. After the evolutionary divergence of monocots from dicots, consecutive duplication of the primordial gene followed by the differential loss of introns occurred in each class to form three (or possibly four in dicots) diverse isozyme genes. In monocots, three ancestral isozyme genes were formed before the divergence of ancestral rice and maize. One of the rice genes, CatA, has an entirely new short intron which was not found in any other plant catalase gene examined. We have investigated the existence of the intron in the CatA homolog in other rice species by polymerase chain reaction (PCR) analysis. One major PCR product was found with the genomic DNAs from O. sativa (indica and japonica types), O. rufipogon and O. glaberrima. DNAs from several accessions of O. longistaminata showed variation in both the number and size of the DNA fragments amplified. PCR analyses and sequencing of the PCR products revealed that there are several CatA homologs having different sequences in some accessions of O. longistaminata. We have extended our study to other species in the Poaceae. The results suggest that the gain of the intron, most likely by insertion of a retroposon, took place in the ancestral genome of rice after its evolutionary divergence from other ancestral cereals such as barley, wheat and oat. Received: 20 November 1997 / Accepted: 5 January 1998     

INTRON REVEALED BY NUCLEOTIDE SEQUENCE OF LARGE SUBUNIT OF RIBULOSE-1,5-BISPHOSPHATE CARBOXYLASE/OXYGENASE FROM CODIUM FRAGILE (CHLOROPHYTA): PHYLOGENETIC ANALYSIS1

The coding sequence for the large subunit of ribulose-1,5-bisphosphate carboxylase/oxygenase (rbcL) from Codium fragile (Suringar) Hariot chloroplast DNA is 1428 bp in length and contains a 1813-bp group II intron. The only other organisms in which introns have been found in the rbcL gene are Euglena and Astasia. The Codium intron likely had a separate origin from the Euglena and Astasia introns, based on comparisons of intron sizes and sequences. Phylogenetic analyses of rbcL nucleotide and amino acid sequences place Codium between Chlorella and two Chlamydomonas spp., indicating that the Chlorophyceae may be polyphyletic.     

Evidence for the Ancestral Origin of Group I Introns in the SSUrDNA of Naegleria spp.

ABSTRACT The sequence variation within the group I intron in five Naegleria spp. was studied and compared with the sequence variation within the flanking small subunit ribosomal DNA. Considerable sequence divergence was observed in the introns as well as in the rDNA. In the intron deletions and insertions are only detected in the sequence contributing to the secondary structure, not in the open reading frame. Most of the sequence variation is detected in the unpaired loops. In the case of nucleotide substitution in helices, compensating base pair changes were observed. The sequence variation does not induce variation in the secondary structure model. The phylogenetic tree based on the intron sequences is similar to the tree based on the flanking rDNA sequences. This observation indicates that the intron might have been acquired at an early stage in evolution, and lost in the majority of Naegleria spp.     

Molecular phylogenetics and taxonomic revision of Habenaria section Pentadactylae (Orchidaceae,Orchidinae)

As a step towards a revision of the sectional classification of Neotropical species of Habenaria, we focus here on section Pentadactylae. In its current delimitation, this is the largest of the 14 New World sections and embraces a group of 34 morphologically heterogeneous species. We expanded the sampling of Neotropical species currently placed in this section and performed Bayesian, maximum likelihood and parsimony analyses using nucleotide sequences from one nuclear (internal transcribed spacer, ITS) and three plastid (matK, trnK intron, rps16–trnK) DNA regions. In addition, morphological features of these species were reassessed. Based on our analyses, we propose that Habenaria section Pentadactylae should be recircumscribed to include only seven species: H. pentadactyla (the type species of the section), H. dutrae, H. ekmaniana, H. exaltata, H. henscheniana, H. megapotamensis and H. montevidensis. Thirty‐two species previously assigned to the section grouped within unrelated clades and are therefore excluded from the section. There are no unambiguous morphological synapomorphies for the section, but the group can be confidently recircumscribed and identified on the basis of a combination of diagnostic morphological vegetative and floral characters. Morphological floral features in Habenaria montevidensis are distinct from those of other species in the section, probably as a result of a shift to diurnal pollinators. Following a taxonomic revision of the group, H. crassipes is placed under the synonymy of H. exaltata and neotypes are designated for H. crassipes, H. montevidensis and H. recta (= H. ekmaniana). All species in the section live in marshes or wet grasslands from northern Argentina to central Brazil; most species are concentrated in southern Brazil. Most species are probably rare, and five may be threatened according to the World Conservation Union (IUCN) criteria. © 2014 The Linnean Society of London, Botanical Journal of the Linnean Society, 2014, 175 , 47–73.     

Genetic polymorphism in Caulerpa taxifolia (Ulvophyceae) chloroplast DNA revealed by a PCR‐based assay of the invasive Mediterranean strain

Abstract An invasive, cold‐tolerant strain of the tropical green alga Caulerpa taxifolia was introduced recently in the Mediterranean Sea and along the Californian coast. We screened 50 aquarium and open‐sea C. taxifolia specimens for the presence/absence of an intron located in the rbcL gene of chloroplast DNA. We also reanalysed a total of 229 sequences of the Internal Transcribed Spacer (ITS) of ribosomal DNA, combining previously published sequences from different studies with 68 new sequences to complement rbcL data. The introduced Mediterranean strain was found to be characterized by the absence of the rbcL intron and by the occurrence of a particular monomorphic ITS type. A PCR assay based on rbcL gene was developed to detect new introductions of the invasive strain of C. taxifolia. This rapid and inexpensive test could be useful to assist environment managers in the preservation of coastal marine ecosystems.     

Phylogenetics of early branching eudicots:Comparing phylogenetic signal across plastid introns,spacers,and genes

<正>Recent phylogenetic analyses revealed a grade with Ranunculales,Sabiales,Proteales,Trochodendrales, and Buxales as first branching eudicots,with the respective positions of Proteales and Sabiales still lacking statistical confidence.As previous analyses of conserved plastid genes remain inconclusive,we aimed to use and evaluate a representative set of plastid introns(group I:trnL;group II:petD,rpll6,trnK) and intergenic spacers(trnL-F,petB-petD, atpB-rbcL,rps3-rpl16) in comparison to the rapidly evolving matK and slowly evolving atpB and rbcL genes. Overall patterns of micro structural mutations converged across genomic regions,underscoring the existence of a general mutational pattern throughout the plastid genome.Phylogenetic signal differed strongly between functionally and structurally different genomic regions and was highest in matK,followed by spacers,then group II and group I introns.The more conserved atpB and rbcL coding regions showed distinctly lower phylogenetic information content.Parsimony,maximum likelihood,and Bayesian phylogenetic analyses based on the combined dataset of non-coding and rapidly evolving regions(14 000 aligned characters) converged to a backbone topology of eudicots with Ranunculales branching first,a Proteales-Sabiales clade second,followed by Trochodendrales and Buxales. Gunnerales generally appeared as sister to all remaining core eudicots with maximum support.Our results show that a small number of intron and spacer sequences allow similar insights into phylogenetic relationships of eudicots compared to datasets of many combined genes.The non-coding proportion of the plastid genome thus can be considered an important information source for plastid phylogenomics.     

Protein Polymorphism Is Negatively Correlated with Conservation of Intronic Sequences and Complexity of Expression Patterns in <Emphasis Type="Italic">Drosophila melanogaster</Emphasis>

We report a significant negative correlation between nonsynonymous polymorphism and intron length in Drosophila melanogaster. This correlation is similar to that between protein divergence and intron length previously reported in Drosophila. We show that the relationship can be explained by the content of conserved noncoding sequences (CNS) within introns. In addition, genes with a high regulatory complexity and many genetic interactions also exhibit larger amounts of CNS within their introns and lower values of nonsynonymous polymorphism. The present study provides relevant evidence on the importance of intron content and expression patterns on the levels of coding polymorphism. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users. [Reviewing Editor: Dr. Dmitri Petrov]