Molecular phylogenetics of moray eels (Muraenidae) demonstrates multiple origins of a shell-crushing jaw (Gymnomuraena, Echidna) and multiple colonizations of the Atlantic Ocean.

Moray eels (Muraenidae) are apex predators on coral reefs around the world, but they are not well studied because their cryptic habitats and occasionally aggressive behaviors make them difficult to collect. We provide a molecular phylogeny of moray eels including 44 species representing two subfamilies, eight genera, and all tropical ocean basins. Phylogenetic relationships among these taxa are estimated from portions of mitochondrial loci cytochrome b (632 bp) and cytochrome oxidase subunit 1 (596 bp), and portions of the nuclear loci RAG-1 (421 bp) and RAG-2 (754 bp). We test four sets of contrasting phylogenetic hypotheses using Bayes Factors, Shimodaira-Hasegawa tests, and Templeton tests. First, our results support the subfamily-level taxonomic distinction between true morays (Muraeninae) and snakemorays (Uropterygiinae), statistically rejecting hypotheses of non-monophyly for each subfamily. Second, we reject a monophyletic grouping of the genera Gymnomuraena and Echidna, which share a durophagous (shell-crushing) cranial morphology and dentition, indicating that the durophagous characters are not homologous. Third, we demonstrate that durophagous feeding habits and associated morphological characters have evolved in parallel in an ancestor of Gymnomuraena and at least three additional times within the genus Echidna. Finally, the tree topology indicates multiple invasions of the Atlantic from the Indo-Pacific, one of these occurring immediately prior to formation of the Isthmus of Panama approximately 2.8 MYA (million years ago) and one or two others occurring in the early to mid Miocene. Cladogenesis occurring within the Atlantic during the mid Miocene and Pliocene also contributed to moray species diversity. These data include a pair of sister species separated by the Isthmus of Panama, allowing a time-calibrated tree with an estimated crown age for Muraenidae at between 41 and 60 MYA, consistent with fossil evidence. Most lineage accumulation within morays occurred from the late Oligocene (24-27 MYA) through the Miocene (5-23 MYA) to the late Pliocene (∼ 2.5 MYA).

Genetic variation increases during biological invasion by a Cuban lizard.

A genetic paradox exists in invasion biology: how do introduced populations, whose genetic variation has probably been depleted by population bottlenecks, persist and adapt to new conditions? Lessons from conservation genetics show that reduced genetic variation due to genetic drift and founder effects limits the ability of a population to adapt, and small population size increases the risk of extinction. Nonetheless, many introduced species experiencing these same conditions during initial introductions persist, expand their ranges, evolve rapidly and become invasive. To address this issue, we studied the brown anole, a worldwide invasive lizard. Genetic analyses indicate that at least eight introductions have occurred in Florida from across this lizard's native range, blending genetic variation from different geographic source populations and producing populations that contain substantially more, not less, genetic variation than native populations. Moreover, recently introduced brown anole populations around the world originate from Florida, and some have maintained these elevated levels of genetic variation. Here we show that one key to invasion success may be the occurrence of multiple introductions that transform among-population variation in native ranges to within-population variation in introduced areas. Furthermore, these genetically variable populations may be particularly potent sources for introductions elsewhere. The growing problem of invasive species introductions brings considerable economic and biological costs. If these costs are to be mitigated, a greater understanding of the causes, progression and consequences of biological invasions is needed.

Niche lability in the evolution of a Caribbean lizard community.

Niche conservatism--the tendency for closely related species to be ecologically similar--is widespread. However, most studies compare closely related taxa that occur in allopatry; in sympatry, the stabilizing forces that promote niche conservatism, and thus inhibit niche shifts, may be countered by natural selection favouring ecological divergence to minimize the intensity of interspecific interactions. Consequently, the relative importance of niche conservatism versus niche divergence in determining community structure has received little attention. Here, we examine a tropical lizard community in which species have a long evolutionary history of ecological interaction. We find that evolutionary divergence overcomes niche conservatism: closely related species are no more ecologically similar than expected by random divergence and some distantly related species are ecologically similar, leading to a community in which the relationship between ecological similarity and phylogenetic relatedness is very weak. Despite this lack of niche conservatism, the ecological structuring of the community has a phylogenetic component: niche complementarity only occurs among distantly related species, which suggests that the strength of ecological interactions among species may be related to phylogeny, but it is not necessarily the most closely related species that interact most strongly.

Phylogeography of two moray eels indicates high dispersal throughout the indo-pacific.

Reef fishes disperse primarily as oceanic "pelagic" larvae, and debate continues over the extent of this dispersal, with recent evidence for geographically restricted (closed) populations in some species. In contrast, moray eels have the longest pelagic larval stages among reef fishes, possibly providing opportunities to disperse over great distances. We test this prediction by measuring mitochondrial DNA (mtDNA) and nuclear DNA variation in 2 species of moray eels, Gymnothorax undulatus (N = 165) and G. flavimarginatus (N = 124), sampled at 14-15 locations across the Indo-Pacific. The mtDNA data comprise 632 bp of cytochrome b and 596 bp of cytochrome oxidase I. Nuclear markers include 2 recombination-activating loci (421 bp of RAG-1 and 754 bp of RAG-2). Analyses of molecular variance and Mantel tests indicate little or no genetic differentiation, and no isolation by distance, across 22 000 km of the Indo-Pacific. We estimate that mitochondrial genetic variation coalesces within the past about 2.3 million years (My) for G. flavimarginatus and within the past about 5.9 My for G. undulatus. Permutation tests of geographic distance on the mitochondrial haplotype networks indicate recent range expansions for some younger haplotypes (estimated within approximately 600 000 years) and episodic fragmentation of populations at times of low sea level. Our results support the predictions that the extended larval durations of moray eels enable ocean-wide genetic continuity of populations. This is the first phylogeographic survey of the moray eels, and morays are the first reef fishes known to be genetically homogeneous across the entire Indo-Pacific.

The role of geography and ecological opportunity in the diversification of day geckos (Phelsuma).

We examine the effects of ecological opportunity and geographic area on rates of species accumulation and morphological evolution following archipelago colonization in day geckos (genus Phelsuma) in the Indian Ocean. Using a newly generated molecular phylogeny for the genus, we present evidence that these geckos likely originated on Madagascar, whereas colonization of three archipelagos in the Indian Ocean, the Seychelles, Mascarene, and Comoros Islands has produced three independent monophyletic radiations. We find that rates of species accumulation are not elevated following colonization but are roughly equivalent on all three isolated archipelagos and on the larger island of Madagascar. However, rates of species accumulation have slowed through time on Madagascar. Rates of morphological evolution are higher in both the Mascarene and Seychelles archipelagos compared to rates on Madagascar. This negative relationship between rate of morphological evolution and island area suggests that ecological opportunity is an important factor in diversification of day gecko species.

In situ genetic differentiation in a Hispaniolan lizard (Ameiva chrysolaema): a multilocus perspective.

A previous phylogeographic study of mitochondrial haplotypes for the Hispaniolan lizard Ameiva chrysolaema revealed deep genetic structure associated with seawater inundation during the late Pliocene/early Pleistocene and evidence of subsequent population expansion into formerly inundated areas. We revisit hypotheses generated by our previous study using increased geographic sampling of populations and analysis of three nuclear markers (alpha-enolase intron 8, alpha-cardiac-actin intron 4, and beta-actin intron 3) in addition to mitochondrial haplotypes (ND2). Large genetic discontinuities correspond spatially and temporally with historical barriers to gene flow (sea inundations). NCPA cross-validation analysis and Bayesian multilocus analyses of divergence times (IMa and MCMCcoal) reveal two separate episodes of fragmentation associated with Pliocene and Pleistocene sea inundations, separating the species into historically separate Northern, East-Central, West-Central, and Southern population lineages. Multilocus Bayesian analysis using IMa indicates asymmetrical migration from the East-Central to the West-Central populations following secondary contact, consistent with expectations from the more pervasive sea inundation in the western region. The West-Central lineage has a genetic signature of population growth consistent with the expectation of geographic expansion into formerly inundated areas. Within each lineage, significant spatial genetic structure indicates isolation by distance at comparable temporal scales. This study adds to the growing body of evidence that vicariant speciation may be the prevailing source of lineage accumulation on oceanic islands. Thus, prior theories of island biogeography generally underestimate the role and temporal scale of intra-island vicariant processes.

Multiple sources, admixture, and genetic variation in introduced anolis lizard populations.

Invasive species are classically thought to suffer from reduced within-population genetic variation compared to their native-range sources due to founder effects and population bottlenecks during introduction. Reduction in genetic variation in introduced species may limit population growth, increase the risk of extinction, and constrain adaptation, hindering the successful establishment and spread of an alien species. Results of recent empirical studies, however, show higher than expected genetic variation, rapid evolution, and multiple native-range sources in introduced populations, which challenge the classical scenario of invasive-species genetics. With mitochondrial DNA (mtDNA) sequence data, we examined the molecular genetics of 10 replicate introductions of 8 species of Anolis lizards. Eighty percent of introductions to Florida and the Dominican Republic were from multiple native-range source populations. MtDNA haplotypes restricted to different geographically distinct populations in the native range of a species commonly occurred as intrapopulation polymorphisms in introduced populations. Two-thirds of introduced populations had two or more sources, and admixture elevated genetic variation in half of the introduced populations above levels typical of native-range populations. The mean pairwise sequence divergence among haplotypes sampled within introduced populations was nearly twice that within native-range populations (2.6% vs. 1.4%). The dynamics of introductions from multiple sources and admixture explained the observed genetic contrasts between native and introduced Anolis populations better than the classical scenario for most introduced populations. Elevated genetic variation through admixture occurred regardless of the mode or circumstances of an introduction. Little insight into the number of sources or amount of genetic variation in introduced populations was gained by knowing the number of physical introductions, the size of a species' non-native range, or whether it was a deliberate or accidental introduction. We hypothesize that elevated genetic variation through admixture of multiple sources is more common in biological invasions than previously thought. We propose that introductions follow a sequential, two-step process involving a reduction in genetic variation due to founder effects and population bottlenecks followed by an increase in genetic variation if admixture of individuals from multiple native-range sources occurs.

Rapid lineage accumulation in a non-adaptive radiation: phylogenetic analysis of diversification rates in eastern North American woodland salamanders (Plethodontidae: Plethodon).

Adaptive radiations have served as model systems for quantifying the build-up of species richness. Few studies have quantified the tempo of diversification in species-rich clades that contain negligible adaptive disparity, making the macroevolutionary consequences of different modes of evolutionary radiation difficult to assess. We use mitochondrial-DNA sequence data and recently developed phylogenetic methodologies to explore the tempo of diversification of eastern North American Plethodon, a species-rich clade of woodland salamanders exhibiting only limited phenotypic disparity. Lineage-through-time analysis reveals a high rate of lineage accumulation, 0.8 species per million years, occurring 11-8 million years ago in the P. glutinosus species group, followed by decreasing rates. This high rate of lineage accumulation is exceptional, comparable to the most rapid of adaptive radiations. In contrast to classic models of adaptive radiation where ecological niche divergence is linked to the origin of species, we propose that phylogenetic niche conservatism contributes to the rapid accumulation of P. glutinosus-group lineages by promoting vicariant isolation and multiplication of species across a spatially and temporally fluctuating environment. These closely related and ecologically similar lineages persist through long-periods of evolutionary time and form strong barriers to the geographic spread of their neighbours, producing a subsequent decline in lineage accumulation. Rapid diversification among lineages exhibiting long-term maintenance of their bioclimatic niche requirements is an under-appreciated phenomenon driving the build-up of species richness.

Phylogenetic analysis of ecomorphological divergence, community structure, and diversification rates in dusky salamanders (Plethodontidae: Desmognathus).

An important dimension of adaptive radiation is the degree to which diversification rates fluctuate or remain constant through time. Focusing on plethodontid salamanders of the genus Desmognathus, we present a novel synthetic analysis of phylogeographic history, rates of ecomorphological evolution and species accumulation, and community assembly in an adaptive radiation. Dusky salamanders are highly variable in life history, body size, and ecology, with many endemic lineages in the southern Appalachian Highlands of eastern North America. Our results show that life-history evolution had important consequences for the buildup of plethodontid-salamander species richness and phenotypic disparity in eastern North America, a global hot spot of salamander biodiversity. The origin of Desmognathus species with aquatic larvae was followed by a high rate of lineage accumulation, which then gradually decreased toward the present time. The peak period of lineage accumulation in the group coincides with evolutionary partitioning of lineages with aquatic larvae into seepage, stream-edge, and stream microhabitats. Phylogenetic simulations demonstrate a strong correlation between morphology and microhabitat ecology independent of phylogenetic effects and suggest that ecomorphological changes are concentrated early in the radiation of Desmognathus. Deep phylogeographic fragmentation within many codistributed ecomorph clades suggests long-term persistence of ecomorphological features and stability of endemic lineages and communities through multiple climatic cycles. Phylogenetic analyses of community structure show that ecomorphological divergence promotes the coexistence of lineages and that repeated, independent evolution of microhabitat-associated ecomorphs has a limited role in the evolutionary assembly of Desmognathus communities. Comparing and contrasting our results to other adaptive radiations having different biogeographic histories, our results suggest that rates of diversification during adaptive radiation are intimately linked to the degree to which community structure persists over evolutionary time.

Phylogenetic analysis of ecological and morphological diversification in Hispaniolan trunk-ground anoles (Anolis cybotes group).

Anolis lizards in the Greater Antilles partition the structural microhabitats available at a given site into four to six distinct categories. Most microhabitat specialists, or ecomorphs, have evolved only once on each island, yet closely related species of the same ecomorph occur in different geographic macrohabitats across the island. The extent to which closely related species of the same ecomorph have diverged to adapt to different geographic macrohabitats is largely undocumented. On the island of Hispaniola, members of the Anolis cybotes species group belong to the trunk-ground ecomorph category. Despite evolutionary stability of their trunk-ground microhabitat, populations of the A. cybotes group have undergone an evolutionary radiation associated with geographically distinct macrohabitats. A combined phylogeographic and morphometric study of this group reveals a strong association between macrohabitat type and morphology independent of phylogeny. This association results from long-term morphological evolutionary stasis in populations associated with mesic-forest environments (A. c. cybotes and A. marcanoi) and predictable morphometric changes associated with entry into new macrohabitat types (i.e., xeric forests, high-altitude pine forest, rock outcrops). Phylogeographic analysis of 73 new mitochondrial DNA sequences (1921 aligned sites) sampled from 68 geographic populations representing 12 recognized species and subspecies diagnoses 16 allopatric or parapatric groupings of populations differing from each other by 5-18% sequence divergence. At least some of these groupings appear to have attained species-level divergence from others. Evolutionary specialization to different macrohabitat types may be a major factor in the evolutionary diversification of Greater Antillean anoles.

Molecular phylogenetic evidence for ancient divergence of lizard taxa on either side of Wallace's Line.

Wallace's Line, separating the terrestrial faunas of South East Asia from the Australia-New Guinea region, is the most prominent and well-studied biogeographical division in the world. Phylogenetically distinct subgroups of major animal and plant groups have been documented on either side of Wallace's Line since it was first proposed in 1859. Despite its importance, the temporal history of fragmentation across this line is virtually unknown and the geological foundation has rarely been discussed. Using molecular phylogenetics and dating techniques, we show that the split between taxa in the South East Asian and the Australian-New Guinean geological regions occurred during the Late Jurassic to Early Cretaceous in two independent lizard clades. This estimate is compatible with the hypothesis of rifting Gondwanan continental fragments during the Mesozoic and strongly rejects the hypothetical origin of various members of the Australian-New Guinean herpetofauna as relatively recent invasions from South East Asia. Our finding suggests an ancient fragmentation of lizard taxa on either side of Wallace's Line and provides further evidence that the composition of modern global communities has been significantly affected by rifting and accretion of Gondwanan continental plates during the Middle to Late Mesozoic.

Partial island submergence and speciation in an adaptive radiation: a multilocus analysis of the Cuban green anoles.

Sympatric speciation is often proposed to account for species-rich adaptive radiations within lakes or islands, where barriers to gene flow or dispersal may be lacking. However, allopatric speciation may also occur in such situations, especially when ranges are fragmented by fluctuating water levels. We test the hypothesis that Miocene fragmentation of Cuba into three palaeo-archipelagos accompanied species-level divergence in the adaptive radiation of West Indian Anolis lizards. Analysis of morphology, mitochondrial DNA (mt DNA) and nuclear DNA in the Cuban green anoles (carolinensis subgroup) strongly supports three pre dictions made by this hypothesis. First, three geographical sets of populations, whose ranges correspond with palaeo-archipelago boundaries, are distinct and warrant recognition as independent evolutionary lineages or species. Coalescence of nuclear sequence fragments sampled from these species and the large divergences observed between their mtDNA haplotypes suggest separation prior to the subsequent unification of Cuba ca. 5 Myr ago. Second, molecular phylogenetic relationships among these species reflect historical geographical relationships rather than morphological similarity. Third, all three species remain distinct despite extensive geographical contact subsequent to island unification, occasional hybridization and introgression of mtDNA haplotypes. Allopatric speciation initiated during partial island submergence may play an important role in speciation during the adaptive radiation of Anolis lizards.

Admixture determines genetic diversity and population differentiation in the biological invasion of a lizard species.

Molecular genetic analyses show that introduced populations undergoing biological invasions often bring together individuals from genetically disparate native-range source populations, which can elevate genotypic variation if these individuals interbreed. Differential admixture among multiple native-range sources explains mitochondrial haplotypic diversity within and differentiation among invasive populations of the lizard Anolis sagrei. Our examination of microsatellite variation supports the hypothesis that lizards from disparate native-range sources, identified using mtDNA haplotypes, form genetically admixed introduced populations. Furthermore, within-population genotypic diversity increases with the number of sources and among-population genotypic differentiation reflects disparity in their native-range sources. If adaptive genetic variation is similarly restructured, then the ability of invasive species to adapt to new conditions may be enhanced.

Differential admixture shapes morphological variation among invasive populations of the lizard Anolis sagrei.

The biological invasion of the lizard Anolis sagrei provides an opportunity to study evolutionary mechanisms that produce morphological differentiation among non-native populations. Because the A. sagrei invasion represents multiple native-range source populations, differential admixture as well as random genetic drift and natural selection, could shape morphological evolution during the invasion. Mitochondrial DNA (mtDNA) analyses reveal seven distinct native-range source populations for 10 introduced A. sagrei populations from Florida, Louisiana and Texas (USA), and Grand Cayman, with 2-5 native-range sources contributing to each non-native population. These introduced populations differ significantly in frequencies of haplotypes from different native-range sources and in body size, toepad-lamella number, and body shape. Variation among introduced populations for both lamella number and body shape is explained by differential admixture of various source populations; mean morphological values of introduced populations are correlated with the relative genetic contributions from different native-range source populations. The number of source populations contributing to an introduced population correlates with body size, which appears independent of the relative contributions of particular source populations. Thus, differential admixture of various native-range source populations explains morphological differences among introduced A. sagrei populations. Morphological differentiation among populations is compatible with the hypothesis of selective neutrality, although we are unable to test the hypothesis of interdemic selection among introductions from different native-range source populations.

Molecular phylogenetics of the Anolis onca series: a case history in retrograde evolution revisited.

Anoles of the Anolis onca series represent a dramatic case of retrograde evolution, exhibiting great reduction (A. annectens) and loss (A. onca) of the subdigital pads considered a key innovation for the evolutionary radiation of anoles in arboreal environments. We present a molecular phylogenetic analysis of these anoles and their closest known relatives (A. auratus, A. lineatus, A. meridionalis, and A. nitens) using new mitochondrial DNA sequence data from the ND2 gene, five tRNA genes (tRNA(Trp), tRNA(Ala), tRNA(Asn), tRNA(Cys), tRNA(Tyr)), the origin of light-strand replication, and a portion of the CO1 gene (1,446 aligned base positions, 612 parsimony informative). Our results confirm monophyly of the A. onca series and suggest an evolutionary separation of approximately 10 million years between A. annectens and A. onca. Evolution of subdigital structure in this series illustrates ectopic expression of developmental programs that replace flexible subdigital lamellae of the toepad with rigid, keeled scales resembling dorsal digital scales. Our phylogenetic results indicate that narrowing of the toepad in A. auratus evolved separately from toepad reduction in the A. onca series. Expansion of the subdigital lamellae along the phalanges in A. auratus appears to compensate constriction of lamellae by digital narrowing, maintaining greater climbing capability in this species. Toepad evolution in the lineage ancestral to A. auratus features changes of the same developmental modules as the A. onca series but in the opposite direction. Large molecular distances between geographic populations of A. auratus indicate that its derived toepad structure is at least 9 million years old.

Gene lineages and eastern North American palaeodrainage basins: phylogeography and speciation in salamanders of the Eurycea bislineata species complex.

Contemporary North American drainage basins are composites of formerly isolated drainages, suggesting that fragmentation and fusion of palaeodrainage systems may have been an important factor generating current patterns of genetic and species diversity in stream-associated organisms. Here, we combine traditional molecular-phylogenetic, multiple-regression, nested clade, and molecular-demographic analyses to investigate the relationship between phylogeographic variation and the hydrogeological history of eastern North American drainage basins in semiaquatic plethodontid salamanders of the Eurycea bislineata species complex. Four hundred forty-two sequences representing 1108 aligned bases from the mitochondrial genome are reported for the five formally recognized species of the E. bislineata complex and three outgroup taxa. Within the in-group, 270 haplotypes are recovered from 144 sampling locations. Geographic patterns of mtDNA-haplotype coalescence identify 13 putatively independent population-level lineages, suggesting that the current taxonomy of the group underestimates species-level diversity. Spatial and temporal patterns of phylogeographic divergence are strongly associated with historical rather than modern drainage connections, indicating that shifts in major drainage patterns played a pivotal role in the allopatric fragmentation of populations and build-up of lineage diversity in these stream-associated salamanders. More generally, our molecular genetic results corroborate geological and faunistic evidence suggesting that palaeodrainage connections altered by glacial advances and headwater erosion occurring between the mid-Miocene and Pleistocene epochs explain regional patterns of biodiversity in eastern North American streams.

Assembly of the eastern North American herpetofauna: new evidence from lizards and frogs.

Darwin first recognized the importance of episodic intercontinental dispersal in the establishment of worldwide biotic diversity. Faunal exchange across the Bering Land Bridge is a major example of such dispersal. Here, we demonstrate with mitochondrial DNA evidence that three independent dispersal events from Asia to North America are the source for almost all lizard taxa found in continental eastern North America. Two other dispersal events across Beringia account for observed diversity among North American ranid frogs, one of the most species-rich groups of frogs in eastern North America. The contribution of faunal elements from Asia via dispersal across Beringia is a dominant theme in the historical assembly of the eastern North American herpetofauna.

A molecular assessment of phylogenetic relationships and lineage accumulation rates within the family Salamandridae (Amphibia, Caudata).

We examine phylogenetic relationships among salamanders of the family Salamandridae using approximately 2700 bases of new mtDNA sequence data (the tRNALeu, ND1, tRNAIle, tRNAGln, tRNAMet, ND2, tRNATrp, tRNAAla, tRNAAsn, tRNACys, tRNATyr, and COI genes and the origin for light-strand replication) collected from 96 individuals representing 61 of the 66 recognized salamandrid species and outgroups. Phylogenetic analyses using maximum parsimony and Bayesian analysis are performed on the new data alone and combined with previously reported sequences from other parts of the mitochondrial genome. The basal phylogenetic split is a polytomy of lineages ancestral to (1) the Italian newt Salamandrina terdigitata, (2) a strongly supported clade comprising the "true" salamanders (genera Chioglossa, Mertensiella, Lyciasalamandra, and Salamandra), and (3) a strongly supported clade comprising all newts except S. terdigitata. Strongly supported clades within the true salamanders include monophyly of each genus and grouping Chioglossa and Mertensiella as the sister taxon to a clade comprising Lyciasalamandra and Salamandra. Among newts, genera Echinotriton, Pleurodeles, and Tylototriton form a strongly supported clade whose sister taxon comprises the genera Calotriton, Cynops, Euproctus, Neurergus, Notophthalmus, Pachytriton, Paramesotriton, Taricha, and Triturus. Our results strongly support monophyly of all polytypic newt genera except Paramesotriton and Triturus, which appear paraphyletic, and Calotriton, for which only one of the two species is sampled. Other well-supported clades within newts include (1) Asian genera Cynops, Pachytriton, and Paramesotriton, (2) North American genera Notophthalmus and Taricha, (3) the Triturus vulgaris species group, and (4) the Triturus cristatus species group; some additional groupings appear strong in Bayesian but not parsimony analyses. Rates of lineage accumulation through time are evaluated using this nearly comprehensive sampling of salamandrid species-level lineages. Rate of lineage accumulation appears constant throughout salamandrid evolutionary history with no obvious fluctuations associated with origins of morphological or ecological novelties.

Codon usage patterns in Nematoda: analysis based on over 25 million codons in thirty-two species.

BACKGROUND: Codon usage has direct utility in molecular characterization of species and is also a arker for molecular evolution. To understand codon usage within the diverse phylum Nematoda,we analyzed a total of 265,494 expressed sequence tags (ESTs) from 30 nematode species. The full genomes of Caenorhabditis elegans and C. briggsae were also examined. A total of 25,871,325 codons ere analyzed and a comprehensive codon usage table for all species was generated. This is the first codon usage table available for 24 of these organisms. RESULTS: Codon usage similarity in Nematoda usually persists over the breadth of a genus but thenrapidly diminishes even within each clade. Globodera, Meloidogyne, Pristionchus, and Strongyloides have the most highly derived patterns of codon usage. The major factor affecting differences in codon usage between species is the coding sequence GC content, which varies in nematodes from 32%to 51%. Coding GC content (measured as GC3) also explains much of the observed variation in the effective number of codons (R = 0.70), which is a measure of codon bias, and it even accounts for differences in amino acid frequency. Codon usage is also affected by neighboring nucleotides(N1 context). Coding GC content correlates strongly with estimated noncoding genomic GC content (R = 0.92). On examining abundant clusters in five species, candidate optimal codons were identified that may be preferred in highly expressed transcripts. CONCLUSION: Evolutionary models indicate that total genomic GC content, probably the product of directional mutation pressure, drives codon usage rather than the converse, a conclusion that is supported by examination of nematode genomes.