Fixed allele frequency differences among Palauan and Okinawan populations of the damselfishes Chrysiptera cyanea and Pomacentrus coelestis

A survey of intraspecific allozymic variation among samples of the damselfishes Chrysiptera cyanea and Pomacentrus coelestis collected from Palau and Okinawa revealed levels of genetic divergence far in excess of estimates reported previously for populations of coral reef fishes. Absolute or nearly fixed differences in allele frequencies were detected at ADA ^*, sMDH-2 ^*, MEP-1 ^*, PEPB ^*, PEP-LT ^* and sSOD ^*, and at sAAT ^*, GPI-A ^*, LDH-1 ^*, and PEPB ^*, among the geographic samples of C. cyanea and P. coelestis, respectively. Examination of allele frequencies at most other loci (26 for C. cyanea; 24 for P. coelestis) revealed slight differentiation or identical fixation among the geographic samples. The observed patterns of allele frequency distribution suggest the existence of localized demes of these fishes: these demes may be cryptic and/or incipient species. Whether or not speciation is incipient, the observed patterns of allozymic and isozymic variability suggest that natural selection is a factor in the maintenance of population substructuring of the study species. Pronounced allelic differences were highly concentrated at a small number of loci: strongly bimodal frequency distributions of loci relative to their genetic identities (I) were observed for among-population comparisons in both study species. Allozymes encoded by the PEPB ^* locus in samples of the noncongeners from Palau exhibited identical electrophoretic mobilities. In Palauan P. coelestis, mSOD ^* is expressed but sSOD ^* apparently is not, whereas in Okinawan P. coelestis, both mSOD ^* and sSOD ^* are expressed.     

Incipient speciation within a subgenus of rockfish (<Emphasis Type="Italic">Sebastosomus</Emphasis>) provides evidence of recent radiations within an ancient species flock

The high frequency of speciation events associated with species flocks (i.e., radiations of closely related species) provides invaluable insight into the speciation process. Investigations of the speciation process in the marine environment are rare, and therefore, the genetic analysis of the rockfish genus Sebastes, considered an ancient marine species flock, provides an opportunity to investigate this process in the sea. Using both mitochondrial and nuclear markers, we analyzed five closely related species within the rockfish subgenus Sebastosomus. Our goal was to understand the evolutionary history and genetic relationships among species within this group and to provide evidence of recent speciation events within the subgenus. In the genetic analysis of the subgenus, we found different stages of the speciation process, with greater genetic divergences among three of the five species, evidence of recent divergence between two of the five species, Sebastes entomelas and S. mystinus, and significant genetic divergence between two lineages within S. mystinus revealing a signature of incipient speciation. We also found frequency differences of the two S. mystinus lineages among sample locations and found no evidence of introgression between the lineages at the location where both coexist. Although Sebastes is an example of an ancient species flock, this study provides evidence of ongoing speciation within the genus and reveals stages of this process from incipient to distinct species.     

Divergence of mitochondrial DNA lineage of the rocky intertidal goby Chaenogobius gulosus around the Japanese Archipelago: reference to multiple Pleistocene isolation events in the Sea of Japan

Geographic variations in the mitochondrial cytochrome b gene sequence (1,062 nucleotides) of the rocky intertidal goby Chaenogobius gulosus were examined and compared against the phylogeographic pattern of its sympatric sister species Chaenogobius annularis, which we studied previously. We collected C. gulosus specimens from 17 localities around the Japanese Archipelago and Korean Peninsula. Thirty-one haplotypes were obtained. The phylogenetic tree showed three geographical lineages, of which one was distributed along the coast of the Sea of Japan, suggesting that a Pleistocene isolation event in the Sea of Japan affected the genetic divergence of this species, similar to that which was observed in C. annularis. The net nucleotide divergence between the Sea of Japan and Pacific Ocean lineages in C. gulosus was much lower than that of C. annularis. Both traditional molecular clock and hierarchical approximate Bayesian computational models suggest that different divergence events affected the genetic divergences of C. gulosus and C. annularis. Our study on the sister species pair strongly supports multiple isolation events in the Sea of Japan, leading to varying levels of genetic divergence between the Pacific Ocean and Sea of Japan lineages in the Japanese coastal marine species.     

Comparison of larval duration and pre- and post-settlement growth in two species of damselfish,Chromis atripectoralis andPomacentrus coelestis (Pisces: Pomacentridae), from the Great Barrier Reef

Light traps were used to capture larval fishes, immediately before settlement, at two localities 500 km apart on the Great Barrier Reef (GBR) in December, 1987. Samples from Lizard Island, in the northern GBR, and Davies Reef, in the central GBR, were dominated by two species of damselfish:Chromis atripectralis andPomacentrus coelestis. Analysis of otoliths revealed significant differences in both size and age at settlement between the two localities forP. coelestis, but not forC. atripectoralis. Growth rates determined for pre- and post-settlementP. coelestis suggested a sigmoidal growth trajectory through the larval life, with growth slowing as fishes approached the time of settlement. Post-settlement growth rates were faster than growth prior to settlement in both species. Growth in both species was, however, similar between localities. The relationship between fish size and otolith size was complex, varying both between pre- and post-settlement fishes, and among localities. This emphasizes the need to validate the relationship between fish size and otolith size before otoliths may be used to back-calculate individual growth trajectories.Contribution No. 500 from the Australian Institute of Marine Science     

Pleistocene speciation of sister taxa in a North Pacific clade of brooding sea stars (<Emphasis Type="Italic">Leptasterias</Emphasis>)

Although numerous coastal marine species show intra-specific lineage splitting and population divergence times that date to the period of glacial cycles during the Pleistocene epoch, reported instances of recent speciation in the coastal marine environment are relatively rare. Marine organisms with brood-protection and other reproductive modes with limited dispersal potential have been suggested to experience more frequent speciation and extinction events than related species with higher dispersal rates, but few studies have actually estimated divergence times of sister species in these organisms. Here, two mitochondrial gene regions (cytochrome oxidase subunit I, putative control region and upstream tRNAs) and a nuclear gene region (Elongation factor 1α subunit intron 4) provide evidence of recent (0.5–1.2 Mya) cladogenetic events in four pairs of putative sister taxa in a predominantly North Pacific brooding subgenus of sea stars (Leptasterias subgenus Hexasterias). Calibration is obtained from a trans-arctic migration in a related clade of sea stars (Leptasterias subgenera Hexasterias and Nesasterias) that is timed to the opening of the Bering Strait at 3.5 ± 0.25 Mya, and uncertainty in the calibration point is accommodated with a normally-distributed Bayesian prior probability. Similar estimates of population splitting times for two of the pairs of putative sister taxa were obtained by a multilocus coalescent analysis. Estimates of mitochondrial mutation rates (0.01/My) were approximately 50% of the values calibrated for sister species pairs in tropical sea stars and sea urchins. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Phylogeography and Limited Genetic Connectivity in the Endangered Boring Giant Clam across the Coral Triangle

Abstract: The Coral Triangle is the global center of marine biodiversity; however, its coral reefs are critically threatened. Because of the bipartite life history of many marine species with sedentary adults and dispersive pelagic larvae, designing effective marine protected areas requires an understanding of patterns of larval dispersal and connectivity among geographically discrete populations. We used mtDNA sequence data to examine patterns of genetic connectivity in the boring giant clam (Tridacna crocea) in an effort to guide conservation efforts within the Coral Triangle. We collected an approximately 485 base pair fragment of mtDNA cytochrome c oxidase 1 (CO1) from 414 individuals at 26 sites across Indonesia. Genetic structure was strong between regions (φ_ST=0.549, p < 0.00001) with 3 strongly supported clades: one restricted to western Sumatra, another distributed across central Indonesia, and a third limited to eastern Indonesia and Papua. Even within the single largest clade, small but significant genetic structure was documented (φ_ST=0.069, p < 0.00001), which indicates limited gene flow within and among phylogeographic regions. Significant patterns of isolation by distance indicated an average dispersal distance of only 25–50 km, which is far below dispersal predictions of 406–708 km derived from estimates of passive dispersal over 10 days via surface currents. The strong regional genetic structure we found indicates potent limits to genetic and demographic connectivity for this species throughout the Coral Triangle and provides a regional context for conservation planning. The recovery of 3 distinct evolutionarily significant units within a well‐studied taxonomic group suggests that biodiversity in this region may be significantly underestimated and that Tridacna taxa may be more endangered than currently recognized.     

Ecological niche similarities of Peponapis bees and non-domesticated Cucurbita species

Peponapis bees are considered specialized pollinators of Cucurbita flowers, a genus that presents several species of economic value (squashes and pumpkins). Both genera originated in the Americas, and their diversity dispersion center is in Mexico. Ten species of Peponapis and ten species of Cucurbita (only non-domesticated species) were analyzed considering the similarity of their ecological niche characteristics with respect to climatic conditions of their occurrence areas (abiotic variables) and interactions between species (biotic variables). The similarity of climatic conditions (temperature and precipitation) was estimated through cluster analyses. The areas of potential occurrence of the most similar species were obtained through ecological niche modeling and summed with geographic information system tools. Three main clusters were obtained: one with species that shared potential occurrence areas mainly in deserts (P. pruinosa, P. timberlakei, C. digitata, C. palmata, C. foetidissima), another in moist forests (P. limitaris, P. atrata, C. lundelliana, C. o. martinezii) and a third mainly in dry forests (C. a. sororia, C. radicans, C. pedatifolia, P. azteca, P. smithi, P. crassidentata, P. utahensis). Some species with similar ecological niche presented potential shared areas that are also similar to their geographical distribution, like those occurring predominantly on deserts. However, some clustered species presented larger geographical areas, such as P. pruinosa and C. foetidissima suggesting other drivers than climatic conditions to shape their distributions. The domestication of Cucurbita and also the natural history of both genera were considered also as important factors.     

Resource-associated divergence in the arctic marine amphipod<Emphasis Type="Italic"> Paramphithoe hystrix</Emphasis>

Host specialization has played an important role in the speciation of groups such as herbivorous insects and parasitic invertebrates. In this study, we provide evidence for its role in the diversification of the amphipod Paramphithoe hystrix Ross 1835 in Arctic marine waters. This species lives and feeds on the tissues of varied invertebrate hosts that are much larger than itself, including poriferans, hydrozoans and echinoderms, acting as a 'micropredator'. We examined the genetic structure of P. hystrix at sites in the Canadian Arctic where it is represented by a white form found on the light-coloured sea star Solaster endeca, a red form on the pink soft coral Gersemia rubriformis and by a spotted form with an unidentified host. These phenotypes occur in microallopatry, sometimes occupying alternate hosts just a few metres apart. Although their variation in host associates might have arisen from ingested host tissue, our results indicate that these morphs are genetically distinct. Mitochondrial (cytochrome c oxidase I) and nuclear (28S ribosomal RNA) gene sequences, as well as allozyme data all show strong genetic divergence between the three colour morphs, demonstrating their reproductive isolation. These results suggest that P. hystrix is a complex of at least three species, which may have arisen as a result of disruptive selection following host switches. Moreover, the depths of genetic divergence indicate that diversification of this complex was complete prior to the Pleistocene.Communicated by R.J. Thompson, St. John's     

Relation between respiration,excretion (ammonia and inorganic phosphorus) and activity of amylase and trypsin in different species of pelagic copepods from an Indian Ocean equatorial area

The relationship between respiration, nitrogen and phosphorus excretion and specific activity of amylase and trypsin was investigated from shipboard experiments with several species of copepods from the equatorial divergence area of the Indian Ocean, in June 1978. Statistical analysis of O:N ratios in 6 of the most common species revealed three main groups: a group displaying low O:N (Pontella fera, Candacia pachydactyla); a group with high O:N (Undinula darwini, Euchaeta marina, Temora discaudata); and a species displaying an intermediate O:N ratio: Scolecithrix bradyi. O:P and N:P did not differ significantly between species. There was a direct relationship between average O:N and the ratios of specific activity of the digestive enzymes amylase and trypsin. Species displaying low O:N and/or A:T ratios such as P. fera, C. pachydactyla and Oncaea venusta probably metabolize proteins more efficiently than they do plant carbohydrates (high nitrogen excretion and low amylase activity). Species displaying high O:N and A:T ratios, such as U. darwini, E. marina and T. dicaudata (low nitrogen excretion and high amylase activity) were assumed to use carbohydrates (starch) and proteins with equal efficiency. S. bradyi showed a large range of variations in trypsin activity and low amylase activity, resulting in a low average A:T ratio, but its O:N ratio was intermediate. Variations in O:N and O:P ratios were related to differences in the nutritional strategy of the different species, based on literature data concerning the anatomy of their mouth parts and their selectivity for animal material in mixed-food experiments.     

Testing a post-copulatory pre-zygotic reproductive barrier in a passerine species pair

Sexual selection may drive speciation, but most research focuses on pre-copulatory sexual selection, overlooking post-copulatory processes. Post-copulatory sexual selection in allopatric populations could drive divergence in post-copulatory pre-zygotic (PCPZ) phenotypes, limiting gene flow upon secondary contact. Here, we performed in vitro experiments examining one potential PCPZ barrier between two closely related passerine species, house sparrows (Passer domesticus) and Spanish sparrows (Passer hispaniolensis). In birds, crossing in the vagina may be particularly challenging for sperm, so we tested the effect of female reproductive tract fluids on sperm swimming speed and motility. If a PCPZ barrier exists at this stage of the fertilization process, heterospecific female fluids are predicted to reduce sperm swimming speed or motility relative to conspecific female fluid. We found that house sparrow female fluids affected the two species’ sperm asymmetrically, depending on the control sperm velocity and male species. Overall, however, sperm performed equally in conspecific and heterospecific female fluids, and the species had similar sperm morphology and sperm swimming performance. Low divergence in PCPZ phenotypes between species, perhaps because post-copulatory sexual selection is stabilizing or only moderately strong in these taxa, may be insufficient to cause an overall PCPZ barrier. Reinforcement may be unlikely to drive PCPZ barriers for this species pair, because relatively effective pre-copulatory barriers exist between the species, and because hybrids can be quite successful. Testing the role of PCPZ barriers in birds with more divergent PCPZ phenotypes will improve our understanding of speciation in passerines.     

Repeated geographic divergence in behavior: a case study employing phenotypic trajectory analyses

Environmental effects on behavior have long been a focus of behavioral ecologists. Among the important drivers of behavior is predation environment, which can include the presence/absence of predators, differences in resource availability, and variation in individual density. Environments with predators are often more ecologically complex and “risky” than those without predators. Populations from these environments are sometimes more active and explorative than populations from low-risk, less complex environments. To date, most comparative studies of behavior are limited to within-species comparisons of populations from divergent environments, but neglect comparisons between species following speciation, thus limiting our understanding of post-speciation behavioral evolution. Brachyrhaphis fishes provide an ideal system for studying correlations between divergent environments and behavior within and between species. Here, we test for differences in two behavioral traits—activity and exploration —between sister species Brachyrhaphis roseni and Brachyrhaphis terrabensis that occur in divergent predation environments. Species differed in activity and exploration, with higher activity and exploration levels in populations that co-occur with predators. Furthermore, we found drainage-by-species interactions, indicating that the nature of divergence varied geographically. Using the recently developed phenotypic trajectory analysis (PTA), we quantified this difference and found that, while the geographically isolated populations of sister species tended to evolve in parallel, the magnitude of divergence between species differed between drainages. Our results highlight the utility of PTA for multivariate behavioral data and corroborate past predictions that complex and risky environments are correlated with increased activity and exploration levels and that divergence continues post-speciation.     

Geographic variation in allele frequency of the gamete recognition protein M7 lysin throughout a mosaic blue mussel hybrid zone

Divergence at gamete recognition loci is hypothesized to result in speciation in broadcast spawning invertebrates. Many gamete recognition loci evolve quickly and show patterns of positive selection, yet the advantage of divergence is rarely known. M7 lysin is a sperm protein in the Mytilus edulis species complex that shows evidence of adaptive evolution. This locus is polymorphic with two distinct clades within Mytilus galloprovincialis, one of which, the D clade, shows the strongest signal of positive selection. We tested whether the geographic patterns in allele frequency were consistent with the hypothesis that positive selection on D clade alleles (G _D ) was due to reinforcement. Populations of M. edulis showed little evidence of introgression of G _D alleles, but there was no consistent evidence of reproductive character displacement. The lack of consistent patterns expected of reinforcement suggests that another mechanism is likely responsible for the adaptive divergence of M7 lysin.     

Eurythoe complanata (Polychaeta: Amphinomidae), the ‘cosmopolitan’ fireworm, consists of at least three cryptic species

Eurythoe complanata (Pallas 1766) has been considered a cosmopolitan species with a great morphological similarity across its geographic range. To elucidate whether E. complanata is actually a single species, genetic (cytochrome oxidase subunit I and allozymes) and morphological differences were compared among specimens from the Pacific, Caribbean, and South Atlantic Oceans. Large levels of COI divergence (10–22%) and diagnostic allozyme loci identified three cryptic species: one in the eastern Pacific and two in the Atlantic, with one being morphologically differentiated and found only in islands. COI sequences between Pacific and Atlantic lineages were much more divergent than those of other transisthmian invertebrates, indicating their split before the Panama Isthmus closure or a faster evolutionary rate of COI for this species. The existence of two Atlantic species may be a consequence of parapatric speciation followed by a secondary invasion or even a sympatric speciation in the Atlantic oceanic islands.     

Biofluorescence as a survey tool for cryptic marine species

As ecosystems come under increasing anthropogenic pressure, rare species face the highest risk of extinction. Paradoxically, data necessary to evaluate the conservation status of rare species are often lacking because of the challenges of detecting species with low abundance. One group of fishes subject to this undersampling bias are those with cryptic body patterns. Twenty‐one percent of cryptic fish species assessed for their extinction risk (International Union for Conservation of Nature [IUCN]) are data deficient. We developed a nondestructive method for surveying cryptically patterned marine fishes based on the presence of biofluorescence (underwater biofluorescence census, UBC). Blue LED torches were used to investigate how widespread biofluorescence was in cryptic reef fishes in the Coral Triangle region. The effectiveness of UBC to generate abundance data was tested on a data‐deficient pygmy seahorse species (Hippocampus bargibanti) and compared with data obtained from standard underwater visual census (UVC) surveys. We recorded 95 reef fish species displaying biofluorescence, 73 of which had not been previously described as biofluorescent. Of those fish with cryptic patterns, 87% were biofluorescent compared with 9% for noncryptic fishes. The probability of species displaying biofluorescence was 70.9 times greater for cryptic species than for noncryptic species. Almost twice the number of H. bargibanti was counted using the UBC compared with UVC. For 2 triplefin species (Ucla xenogrammus, Enneapterygius tutuilae), the abundance detected with UBC was triple that detected with UVC. The UBC method was effective at finding cryptic species that would otherwise be difficult to detect and thus will reduce interobserver variability inherent to UVC surveys. Biofluorescence is ubiquitous in cryptic fishes, making this method applicable across a wide range of species. Data collected using UBC could be used with multiple IUCN criteria to assess the extinction risk of cryptic species. Adopting this technique will enhance researchers’ ability to survey cryptic species and facilitate management and conservation of cryptic marine species.     

Krill evolution and the Antarctic ocean currents: evidence of vicariant speciation as inferred by molecular data

The phylogenetic relationships of the Antarctic krill Euphausia superba, the key species in the Antarctic food web, and other Antarctic and sub-Antarctic cuphausiids have been investigated using the 16S ribosomal mitochondrial gene. The phylogenetic reconstructions indicated that the Antarctic species form a monophyletic clade separated by the non-Antarctic species. The results revealed a large genetic divergence between the Antarctic (E. superba and E. crystallorophias) and sub-Antarctic species (E. vallentini). The time of separation between these species, estimated from the molecular data, is around 20 million years ago, which is comparable with the geological time of the formation of a circum-Antarctic water circulation and the Antarctic Polar Frontal Zone. The euphausiid molecular phylogeny therefore represents evidence for vicariant speciation.     

Reproductive isolation among morphotypes of the Atlantic seastar species<Emphasis Type="Italic"> Zoroaster fulgens</Emphasis> (Asteroidea: Echinodermata)

Zoroaster fulgens is a slope-dwelling seastar species that is distributed throughout the Atlantic Ocean. Studies into the population structure and systematics of marine animals have increasingly found that species with a reported cosmopolitan distribution are, in fact, collections of closely related cryptic or sibling species. In the Porcupine Seabight, three morphotypes of Z. fulgens can be found that have a distribution that is stratified by depth. This study investigates the genetic divergence between these morphotypes using sections of the cytochrome oxidase 1 (COI) and 16S regions of the mitochondrial genome. Bathymetrically separated morphotypes of Z. fulgens are reproductively isolated over distances of approximately 1 km while gene flow occurs among morphotypes, along isobaths, over distances of approximately 900 km. Reproductive isolation on the continental slope may have occurred as a result of selection exerted by gradients of depth-correlated physical factors, such as pressure and temperature. However, allopatric speciation with subsequent range expansion may also explain the observed patterns of genetic divergence. Further investigation of radiation within this group may provide important information on the evolution of slope species. Taxonomic revision of the genus is required.Communicated by J.P. Thorpe, Port Erin     

Cryptic species within the commercially most important lobster in the tropical Atlantic, the spiny lobster Panulirus argus

Panulirus argus (Latreille in Ann Mus Hist Nat Paris 3:388–395, 1804) is the lobster of greatest economic importance throughout its distribution. In this study, mitochondrial (Cytochrome Oxidase I and 16S ribosomal genes) and nuclear (Adenine Nucleotide Transporter gene) sequences were used to evaluate the taxonomic status of P. argus sampled from five sites in the Caribbean Sea and nine sites in the Southwest Atlantic. Phylogenetic analyses indicate that lobsters from the two regions form two monophyletic groups with a molecular divergence similar to that observed between distinct congeneric lobster species and much larger than that found between conspecific lobster populations. Therefore, the Caribbean and the Southwest Atlantic lobster populations originally attributed to P. argus belong to different species, with an estimated time of isolation of around 16 Million years. An important consequence of these findings is that the fisheries of spiny lobsters from the Caribbean and the Southwest Atlantic species must be managed separately.     

Phylogeography of the European stalked barnacle (Pollicipes pollicipes): identification of glacial refugia

The alternation of glacial and interglacial events during the Pleistocene has produced changes in species distribution ranges leading to bottlenecks and alterations of patterns of gene flow. The European stalked barnacle, Pollicipes pollicipes, is a sessile pedunculate cirripede that inhabits the rocky intertidal frame, from Senegal to the northwestern coast of France. In this work, we have analyzed a fragment of the mitochondrial gene cytochrome c oxidase subunit I for 569 individuals of P. pollicipes in order to investigate whether the shifts in climatic conditions that occurred during the Pleistocene influenced the current pattern of distribution of genetic variation of P. pollicipes. A pre-last glacial maximum pattern of demographic expansion was found, in concordance with many other North Atlantic marine species. On the other hand, three potential glacial refugia were identified: North African coasts, northwestern Iberian Peninsula and English Channel/Brittany.     

Population genetic structure of the neon damselfish (<Emphasis Type="Italic">Pomacentrus coelestis</Emphasis>) in the northwestern Pacific Ocean

The population genetic structure of the neon damselfish (Pomacentrus coelestis) in the northwestern Pacific Ocean was revealed by the hypervariable control region of the mitochondrial gene (343 bp). In total, 170 individuals were sampled from 8 localities distributed between Taiwan and Japan, and 71 haplotypes were obtained through sequence alignment. High haplotype diversity (h = 0.956 ± 0.008) with low nucleotide diversity (π = 0.010 ± 0.006) was observed, and the results of the mismatch distribution test suggested that a historical population expansion after a period of population bottleneck might have occurred among P. coelestis populations. Based on the results of the UPGMA tree and AMOVA (Φct = 0.193, P < 0.05) analyses, fish populations from eight localities could be divided into two groups: one includes populations from localities around mainland Japan, and the other includes those from Okinawa and southern Taiwan. A genetic break was found between populations from mainland Japan and Okinawa, and this break was congruent with the pattern of phenotypic variations documented in previous studies. This evidence supports the latitudinal variation of reproductive traits among P. coelestis populations likely being genetically based. It is suggested that the changes in sea level and sea surface temperatures during past glaciations might have resulted in population bottlenecks in P. coelestis and the modern populations in the northern West Pacific are likely the results of recolonization after such events. The Kuroshio Current acts not only as a vehicle for larval transport along its pathway (between populations in southern Taiwan and Okinawa) but also as a barrier for larval dispersal across the Kuroshio axis (between populations in mainland Japan and Okinawa). Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Seasonal variations in leaf capturing of particulate matter, surface wettability and micromorphology in urban tree species

The seasonal changes in leaf particulate matter (PM) accumulation, surface wettability and micromorphology in urban tree species, including Sophora japonica (S. japonica), Platanus acerifolia (P. acerifolia) and Cedrus deodara (C. deodara), were studied during a single growing season. The three species showed distinct seasonal trends in PM accumulation, increasing from spring to autumn (or winter) even during the rainy season, but at different rates. During the study, the leaf PM retention amount of P. acerifolia, a species with ridged leaf surfaces, was significantly higher than that of S. japonica and C. deodara, species with waxy leaf surfaces. The contact angles of water droplets on leaves decreased with leaf age except on the abaxial surface of S. japonica, which remained non-wettable or highly non-wettable throughout the growing season; the decrease in the contact angle on adaxial surface of S. japonica was greater when compared with P. acerifolia and C. deodara. A significant and negative relationship existed between leaf PM retention amounts and contact angles on adaxial surface of leaves of all three species. The increase in wettability, probably caused when epicuticular wax was destroyed by mechanical and chemical abrasion, seemed to be the main factor leading to seasonal variations in leaf PM accumulation.