Hirokazu Toju

High-coverage ITS primers for the DNA-based identification of ascomycetes and basidiomycetes in environmental samples.

The kingdom Fungi is estimated to include 1.5 million or more species, playing key roles as decomposers, mutualists, and parasites in every biome on the earth. To comprehensively understand the diversity and ecology of this huge kingdom, DNA barcoding targeting the internal transcribed spacer (ITS) region of the nuclear ribosomal repeat has been regarded as a prerequisite procedure. By extensively surveying ITS sequences in public databases, we designed new ITS primers with improved coverage across diverse taxonomic groups of fungi compared to existing primers. An in silico analysis based on public sequence databases indicated that the newly designed primers matched 99% of ascomycete and basidiomycete ITS taxa (species, subspecies or varieties), causing little taxonomic bias toward either fungal group. Two of the newly designed primers could inhibit the amplification of plant sequences and would enable the selective investigation of fungal communities in mycorrhizal associations, soil, and other types of environmental samples. Optimal PCR conditions for the primers were explored in an in vitro investigation. The new primers developed in this study will provide a basis for ecological studies on the diversity and community structures of fungi in the era of massive DNA sequencing.

Imbalance of Predator and Prey Armament: Geographic Clines in Phenotypic Interface and Natural Selection

The escalation of defensive/offensive arms is ubiquitous in prey‐predator evolutionary interactions. However, there may be a geographically varying imbalance in the armaments of participating species that affects the outcome of local interactions. In a system involving the Japanese camellia (Camellia japonica) and its obligate seed predator, the camellia weevil (Curculio camelliae), we investigated the geographic variation in physical defensive/offensive traits and that in natural selection on the plant’s defense among 17 populations over a 700‐km‐wide area in Japan. The sizes of the plant defensive apparatus (pericarp thickness) and the weevil offensive apparatus (rostrum length) clearly correlated with each other across populations. Nevertheless, the balance in armaments between the two species was geographically structured. In the populations for which the balance was relatively advantageous for the plant’s defense, natural selection on the trait was stronger because in the other populations, most plant individuals were too vulnerable to resist the attacks of the weevil, and their seeds were infested independent of pericarp thickness. We also found that the imbalance between the defensive/offensive armaments and the intensity of natural selection showed clear latitudinal clines. Overall, our results suggest that the imbalance of armament between sympatric prey and predator could determine the strength of local selection and that climatic conditions could affect the local and overall trajectory of coevolutionary arms races.

Diversity and infection prevalence of endosymbionts in natural populations of the chestnut weevil: relevance of local climate and host plants.

Many insects are ubiquitously associated with multiple endosymbionts, whose infection patterns often exhibit spatial and temporal variations. How such endosymbiont variations are relevant to local adaptation of the host organisms is of ecological interest. Here, we report a comprehensive survey of endosymbionts in natural populations of the chestnut weevil Curculio sikkimensis, whose larvae are notorious pests of cultivated chestnuts and also infest acorns of various wild oaks. From 968 insects representing 55 localities across the Japanese Archipelago and originating from 10 host plant species, we identified six distinct endosymbiont lineages, namely Curculioniphilus, Sodalis, Serratia, Wolbachia, Rickettsia and Spiroplasma, at different infection frequencies (96.7%, 12.8%, 82.3%, 82.5%, 28.2% and 6.8%, respectively) and with different geographical distribution patterns. Multiple endosymbiont infections were very common; 3.18 ± 0.61 (ranging from 1.74 to 5.50) endosymbionts per insect on average in each of the local populations. Five pairs of endosymbionts (Curculioniphilus–Serratia, Curculioniphilus–Wolbachia, Sodalis–Rickettsia, Wolbachia–Rickettsia and Rickettsia–Spiroplasma) co‐infected the same host individuals more frequently than expected, while infections with Serratia and Wolbachia were negatively correlated to each other. Infection frequencies of the endosymbionts were significantly correlated with climatic and ecological factors: for example, higher Sodalis, Wolbachia and Rickettsia infections at localities of higher temperature; lower Wolbachia and Rickettsia infections at localities of greater snowfall; and higher Curculioniphilus, Sodalis, Serratia, Wolbachia and Rickettsia infections on acorns than on chestnuts. These patterns are discussed in relation to potential host–endosymbiont co‐evolution via local adaptation across geographical populations.

Assembly of complex plant–fungus networks

Species in ecological communities build complex webs of interaction. Although revealing the architecture of these networks is fundamental to understanding ecological and evolutionary dynamics in nature, it has been difficult to characterize the structure of most species-rich ecological systems. By overcoming this limitation through next-generation sequencing technology, we herein uncover the network architecture of below-ground plant–fungus symbioses, which are ubiquitous to terrestrial ecosystems. The examined symbiotic network of a temperate forest in Japan includes 33 plant species and 387 functionally and phylogenetically diverse fungal taxa, and the overall network architecture differs fundamentally from that of other ecological networks. In contrast to results for other ecological networks and theoretical predictions for symbiotic networks, the plant–fungus network shows moderate or relatively low levels of interaction specialization and modularity and an unusual pattern of ‘nested’ network architecture. These results suggest that species-rich ecological networks are more architecturally diverse than previously recognized.

Two New Computational Methods for Universal DNA Barcoding: A Benchmark Using Barcode Sequences of Bacteria, Archaea, Animals, Fungi, and Land Plants

Taxonomic identification of biological specimens based on DNA sequence information (a.k.a. DNA barcoding) is becoming increasingly common in biodiversity science. Although several methods have been proposed, many of them are not universally applicable due to the need for prerequisite phylogenetic/machine-learning analyses, the need for huge computational resources, or the lack of a firm theoretical background. Here, we propose two new computational methods of DNA barcoding and show a benchmark for bacterial/archeal 16S, animal COX1, fungal internal transcribed spacer, and three plant chloroplast (rbcL, matK, and trnH-psbA) barcode loci that can be used to compare the performance of existing and new methods. The benchmark was performed under two alternative situations: query sequences were available in the corresponding reference sequence databases in one, but were not available in the other. In the former situation, the commonly used “1-nearest-neighbor” (1-NN) method, which assigns the taxonomic information of the most similar sequences in a reference database (i.e., BLAST-top-hit reference sequence) to a query, displays the highest rate and highest precision of successful taxonomic identification. However, in the latter situation, the 1-NN method produced extremely high rates of misidentification for all the barcode loci examined. In contrast, one of our new methods, the query-centric auto-k-nearest-neighbor (QCauto) method, consistently produced low rates of misidentification for all the loci examined in both situations. These results indicate that the 1-NN method is most suitable if the reference sequences of all potentially observable species are available in databases; otherwise, the QCauto method returns the most reliable identification results. The benchmark results also indicated that the taxon coverage of reference sequences is far from complete for genus or species level identification in all the barcode loci examined. Therefore, we need to accelerate the registration of reference barcode sequences to apply high-throughput DNA barcoding to genus or species level identification in biodiversity research.

Community composition of root-associated fungi in a Quercus-dominated temperate forest: “codominance” of mycorrhizal and root-endophytic fungi

In terrestrial ecosystems, plant roots are colonized by various clades of mycorrhizal and endophytic fungi. Focused on the root systems of an oak-dominated temperate forest in Japan, we used 454 pyrosequencing to explore how phylogenetically diverse fungi constitute an ecological community of multiple ecotypes. In total, 345 operational taxonomic units (OTUs) of fungi were found from 159 terminal-root samples from 12 plant species occurring in the forest. Due to the dominance of an oak species (Quercus serrata), diverse ectomycorrhizal clades such as Russula, Lactarius, Cortinarius, Tomentella, Amanita, Boletus, and Cenococcum were observed. Unexpectedly, the root-associated fungal community was dominated by root-endophytic ascomycetes in Helotiales, Chaetothyriales, and Rhytismatales. Overall, 55.3% of root samples were colonized by both the commonly observed ascomycetes and ectomycorrhizal fungi; 75.0% of the root samples of the dominant Q. serrata were so cocolonized. Overall, this study revealed that root-associated fungal communities of oak-dominated temperate forests were dominated not only by ectomycorrhizal fungi but also by diverse root endophytes and that potential ecological interactions between the two ecotypes may be important to understand the complex assembly processes of belowground fungal communities.

Diversification of endosymbiosis: replacements, co-speciation and promiscuity of bacteriocyte symbionts in weevils

The processes and mechanisms underlying the diversification of host–microbe endosymbiotic associations are of evolutionary interest. Here we investigated the bacteriocyte-associated primary symbionts of weevils wherein the ancient symbiont Nardonella has experienced two independent replacement events: once by Curculioniphilus symbiont in the lineage of Curculio and allied weevils of the tribe Curculionini, and once by Sodalis-allied symbiont in the lineage of grain weevils of the genus Sitophilus. The Curculioniphilus symbiont was detected from 27 of 36 Curculionini species examined, the symbiont phylogeny was congruent with the host weevil phylogeny, and the symbiont gene sequences exhibited AT-biased nucleotide compositions and accelerated molecular evolution. These results suggest that the Curculioniphilus symbiont was acquired by an ancestor of the tribe Curculionini, replaced the original symbiont Nardonella, and has co-speciated with the host weevils over evolutionary time, but has been occasionally lost in several host lineages. By contrast, the Sodalis-allied symbiont of Sitophilus weevils exhibited no host–symbiont co-speciation, no AT-biased nucleotide compositions and only moderately accelerated molecular evolution. These results suggest that the Sodalis-allied symbiont was certainly acquired by an ancestor of the Sitophilus weevils and replaced the original Nardonella symbiont, but the symbiotic association must have experienced occasional re-associations such as new acquisitions, horizontal transfers, replacements and/or losses. We detected Sodalis-allied facultative symbionts in populations of the Curculionini weevils, which might represent potential evolutionary sources of the Sodalis-allied primary symbionts. Comparison of these newcomer bacteriocyte-associated symbiont lineages highlights potential evolutionary trajectories and consequences of novel symbionts after independent replacements of the same ancient symbiont.

How are plant and fungal communities linked to each other in belowground ecosystems? A massively parallel pyrosequencing analysis of the association specificity of root-associated fungi and their host plants.

In natural forests, hundreds of fungal species colonize plant roots. The preference or specificity for partners in these symbiotic relationships is a key to understanding how the community structures of root-associated fungi and their host plants influence each other. In an oak-dominated forest in Japan, we investigated the root-associated fungal community based on a pyrosequencing analysis of the roots of 33 plant species. Of the 387 fungal taxa observed, 153 (39.5%) were identified on at least two plant species. Although many mycorrhizal and root-endophytic fungi are shared between the plant species, the five most common plant species in the community had specificity in their association with fungal taxa. Likewise, fungi displayed remarkable variation in their association specificity for plants even within the same phylogenetic or ecological groups. For example, some fungi in the ectomycorrhizal family Russulaceae were detected almost exclusively on specific oak (Quercus) species, whereas other Russulaceae fungi were found even on “non-ectomycorrhizal” plants (e.g., Lyonia and Ilex). Putatively endophytic ascomycetes in the orders Helotiales and Chaetothyriales also displayed variation in their association specificity and many of them were shared among plant species as major symbionts. These results suggest that the entire structure of belowground plant–fungal associations is described neither by the random sharing of hosts/symbionts nor by complete compartmentalization by mycorrhizal type. Rather, the colonization of multiple types of mycorrhizal fungi on the same plant species and the prevalence of diverse root-endophytic fungi may be important features of belowground linkage between plant and fungal communities.

Phylogeography and the geographic cline in the armament of a seed-predatory weevil: effects of historical events vs. natural selection from the host plant

Japanese camellia (Camellia japonica) and its seed predator, the camellia weevil (Curculio camelliae), provide a notable example of a geographic mosaic of coevolution. In the species interaction, the offensive trait of the weevil (rostrum length) and the defensive trait of the plant (pericarp thickness) are involved in a geographically‐structured arms race, and these traits and selective pressures acting on the plant defence vary greatly across a geographical landscape. To further explore the geographical structure of this interspecific interaction, we tested whether the geographical variation in the weevil rostrum over an 800‐km range along latitude is attributed to local natural selection or constrained by historical (phylogeographical) events of local populations. Phylogeographical analyses of the mitochondrial DNA sequences of the camellia weevil revealed that this species has experienced differentiation into two regions, with a population bottleneck and subsequent range and/or population expansion within each region. Although these phylogeographical factors have affected the variation in rostrum length, analyses of competing factors for the geographical variation revealed that this pattern is primarily determined by the defensive trait of the host plant rather than by the effects of historical events of populations and a climatic factor (annual mean temperature). Thus, our study suggests the overwhelming strength of coevolutionary selection against the effect of historical events, which may have limited local adaptation.

Adaptive divergence of scaling relationships mediates the arms race between a weevil and its host plant

Coevolution of exaggerated morphologies between insects and plants is a well-known but poorly understood phenomenon in evolutionary biology. In the antagonistic interaction between a seed-predatory insect, the camellia weevil (Curculio camelliae), and its host plant, Japanese camellia (Camellia japonica), we examined the evolutionary trajectory of an exaggerated offensive trait of the weevil (rostrum length) in terms of scaling relationship. Sampling throughout Japan revealed that the ratio of the rostrum length to overall body size was correlated with the ratio of the pericarp thickness to overall fruit size across the localities. We found a geographical interpopulation divergence in a parameter pertaining to the allometric equation of rostrum length (the coefficient a in y=axb, where y and x denote rostrum and body lengths, respectively), and the pattern of geographical differentiation in the allometric coefficient was closely correlated with the variation in the pericarp thickness of Japanese camellia. Our results provide a novel example of a geographically diverged scaling relationship in an insect morphology resulting from a coevolutionary arms race with its host plant.