Macroecology of Aquatic Insects: A Quantitative Analysis of Taxonomic Richness and Composition in the Andes Mountains of Northern Ecuador

To determine the spatial dynamics of Neotropical lotic insect species, specimens were collected from 41 streams on the eastern and western flanks of the Andes Mountains in Ecuador. We examined the manner in which taxonomic richness and composition differed with elevation, latitude, and versant. Statistical analyses were limited to 5 families (Ephemeroptera: Baetidae, Leptohyphidae, Oligoneuriidae; Heteroptera: Naucoridae; Coleoptera: Elmidae), comprising 32 genera and 85 species, for which identifications or morphospecies assignments were reliable. Assessment of taxonomic diversity was based on the richness of genera and species at each site. In addition, each site was characterized by species richness within each of 4 families with more than 10 species. The effects of versant and transect on composition and structure were family‐specific. Mean site differences between versants in elmid richness at generic and specific levels depended on transect. Only baetid richness was affected by versant and transect in a consistent manner. Variation among sites in composition based on all genera and species was captured using multidimensional scaling (MDS). Subsequent MANOVAs based on scores from MDS axes revealed that differences between versants were much stronger in the southern transect, although transects differed from each other at specific and generic levels. A Jaccard's similarity matrix was computed for each family to reflect the spatial organization of taxonomic composition. Mesogeographic patterns of species composition for each of the four families were correlated (Mantel analysis) at both the regional level and at the level of the entire study area. At the regional level, the only pair of families to exhibit correlated patterns of species composition was elmids and naucorids in the southwestern region. The pattern of species composition for each family was correlated with the patterns for one or more other families at the level of the entire study area. Thus, spatial dynamics of species composition was similar for the families examined, suggesting that the Andes exert a consistent influence on species distributions within families, regardless of ordinal affiliation. At a local scale, however, the way in which taxonomic composition changed with latitude and versant was family‐specific. Mayflies, the most vagile of the taxa studied, had the highest percentage of species overlap between versants. Of three genera of Naucoridae collected, species of Ambrysus, of probable Mexican origin, were found only on the eastern versant, corroborating other evidence that the genus is recent in South America. Moreover, dispersion by Ambrysus across the Andes Mountains may not have occurred, as it has for Cryphocricos and Limnocoris, which are of probable South American origin.     

Taxonomic and structural similarities in soil oribatid communities

We examined the taxonomic and structural concordance between oribatid mite communities at different sites from around the world using meta-analysis. The study explored similarities in the taxonomie composition and species abundance distributions among oribatid mite communities using 25 published species lists from Africa, Australia, Europe. Japan. Malaysia, New Zealand and the United States, The final data set contained 253 genera from 97 families in 43 superfamilies, A second data set, which included only studies where trees were present at the sites, contained 234 genera, 93 families and 41 superfamilies. Composition was analysed at taxonomic levels higher than species. Only two families and four superfamilies were present in > 90% of all studies. Two genera, six families and seven superfamilies were found in all sites with trees. Common families were different among habitats suggesting that habitat preferences may be expressed at the family level. There were very strong relationships between oribatid mite species richness and richness at higher taxonomic levels with r² values for regression relationships > 0.8. Few species at a site were from the same genus, family or superfamily. Species abundance distributions were remarkably consistent between studies with the three most dominant species each constituting > 10% of the oribatid individuals. Overall, similar factors may determine the relative proportions of soil-dwelling oribatid species at a site.     

Latitudinal variation in the shape of the species body size distribution: an analysis using freshwater fishes

Many taxonomic and ecological assemblages of species exhibit a right-skewed body size-frequency distribution when characterized at a regional scale. Although this distribution has been frequently described, factors influencing geographic variation in the distribution are not well understood, nor are mechanisms responsible for distribution shape. In this study, variation in the species body size-frequency distributions of 344 regional communities of North American freshwater fishes is examined in relation to latitude, species richness, and taxonomic composition. Although the distribution of all species of North American fishes is right-skewed, a negative correlation exists between latitude and regional community size distribution skewness, with size distributions becoming left-skewed at high latitudes. This relationship is not an artifact of the confounding relationship between latitude and species richness in North American fishes. The negative correlation between latitude and regional community size distribution skewness is partially due to the geographic distribution of families of fishes and apparently enhanced by a nonrandom geographic distribution of species within families. These results are discussed in the context of previous explanations of factors responsible for the generation of species size-frequency distributions related to the fractal nature of the environment, energetics, and evolutionary patterns of body size in North American fishes.     

Taxonomic Composition of Hymenopteran Parasitoid Assemblages from Agromyzid Leaf-miners Sampled in Central Argentina

This paper deals with the taxonomic composition of parasitoid assemblages reared from agromyzid leaf-miners of 28 species collected in Cordoba, Central Argentina. The relative representation of each parasitic family was estimated, taking into consideration their species richness and abundance. Sample size affected the total species richness in each complex, and also that of Eulophidae and Pteromalidae, indicating that not only size but also taxonomic structure of parasitic assemblages are influenced by sampling effort. 3 superfamilies and 5 families of Hymenoptera Parasitica were represented, with Eulophidae (Chalcidoidea) specifically and numerically dominant both in the total community of parasitoids and in the parasitic assemblages of individual hosts. Chrysonotomyia, Chrysocharis (Eulophidae), Halticoptera (Pteromalidae) and Phaedrotoma (Braconidae) were the most relevant genera in terms of number of species and number of specimens reared. The parasitoids obtained from samples taken in agroecosystems differed from the global spectrum of taxa mainly by showing a numerical dominance of Braconidae, despite a noticeable specific impoverishment of this family.     

Traits explain community disassembly and trophic contraction following experimental environmental change

Many ecosystems are currently undergoing dramatic changes in biodiversity due to habitat loss and climate change. Responses to global change at the community level are poorly understood, as are the impacts of community disassembly on ecosystem‐level processes. Uncertainties remain regarding the patterns of extirpation and persistence under single vs. multiple forms of environmental change. Here, we use a trait‐based and food web approach to examine the effects of experimentally changing moisture, temperature and habitat ‘openness’ on a functionally important group of microarthropods associated with a boreal forest floor bryosphere (detrital moss) system. Overall, the outcome of community disassembly was mediated by the correlation between our environmental factors and species traits, particularly body size. Minor increases in summer temperatures maintained greater species richness, whereas drought stress had a significant negative effect on community‐level abundance and richness. These effects were reflected in modifications to the community‐wide body‐size spectra. Habitat openness alleviated biodiversity loss in the larger‐bodied species of the most abundant taxonomic group, but did not fully mitigate the effects of drought. The most striking result of this experiment was an overall contraction of the food web among persistent species under drought stress (i.e. those not extirpated by environmental change). These results suggest that major changes in boreal microarthropod community structure are likely to occur in response to common forms of global change. Moreover, the contraction in trophic structure even amongst tolerant species suggests that ecosystem function within the bryosphere can be altered by environmental change.     

A unified model explains commonness and rarity on coral reefs

Abundance patterns in ecological communities have important implications for biodiversity maintenance and ecosystem functioning. However, ecological theory has been largely unsuccessful at capturing multiple macroecological abundance patterns simultaneously. Here, we propose a parsimonious model that unifies widespread ecological relationships involving local aggregation, species‐abundance distributions, and species associations, and we test this model against the metacommunity structure of reef‐building corals and coral reef fishes across the western and central Pacific. For both corals and fishes, the unified model simultaneously captures extremely well local species‐abundance distributions, interspecific variation in the strength of spatial aggregation, patterns of community similarity, species accumulation, and regional species richness, performing far better than alternative models also examined here and in previous work on coral reefs. Our approach contributes to the development of synthetic theory for large‐scale patterns of community structure in nature, and to addressing ongoing challenges in biodiversity conservation at macroecological scales.     

A metagenetic approach for revealing community structure of marine planktonic copepods

Marine planktonic copepods are an ecologically important group with high species richness and abundance. Here, we propose a new metagenetic approach for revealing the community structure of marine planktonic copepods using 454 pyrosequencing of nuclear large subunit ribosomal DNA. We determined an appropriate similarity threshold for clustering pyrosequencing data into molecular operational taxonomic units (MOTUs) using an artificial community containing 33 morphologically identified species. The 99% similarity threshold had high species‐level resolution for MOTU clustering but overestimated species richness. The artificial community was appropriately clustered into MOTUs at 97% similarity, with little inflation in MOTU numbers and with relatively high species‐level resolution. The number of sequence reads of each MOTU was correlated with dry weight of that taxon, suggesting that sequence reads could be used as a proxy for biomass. Next, we applied the method to field‐collected samples, and the results corresponded reasonably well with morphological analysis of these communities. Numbers of MOTUs were well correlated with species richness at 97% similarity, and large numbers of sequence reads were generally observed in MOTUs derived from species with large biomass. Further, MOTUs were successfully classified into taxonomic groups at the family level at 97% similarity; similar patterns of species richness and biomass were revealed within families with metagenetic and morphological analyses. At the 99% similarity threshold, MOTUs with high proportions of sequence reads were identified as biomass‐dominant species in each field‐collected sample. The metagenetic approach reported here can be an effective tool for rapid and comprehensive assessment of copepod community structure.     

Non‐neutrality in forest communities: evolutionary and ecological determinants of tree species abundance distributions

The unified neutral theory of biodiversity and biogeography provides a promising framework that can be used to integrate stochastic and ecological processes operating in ecological communities. Based on a mechanistic non‐neutral model that incorporates density‐dependent mortality, we evaluated the deviation from a neutral pattern in tree species abundance distributions and explored the signatures of historical and ecological processes that have shaped forest biomes. We compiled a dataset documenting species abundance distributions in 1168 plots encompassing 16 973 tree species across tropical, temperate, and boreal forests. We tested whether deviations from neutrality of species abundance distributions vary with climatic and historical conditions, and whether these patterns differ among regions. Non‐neutrality in species abundance distributions was ubiquitous in tropical, temperate, and boreal forests, and regional differences in patterns of non‐neutrality were significant between biomes. Species abundance evenness/unevenness caused by negative density‐dependent or abiotic filtering effects had no clear macro‐scale climatic drivers, although temperature was non‐linearly correlated with species abundance unevenness on a global scale. These findings were not significantly biased by heterogeneity of plot data (the differences of plot area, measurement size, species richness, and the number of individuals sampled). Therefore, our results suggest that environmental filtering is not universally increasing from warm tropical to cold boreal forests, but might affect differently tree species assembly between and within biomes. Ecological processes generating particularly dominant species in local communities might be idiosyncratic or region‐specific and may be associated with geography and climate. Our study illustrates that stochastic dynamical models enable the analysis of the interplay of historical and ecological processes that influence community assemblies and the dynamics of biodiversity.     

Metacommunity‐scale biodiversity regulation and the self‐organised emergence of macroecological patterns

There exist a number of key macroecological patterns whose ubiquity suggests that the spatio‐temporal structure of ecological communities is governed by some universal mechanisms. The nature of these mechanisms, however, remains poorly understood. Here, we probe spatio‐temporal patterns in species richness and community composition using a simple metacommunity assembly model. Despite making no a priori assumptions regarding biotic spatial structure or the distribution of biomass across species, model metacommunities self‐organise to reproduce well‐documented patterns including characteristic species abundance distributions, range size distributions and species area relations. Also in agreement with observations, species richness in our model attains an equilibrium despite continuous species turnover. Crucially, it is in the neighbourhood of the equilibrium that we observe the emergence of these key macroecological patterns. Biodiversity equilibria in models occur due to the onset of ecological structural instability, a population‐dynamical mechanism. This strongly suggests a causal link between local community processes and macroecological phenomena.     

Dynamic relationships between body size,species richness,abundance, and energy use in a shallow marine epibenthic faunal community

We study the temporal variation in the empirical relationships among body size (S), species richness (R), and abundance (A) in a shallow marine epibenthic faunal community in Coliumo Bay, Chile. We also extend previous analyses by calculating individual energy use (E) and test whether its bivariate and trivariate relationships with S and R are in agreement with expectations derived from the energetic equivalence rule. Carnivorous and scavenger species representing over 95% of sample abundance and biomass were studied. For each individual, body size (g) was measured and E was estimated following published allometric relationships. Data for each sample were tabulated into exponential body size bins, comparing species‐averaged values with individual‐based estimates which allow species to potentially occupy multiple size classes. For individual‐based data, both the number of individuals and species across body size classes are fit by a Weibull function rather than by a power law scaling. Species richness is also a power law of the number of individuals. Energy use shows a piecewise scaling relationship with body size, with energetic equivalence holding true only for size classes above the modal abundance class. Species‐based data showed either weak linear or no significant patterns, likely due to the decrease in the number of data points across body size classes. Hence, for individual‐based size spectra, the SRA relationship seems to be general despite seasonal forcing and strong disturbances in Coliumo Bay. The unimodal abundance distribution results in a piecewise energy scaling relationship, with small individuals showing a positive scaling and large individuals showing energetic equivalence. Hence, strict energetic equivalence should not be expected for unimodal abundance distributions. On the other hand, while species‐based data do not show unimodal SRA relationships, energy use across body size classes did not show significant trends, supporting energetic equivalence.     

Size-dependent species richness: trends within plant communities and across latitude

We examine how species richness and species‐specific plant density (number of species and number of individuals per species, respectively) vary within community size frequency distributions and across latitude. Communities from Asia, Africa, Europe, and North, Central and South America were studied (60°4′N–41°4′S latitude) using the Gentry data base. Log–log linear stem size (diameter) frequency distributions were constructed for each community and the species richness and species‐specific plant density within each size class were determined for each frequency distribution. Species richness in the smallest stem size class correlated with the Y‐intercepts (β‐values) of the regression curves describing each log–log linear size distributions. Two extreme community types were identified (designated as type A and type B). Type A communities had steep size distributions (i.e. large β‐values), log–log linear species‐richness size distributions, low species‐specific plant density distributions, and a small size class (2–4 cm) containing the majority of all species but rarely conspecifics of the dominant tree species. Type B communities had shallow size distributions (i.e. small β‐values), more or less uniform (and low) size class species‐ richness and species‐specific density distributions and size‐dominant species resident in the smallest size class. Type A communities were absent in the higher latitudes but increased in number towards the equator, i.e. in the smallest size class, species richness increased (and species‐specific density decreased) towards the tropics. Based on our survey of type A and type B communities (and their intermediates), species richness evinces size‐dependent and latitudinal trends, i.e. species richness increased with decreasing body size and most species increasingly reside in the smallest plant size class towards the tropics. Across all latitudes, a trade‐off exists between the number of species and the number of individuals per species residing in the smaller size classes.     

Abundance, diversity and body size: patterns from a grassland arthropod community

1. The empirical relationships among body size, species richness and number of individuals may give insight into the factors controlling species diversity and the relative abundances of species. To determine these relationships, we sampled the arthropods of grasslands and savannahs at Cedar Creek, MN using sweep nets (90 525 individuals of 1225 species) and pitfall traps (12 721 individuals of 92 species). Specimens were identified, enumerated and measured to determine body size.
2. Both overall and within abundant taxonomic orders, species richness and numbers of individuals peaked at body sizes intermediate for each group. Evolution could create unimodal diversity patterns by random diversification around an ancestral body size or from size-dependent fitness differences. Local processes such as competition or predation could also create unimodal diversity distributions.
3. The average body size of a species depended significantly on its taxonomic order, but on contemporary trophic role only within the context of taxonomic order.
4. Species richness ( S _i) within size classes was related to the number of individuals ( I _i) as S _i =  I i^0·5. This relationship held across a 100 000-fold range of body sizes. Within size classes, abundance distributions of size classes were all similar power functions. A general rule of resource division, together with similar minimum population sizes, is sufficient to generate the relationship between species richness and number of individuals.
5. Smaller bodied species had slightly shallower abundance distributions and may, in general, persist at lower densities than larger species.
6. Our results suggest there may be fewer undescribed small arthropod species than previously thought and that most undescribed species will be smaller than arthropods.     

Multitaxonomic diversity patterns along a desert riparian-upland gradient

Riparian areas are noted for their high biodiversity, but this has rarely been tested across a wide range of taxonomic groups. We set out to describe species richness, species abundance, and community similarity patterns for 11 taxonomic groups (forbs & grasses, shrubs, trees, solpugids, spiders, scarab beetles, butterflies, lizards, birds, rodents, and mammalian carnivores) individually and for all groups combined along a riparian-upland gradient in semiarid southeastern Arizona, USA. Additionally, we assessed whether biological characteristics could explain variation in diversity along the gradient using five traits (trophic level, body size, life span, thermoregulatory mechanism, and taxonomic affiliation). At the level of individual groups diversity patterns varied along the gradient, with some having greater richness and/or abundance in riparian zones whereas others were more diverse and/or abundant in upland zones. Across all taxa combined, riparian zones contained significantly more species than the uplands. Community similarity between riparian and upland zones was low, and beta diversity was significantly greater than expected for most taxonomic groups, though biological traits explained little variance in diversity along the gradient. These results indicate heterogeneity amongst taxa in how they respond to the factors that structure ecological communities in riparian landscapes. Nevertheless, across taxonomic groups the overall pattern is one of greater species richness and abundance in riparian zones, coupled with a distinct suite of species.     

The saproxylic beetle assemblage associated with different host trees in Southwest China

Dead wood is a habitat for many insects and other small animals, some of which may be rare or endangered and in need of effective protection. In this paper, saproxylic beetle assemblages associated with different host trees in the subtropical forests in southwestern China were investigated. A total of 277 species (1 439 specimens) in 36 beetle families were collected from 117 dead wood samples, of which 101 samples were identified and respectively belonged to 12 tree genera. The number of saproxylic beetle species varied greatly among logs of different tree genera, with the highest diversity on logs of Juglans. Generally, broad‐leaved trees had a higher richness and abundance of saproxylic species than coniferous trees. Cluster analysis revealed that assemblages from broad‐leaved tree genera were generally similar (except for Betula) and assemblages from coniferous trees formed another distinct cluster. The subsequent indicator analysis proposed that there are different characteristic species for different cluster groups of host tree genera. In our study, log diameter has no positive influence on beetle species density. Conversely, comparisons of individual‐based rarefaction curves suggested that beetle species richness was highest in the small diameter class both in coniferous and broad‐leaved tree genera. With increased wood decay, proportion of habitat specialists (saproxylic beetles living on one tree genus) decreased, whereas proportion of habitat generalists (living on more than three tree genera) increased. The beetle species density was found to be higher in early stages, and decreased in later stages as well. A negative influence of altitude on saproxylic beetle species richness and abundance was detected. It was indicated that different tree genera and altitudes possibly display cross effects in modulating the altitudinal distribution and host preference of the beetles.     

陕南秦巴山区药用植物群落物种多样性研究

采用样方法对陕南秦巴山区药用植物群落物种多样性进行了研究。结果显示,该地区药用植物共计105科268属361种,其中药用裸子植物5科6属7种,药用被子植物92科254属344种,药用蕨类植物7科8属10种;含1属的科、含2~5属的科和含1种的属的药用植物数量较多,所占比例较高;该地区各植被类型的物种多样性变化各具特点,其中,物种丰富度指数(S)、Shannon-Wiener多样性指数(H)及Simpson多样性指数(D)表现为:乔木层最小、灌木层较大、草本层最大,但Pielou物种均匀度指数(J)的变化较为复杂;药用植物群落总的物种丰富度、乔木层物种丰富度和草本层物种丰富度均随海拔的升高表现出先升高后下降的趋势,且物种丰富度最大值出现在海拔1000~1200 m的山地常绿落叶阔叶林和温性针阔混交林的过渡区域,但灌木层物种丰富度与海拔之间无相关关系。本研究结果为阐明山地物种多样性的分布格局提供了基础资料,也为该区域药用植物的保护和管理提供参考。     

Seasonal drought limits tree species across the Neotropics

Within the tropics, the species richness of tree communities is strongly and positively associated with precipitation. Previous research has suggested that this macroecological pattern is driven by the negative effect of water‐stress on the physiological processes of most tree species. This implies that the range limits of taxa are defined by their ability to occur under dry conditions, and thus in terms of species distributions predicts a nested pattern of taxa distribution from wet to dry areas. However, this ‘dry‐tolerance’ hypothesis has yet to be adequately tested at large spatial and taxonomic scales. Here, using a dataset of 531 inventory plots of closed canopy forest distributed across the western Neotropics we investigated how precipitation, evaluated both as mean annual precipitation and as the maximum climatological water deficit, influences the distribution of tropical tree species, genera and families. We find that the distributions of tree taxa are indeed nested along precipitation gradients in the western Neotropics. Taxa tolerant to seasonal drought are disproportionally widespread across the precipitation gradient, with most reaching even the wettest climates sampled; however, most taxa analysed are restricted to wet areas. Our results suggest that the ‘dry tolerance' hypothesis has broad applicability in the world's most species‐rich forests. In addition, the large number of species restricted to wetter conditions strongly indicates that an increased frequency of drought could severely threaten biodiversity in this region. Overall, this study establishes a baseline for exploring how tropical forest tree composition may change in response to current and future environmental changes in this region.     

Macroecology of local insect communities

The inter-relationships between animal body weight, range size, species richness and abundance are currently the basis of macroecology. Using 41 099 insects sampled from 31 Acacia tree canopies in north-east Tanzania, we first documented the basic macroecological patterns. The relationship between body weight and both species richness and abundance was polygonal with the highest insect species richness and abundance occurring at intermediate body weights. Across individual tree communities, the most statistically robust relationships were found between range size, abundance and species richness and they were all linear. In a second part, we focused on the positive abundance-range size relationship and we could test predictions of six of the eight proposed hypotheses to explain this widely documented pattern of community structure. The relationship is most likely explained by the metapopulation hypothesis stating that with more patches being occupied, local abundance in a given patch increases due to a higher rate of immigration from nearby patches. In addition, we found high slopes for the species-area relationship, typical of island systems and thus it seems reasonable to characterise Acacia trees in the savannah as habitat islands for insects.     

Seasonal Incidence and Diversity Pattern of Avian Communities in the Bangalore University Campus,India

The present study deals with the species abundance, diversity and species richness of avian communities in the Bangalore University Campus (BUC), Bengaluru, India. One hundred and six species of birds belonging to 42 families under 68 genera were recorded. Shannon–Wiener’s and Fisher’s alpha diversities, species evenness, species richness of bird communities, number of bird species and percentage of population density of birds between various seasons in the BUC differed significantly between the study years. Of these bird species, the relative abundance (6.96 %) and species distribution ratio (0.070) of Psittacula krameri were highest, whereas relative abundance (0.04 %) and species distribution ratio (0.002) of Coracias benghalensis were lowest. The existing 32 species of flowering plants/trees belonging to 29 genera under 14 families in the campus are used for perching by birds. Moreover 29 species of flowering plants/trees belonging to 24 genera under 16 families depend on birds for pollination and/or seed dispersal. Occurrence of greater bird diversity and abundance of avian communities were recorded highest in the winter season in the BUC premises. In the different seasons, the BUC had varying community structure of birds between the study years. BUC suffers from numerous threats namely grass cutting, fire and grazing of domestic animals. Conservation methods needed for habitat management are restoration of vegetation and wetlands, and increase plant and tree diversity to protect the ecosystem of BUC habitat and to preserve its diversity of avifauna.     

Effects of Tree Genotypic Diversity and Species Diversity on the Arthropod Community Associated with Big‐leaf Mahogany

Despite potential interactive effects of plant species and genotypic diversity (SD and GD, respectively) on consumers, studies have usually examined these effects separately. We evaluated the individual and combined effects of tree SD and mahogany (Swietenia macrophylla) GD on the arthropod community associated with mahogany. We conducted this study within the context of a tree diversity experiment consisting of 74 plots with 64 saplings/plot. We sampled 24 of these plots, classified as monocultures of mahogany or polycultures of four species (including mahogany). Within each plot type, mahogany was represented by either one or four maternal families. We surveyed arthropods on mahogany and estimated total arthropod abundance and species richness, as well as abundance and richness separately for herbivorous and predatory arthropods. Overall tree SD and mahogany GD had positive effects on total arthropod species richness and abundance on mahogany, and also exerted interactive effects on total species richness (but not abundance). Analyses conducted by trophic level group showed contrasting patterns; SD positively influenced herbivore species richness but not abundance, and did not affect either predator richness or abundance. GD influenced predator species richness but not abundance, and did not influence herbivore abundance or richness. There were interactive effects of GD and SD only for predator species richness. These results provide evidence that intra‐ and inter‐specific plant diversity exert interactive controls on associated consumer communities, and that the relative importance of SD and GD may vary among higher trophic levels, presumably due to differences in the underlying mechanisms or consumer traits.     

Host effects on ectomycorrhizal fungal communities: insight from eight host species in mixed conifer-broadleaf forests

To advance our understanding of host effects on the community structure of ectomycorrhizal fungi (EMF), EMF communities were compared among different host species, genera and families in two mixed conifer-broadleaf forests in Japan. Using molecular identification methods we examined EMF root tips of eight coexisting species belonging to six genera (three families): Abies and Tsuga (Pinaceae), Betula and Carpinus (Betulaceae) and Fagus and Quercus (Fagaceae). In total, 205 EMF species were detected, and the total richness was estimated to exceed 300 species using major estimators. Of the 55 EMF species occurring three or more times, eight showed significantly biased host preference. A Mantel test showed a significantly negative correlation between EMF community similarity and host taxonomic distance. Detrended correspondence analysis separated EMF communities mainly by host taxonomic and successional status. Thus, EMF communities are similar on hosts with similar taxonomic and successional status. A significant proportion of EMF exhibited host specificity, which may contribute to the extremely diverse EMF community in conifer-broadleaf forests.