Phylogeography of Phytophagous Weevils and Plant Species in Broadleaved Evergreen Forests: A Congruent Genetic Gap between Western and Eastern Parts of Japan

Evolution of secondary metabolites, morphological structures and associated gene expression patterns in galls induced by four closely related aphid species on a host plant species

Gall-forming insects induce various types of galls on their host plants by altering gene expression in host plant organs, and recent studies have been conducted for gene expression in galls. However, the evolutionary trajectories of gene expression patterns and the resulting phenotypes have not yet been studied using multiple related species. We investigated the speciation and the diversification process of galls induced by four closely related aphid species (Hormaphidini) on a host plant species (Hamamelis japonica) by examining the phylogenetic congruence between the geographical divergences of aphids and the host plant, and by comparing their gene expression patterns and resulting phenotypes. Phylogenetic analysis of aphids and the host plant showed that geographical isolation among host plant populations has interrupted gene flow in aphids and accelerated the speciation process. The concentration of phenolics and the complexity of the internal structure of galls were correlated with the expression levels of genes for the biosynthesis of phenolics and morphogenesis respectively. These results suggest that the expression levels of genes for the biosynthesis of phenolics and morphogenesis have evolutionarily increased in galls accelerated by the speciation process of aphids due to the distribution change of the host plant, leading to the related phenotypic evolution. Our study showed the evolutionary process of phenotypic traits in galls in the wild from both gene expression and actual phenotype levels.

Phylogeographical evidence for historical long-distance dispersal in the flightless stick insect Ramulus mikado

Exploring how organisms overcome geographical barriers to dispersal is a fundamental question in biology. Passive long-distance dispersal events, although infrequent and unpredictable, have a considerable impact on species range expansions. Despite limited active dispersal capabilities, many stick insect species have vast geographical ranges, indicating that passive long-distance dispersal is vital for their distribution. A potential mode of passive dispersal in stick insects is via the egg stage within avian digestive tracts, as suggested by experimental evidence. However, detecting such events under natural conditions is challenging due to their rarity. Therefore, to indirectly assess the potential of historical avian-mediated dispersal, we examined the population genetic structure of the flightless stick insect Ramulus mikado across Japan, based on a multifaceted molecular approach [cytochrome oxidase subunit I (COI) haplotypes, nuclear simple sequence repeat markers and genome-wide single nucleotide polymorphisms]. Subsequently, we identified unique phylogeographic patterns, including the discovery of identical COI genotypes spanning considerable distances, which substantiates the notion of passive long-distance genotypic dispersal. Overall, all the molecular data revealed the low and mostly non-significant genetic differentiation among populations, with identical or very similar genotypes across distant populations. We propose that long-distance dispersal facilitated by birds is the plausible explanation for the unique phylogeographic pattern observed in this flightless stick insect.

Genetic diversity and demography of Bufo japonicus and B. torrenticola (Amphibia: Anura: Bufonidae) influenced by the Quaternary climate

The Quaternary climate affected the present species richness and geographic distribution patterns of amphibians by limiting their activities during the glacial period. The present study examined the phylogenetic relationships of Japanese toads (Bufo japonicus and B. torrenticola) and the demography of each lineage from the past to the present based on mitochondrial sequences and ecological niche models. Japanese toads are a monophyletic group with two main clades (clades A and B). Clade A represents B. j. formosus, including three clades (clades A1, A2, and A3). Clade B contains three clades, two of which corresponded to B. j. japonicus (clades B1 and B2) and the other to B. torrenticola. Clade B2 and B. torrenticola made a sister group, and, thus, B. j. japonicus is paraphyletic. Clades A and B diverged in the late Miocene 5.7 million years ago (Mya) during the period when the Japanese archipelago was constructed. The earliest divergence between the three clades of clade A was estimated at 1.8 Mya. Clades A1 and A2 may have diverged at 0.8 Mya, resulting from the isolation in the multiple different refugia; however, the effects of the glacial climate on the divergence events of clade A3 are unclear. Divergences within clade B occurred from the late Pliocene to the early Pleistocene (3.2-2.2 Mya). Niche similarity between the parapatric clade in clade B (clades B1 and B2) indicated their allopatric divergence. It was suggested that niche segregation between B. japonicus and B. torrenticola contributed to a rapid adaptation of B. torrenticola for lotic breeding. All clade of Japanese toads retreated to each refugium at a low elevation in the glacial period, and effective population sizes increased to construct the current populations after the Last Glacial Maximum. Furthermore, we highlight the areas of climate stability from the last glacial maximum to the present that have served as the refugia of Japanese toads and, thus, affected their present distribution patterns.

Timing Is Everything. Temporal and Spatial Niche Segregation in Curculio spp. (Coleoptera: Curculionidae) Associated with Oak Trees

Oak seed predatory weevils occurring in Poland are prone to increased interspecific competition due to the limited number of Quercus species, compared to southern Europe, in which they can develop. Therefore, analyses on the preferences of three weevil species for acorn sizes chosen for reproduction, as well as on reproductive period duration, were performed. Cafeteria-type experiments were set for females of three species associated with one oak species. Females were allowed to choose and oviposit in acorns of different sizes and growth stages. Research revealed statistically significant differences between the masses of acorns chosen for oviposition by females of Curculio glandium (the biggest), C. pellitus (medium), and C. venosus (the smallest). Studied weevils also differed in terms of the beginning of the reproductive period, which corresponded with the increasing mass of growing acorns. Moreover, C. glandium was observed to be the only species to perform radial egg galleries and lay a considerably higher and varied number of eggs. The results support the hypothesis of a strategy aimed at reducing interspecific competition between Curculio spp. in terms of limited host plant species number.

Independent allopatric polyploidizations shaped the geographical structure and initial stage of reproductive isolation in an allotetraploid fern, Lepisorus nigripes (Polypodiaceae)

Although polyploidy is pervasive and its evolutionary significance has been recognized, it remains unclear how newly formed polyploid species become established. In particular, the impact of multiple origins on genetic differentiation among populations of a polyploid species and whether lineages of independent origins have different evolutionary potentials remain open questions. We used population genetic and phylogenetic approaches to identify genetic differentiation between lineages with independent origins within an allotetraploid fern, Lepisorus nigripes. A total of 352 individuals from 51 populations were collected throughout the distribution range. To examine the genetic structure, multilocus genotyping, Bayesian population structure analysis, and neighbor-net analysis were carried out using single-copy nuclear genes. Phylogenetic trees were constructed to detect recurrent polyploid origins. Proportions of abortive spores were analysed as the measure of postzygotic reproductive isolation. Two genetically distinct lineages, the East-type and the West-type, were distributed mainly in the eastern and western parts, respectively, of the Japanese archipelago. Phylogenetic analyses indicated independent origins of these types and detected additional independent origins within each type. We also revealed limited genetic recombination between both types, even in their sympatric regions. F₁ hybrids between the East- and West-types showed a reduction in fertility. It is likely that the East- and West-types formed independently in the eastern and western parts of Japan, respectively. The limited genetic recombination and reduced fertility of hybrids suggest that the two types are at an incipient stage of speciation. Two polyploid lineages with independent geographic origins could develop reproductive isolation barrier(s).

Phylogenetic analyses reveal that simian foamy virus isolated from Japanese Yakushima macaques (Macaca fuscata yakui) is distinct from most of Japanese Hondo macaques (Macaca fuscata fuscata)

Japanese macaque (Macaca fuscata) is an indigenous Old World monkey (OWM) species that inhabits the Japanese archipelago. There are two subspecies of Japanese macaque: Yakushima macaque (M. f. yakui) which inhabits Yakushima Island exclusively, and Hondo macaque (M. f. fuscata) which inhabits the mainland of Japan. Yakushima macaque is considered to be branched off from a certain parental macaque group that had inhabited the mainland of Japan. However, the process of sub-speciation of the Yakushima macaque is still unclear at present. In this study, to gain new insight into the process of sub-speciation of Japanese macaque, we utilized the simian foamy virus (SFV) as a marker. SFVs are found in virtually all primates except humans and undergo species-specific cospeciation with the hosts. The phylogenetic analysis of conserved regions of the env gene in SFVs remarkably resembled that of the OWMs with high statistical confidence. The phylogenetic analyses also indicated that there are four (1-4) genotypes among Asian OWMs investigated. SFVs derived from Asian OWMs except Yakushima macaque were classified as genotypes 1-3, whereas SFVs isolated from all Yakushima macaques and one Hondo macaque were classified as genotype 4. Interestingly, genotype 4 was firstly branched off from the rest of the genotypes, which might indicate that the macaques infected with genotype 4 SFV were derived from the "older" population of Japanese macaques. The high prevalence of genotype 4 SFVs among Yakushima macaque might reflect the possibility that they are a descendant of the population settled earlier, which has been geographically isolated in Yakushima Island.

Genetic data reveals a complex history of multiple admixture events in presently allopatric wild gingers (Asarum spp.) showing intertaxonomic clinal variation in calyx lobe length

Clinal variation is a major pattern of observed phenotypic diversity and identifying underlying demographic processes is a necessary step to understand the establishment of clinal variation. The wild ginger series Sakawanum (genus Asarum) comprises four taxa, which exhibit intertaxonomic clinal variation in calyx lobe length across two continental islands isolated by a sea strait. To test alternative hypotheses of the evolutionary history and to determine the implications for the formation of clinal variation, we conducted approximate Bayesian computation (ABC) analysis and ecological niche modeling (ENM). ABC analysis indicated that the scenario assuming multiple admixture events was strongly supported. This scenario assumed two admixture events occurred between morphologically distinct taxa, likely leading to the generation of intermediate taxa. One of the admixture events was estimated to have occurred during the last glacial maximum (LGM), during which the taxa were estimated to have formed a common refugia in southern areas by ENM analysis. Although four taxa are currently distributed allopatrically on different islands and trans-oceanic dispersal appears unlikely, the formation of a land bridge and the geographic range shift to refugia would have allowed secondary contact between previously isolated taxa. This study suggests that clinal variation can be shaped by demographic history including multiple admixtures due to climatic oscillations.

Approximate Bayesian computation analysis of EST-associated microsatellites indicates that the broadleaved evergreen tree Castanopsis sieboldii survived the Last Glacial Maximum in multiple refugia in Japan

Climatic changes have played major roles in plants' evolutionary history. Glacial oscillations have been particularly important, but some of their effects on plants' populations are poorly understood, including the numbers and locations of refugia in Asian warm temperate zones. In the present study, we investigated the demographic history of the broadleaved evergreen tree species Castanopsis sieboldii (Fagaceae) during the last glacial period in Japan. We used approximate Bayesian computation (ABC) for model comparison and parameter estimation for the demographic modeling using 27 EST-associated microsatellites. We also performed the species distribution modeling (SDM). The results strongly support a demographic scenario that the Ryukyu Islands and the western parts in the main islands (Kyushu and western Shikoku) were derived from separate refugia and the eastern parts in the main islands and the Japan Sea groups were diverged from the western parts prior to the coldest stage of the Last Glacial Maximum (LGM). Our data indicate that multiple refugia survived at least one in the Ryukyu Islands, and the other three regions of the western and eastern parts and around the Japan Sea of the main islands of Japan during the LGM. The SDM analysis also suggests the potential habitats under LGM climate conditions were mainly located along the Pacific Ocean side of the coastal region. Our ABC-based study helps efforts resolve the demographic history of a dominant species in warm temperate broadleaved forests during and after the last glacial period, which provides a basic model for future phylogeographical studies using this approach.

Evolution of Trichobaris (Curculionidae) in relation to host plants: Geometric morphometrics, phylogeny and phylogeography

The family Curculionidae (Coleoptera), the "true" weevils, have diversified tightly linked to the evolution of flowering plants. Here, we aim to assess diversification at a lower taxonomic level. We analyze the evolution of the genus Trichobaris in association with their host plants. Trichobaris comprises eight to thirteen species; their larvae feed inside the fruits of Datura spp. or inside the stem of wild and cultivated species of Solanaceae, such as potato, tobacco and tomato. We ask the following questions: (1) does the rostrum of Trichobaris species evolve according to the plant tissue used to oviposit, i.e., shorter rostrum to dig in stems and longer to dig in fruits? and (2) does Trichobaris diversify mainly in relation to the use of Datura species? For the first question, we estimated the phylogeny of Trichobaris based on four gene sequences (nuclear 18S and 28S rRNA genes and mitochondrial 16S rRNA and COI genes). Then, we carried out morphogeometric analyses of the Trichobaris species using 75 landmarks. For the second question, we calibrated a COI haplotype phylogeny using a constant rate of divergence to infer the diversification time of Trichobaris species, and we traced the host plant species on the haplotype network. We performed an ancestral state reconstruction analysis to infer recent colonization events and conserved associations with host plant species. We found that ancestral species in the Trichobaris phylogeny use the stem of Solanum plants for oviposition and display weak sexual dimorphism of rostrum size, whereas other, more recent species of Trichobaris display sexual dimorphism in rostrum size and use the fruits of Datura species, and a possible reversion to use the stem of Solanaceae was detected in one Trichobaris species. The use of Datura species by Trichobaris species is widely distributed on haplotype networks and restricted to Trichobaris species that originated ca. 5 ± 1.5 Ma. Given that the origin of Trichobaris is estimated to be ca. 6 ± 1.5 Ma, it is likely that Datura has played a role in its diversification.

Environmental factors contribute to the formation and maintenance of the contact zone observed in deciduous broad-leaved tree species in Japan

Contact zones are defined as areas where populations from different refugia meet during a postglacial expansion and distinct DNA lineages are mixedly distributed. In Japan, contact zones of various plants and animals were reported from the Kinki-Chugoku region. These contact zones appear to be maintained without any drastic topographic barriers such as those observed in the Alps and Pyrenees Mountains. In this study, the mechanisms underlying the formation and/or maintenance of these contact zones were investigated using six deciduous broad-leaved tree species (Carpinus laxiflora, C. tschonoskii, C. japonica, Magnolia obovata, Padus grayana, and Euonymus oxyphyllus). First, the precise location of the contact zones was examined by intensive genetic analysis of the six species. Second, the relationships between the geographic location of the contact zone and various environmental factors, including climate and topography, were investigated by generalized additive models to reveal the mechanisms of the formation and maintenance of the contact zones. As a result, four of the six examined plant species clearly showed a geographically common contact zone in Hyogo Prefecture and its adjacent areas. The results of the generalized additive models indicate that the pattern of low habitat suitability estimated by ecological niche modeling was the most important factor for determining the location of the common contact zone. These results suggest that areas with low habitat suitability in Hyogo Prefecture restrict the migration and gene flow of the four species in this region, and thus, they maintain the pattern of the contact zones. This study suggests that there are major effects of habitat suitability on the formation and maintenance of the contact zones.

A comparative analysis of genetic differentiation across six shared willow host species in leaf- and bud-galling sawflies

Genetic divergence and speciation in plant-feeding insects could be driven by contrasting selection pressures imposed by different plant species and taxa. While numerous examples of host-associated differentiation (HAD) have been found, the overall importance of HAD in insect diversification remains unclear, as few studies have investigated its frequency in relation to all speciation events. One promising way to infer the prevalence and repeatability of HAD is to estimate genetic differentiation in multiple insect taxa that use the same set of hosts. To this end, we measured and compared variation in mitochondrial COI and nuclear ITS2 sequences in population samples of leaf-galling Pontania and bud-galling Euura sawflies (Hymenoptera: Tenthredinidae) collected from six Salix species in two replicate locations in northern Fennoscandia. We found evidence of frequent HAD in both species complexes, as individuals from the same willow species tended to cluster together on both mitochondrial and nuclear phylogenetic trees. Although few fixed differences among the putative species were found, hierarchical AMOVAs showed that most of the genetic variation in the samples was explained by host species rather than by sampling location. Nevertheless, the levels of HAD measured across specific pairs of host species were not correlated in the two focal galler groups. Hence, our results support the hypothesis of HAD as a central force in herbivore speciation, but also indicate that evolutionary trajectories are only weakly repeatable even in temporally overlapping radiations of related insect taxa.

Range expansion and lineage admixture of the Japanese evergreen tree Machilus thunbergii in central Japan

We investigated the range expansion histories of Machilus thunbergii populations in the Kinki region of central Japan on the basis of nuclear microsatellite data. In the Kinki region, M. thunbergii is typically found in the coastal area, with some fragmented populations inland, around Lake Biwa. Phylogenetic and Bayesian clustering analysis (STRUCTURE analysis) revealed that the inland populations have different genetic components between the west and east sides of Lake Biwa. The population located on the north side of the lake has an admixture of the two genetically differentiated lineages, contributing to an increase in the genetic diversity of the population. Populations around Lake Biwa had lost rare alleles and the F value obtained from STRUCTURE analysis was lower in the coastal populations than in the lake populations. These results suggest that populations around Lake Biwa experienced a bottleneck due to a founder effect during the initial migration to the lake and that glacial refugia of M. thunbergii in the Kinki region existed along the coast.

Genetic differentiation and genetic diversity of Castanopsis (Fagaceae), the dominant tree species in Japanese broadleaved evergreen forests, revealed by analysis of EST-associated microsatellites

The broadleaved evergreen forests of the East Asian warm temperate zone are characterised by their high biodiversity and endemism, and there is therefore a need to extend our understanding of its genetic diversity and phylogeographic patterns. Castanopsis (Fagaceae) is one of the dominant tree species in the broadleaved evergreen forests of Japan. In this study we investigate the genetic diversity, genetic structure and leaf epidermal morphology of 63 natural populations of C. sieboldii and C. cuspidata, using 32 Expressed Sequence Tag associated microsatellites. The overall genetic differentiation between populations was low (G_ST = 0.069 in C. sieboldii and G_ST = 0.057 in C. cuspidata). Neighbor-joining tree and Bayesian clustering analyses revealed that the populations of C. sieboldii and C. cuspidata were genetically clearly differentiated, a result which is consistent with the morphology of their epidermal cell layers. This suggests that C. sieboldii and C. cuspidata should be treated as independent species, although intermediate morphologies are often observed, especially at sites where the two species coexist. The higher level of genetic diversity observed in the Kyushu region (for both species) and the Ryukyu Islands (for C. sieboldii) is consistent with the available fossil pollen data for Castanopsis-type broadleaved evergreen trees during the Last Glacial Maximum and suggests the existence of refugia for Castanopsis forests in southern Japan. Within the C. sieboldii populations, Bayesian clustering analyses detected three clusters, in the western and eastern parts of the main islands and in the Ryukyu Islands. The west-east genetic differentiation observed for this species in the main islands, a pattern which is also found in several plant and animal species inhabiting Castanopsis forests in Japan, suggests that they have been isolated from each other in the western and eastern populations for an extended period of time, and may imply the existence of eastern refugia.

Geographic distributions of sexual and parthenogenetic forms of the weevil Scepticus insularis (Coleoptera: Curculionidae) in Japan

Geographic ranges of parthenogenetic animals often differ from those of their sexual conspecifics, providing valuable information about evolutionary histories. We studied the geographic distributions of sexual and parthenogenetic forms of the weevil Scepticus insularis (Roelofs) (Coleoptera, Curculionidae) on Honshu, Shikoku, and Kyushu Islands, Japan. We collected Scepticus insularis at 18 sites in the northern half of Honshu; sequenced 633 bp of the mitochondrial ND2 (NADH dehydrogenase subunit 2) gene for 102 Honshu beetles; and constructed a ML gene tree based on these sequences and sequences previously reported for 104 conspecific beetles from Hokkaido of known reproductive mode. The ND2 haplotypes fell into two distinct clades (I and II) separated by 9.5% nucleotide sequence divergence. The haplotypes from all parthenogenetic females on Hokkaido and those from all-female populations on Honshu invariably belonged to clade I, whereas sexual females on Hokkaido and some of the females in male/female populations on Honshu belonged to clade II. Haplotypes of males were in clade II, except for two Hokkaido males having a clade-I haplotype. Excluding these two enigmatic males, we treated individuals with clade-I haplotypes as parthenogenetic, and those with clade-II haplotypes as sexual. Parthenogenetic beetles were widespread throughout Hokkaido and central and northern Honshu, whereas sexuals were found only in three separate areas on Hokkaido and one on Honshu, often sympatrically with parthenogens. The distribution of S. insularis is congruent with the pattern termed geographic parthenogenesis.

Phylogeography of Saproxylic and Forest Floor Invertebrates from Tallaganda, South-eastern Australia

The interaction between physiogeographic landscape context and certain life history characteristics, particularly dispersal ability, can generate predictable outcomes for how species responded to Pleistocene (and earlier) climatic changes. Furthermore, the extent to which impacts of past landscape-level changes ‘scale-up’ to whole communities has begun to be addressed via comparative phylogeographic analyses of co-distributed species. Here we present an overview of a body of research on flightless low-mobility forest invertebrates, focusing on two springtails and two terrestrial flatworms, from Tallaganda on the Great Dividing Range of south-eastern Australia. These species are distantly-related, and represent contrasting trophic levels (i.e., slime-mold-grazers vs. higher-level predators). However, they share an association with the dead wood (saproxylic) habitat. Spatial patterns of intraspecific genetic diversity partly conform to topography-based divisions that circumscribe five ‘microgeographic regions’ at Tallaganda. In synthesizing population processes and past events that generated contemporary spatial patterns of genetic diversity in these forest floor invertebrates, we highlight cases of phylogeographic congruence, pseudo-congruence, and incongruence. Finally, we propose conservation-oriented recommendations for the prioritisation of areas for protection.