MIG-seq: an effective PCR-based method for genome-wide single-nucleotide polymorphism genotyping using the next-generation sequencing platform

Detection of reproductive interference between closely related Salvia species with small-scale separated distributions by multifaceted pollination and molecular analyses

Reproductive interference, an interspecific interaction in reproductive process that exerts an adverse effect, has gained attention as a contributing factor in promoting exclusive distributions between closely related species. However, detailed studies on the possibility of reproductive interference between native plants are still lacking, presumably because strong reproductive interference can rapidly realize exclusive distributions, leaving the two species apparently independent. Salvia japonica and S. lutescens are found in separate localities at a small scale, although their distributions overlap at a large scale. We investigated the possibility of reproductive interference between them through field surveys, hand-pollination experiments, evaluation of hybrid fertility, cpDNA and nrDNA genotyping, and genome-wide DNA analysis. The field survey results did not reveal apparent negative interaction in competition for pollinator services. Mixed pollination with conspecific pollen and counterpart pollen reduced seed set in S. japonica, and hybrid progeny produced by mixed pollination were less than 20% as fertile compared to the pure species. The DNA genotyping results suggested the possibility of hybridization where their distributions overlap, and the genome-wide DNA analysis results showed clear genetic differentiation between the two species as well as the existence of hybrids. These results suggest that bi-directional reproductive interference between S. japonica and S. lutescens may have led to their present separated distributions at a small scale.

Genetic diversification of allohexaploid Brassica hybrids (AABBCC) using a fertile octoploid with excessive C genome set (AABBCCCC)

MAIN CONCLUSION

Using octoploid somatic hybrids with excessive C genome sets, AABBCCCC, a diverse allohexaploid, AABBCC, was produced by C genome reduction through subsequent crossing with various AABB cultivars. Even when somatic hybrids are produced, the plants that are produced are rarely in themselves an innovative crop. In this study, we used somatic hybrids of Brassica juncea (AABB) and B. oleracea (CC) as model cases for the genetic diversification of the somatic hybrids. One cell of 'Akaoba Takana' (B. juncea) and two cells of 'Snow Crown' (B. oleracea) were fused to create several somatic hybrids with excessive C genomes, AABBCCCC. Using AABBCCCC somatic hybrids as mother plants and crossing with 'Akaoba Takana', the AABBCC progenies were generated. When these AABBCC plants were self-fertilized, and flow cytometric (FCM) analysis was performed on the next generations, differences in the relative amount of genome size variation were observed, depending on the different AABBCCCC parents used for AABBCC creation. Further self-progeny was obtained for AABBCC plants with a theoretical allohexaploid DNA index by FCM. However, as the DNA indices of the progeny populations varied between plants used and aneuploid individuals still occurred in the progeny populations, it was difficult to say that the allohexaploid genome was fully stabilized. Next, to obtain genetic diversification of the allohexaploid, different cultivars of B. juncea were crossed with AABBCCCC, resulting in diverse AABBCC plants. Genetic diversity can be further expanded by crossbreeding plants with different AABBCC genome sets. Although genetic stability is necessary to ensure in the later generations, the results obtained in this study show that the use of somatic hybrids with excess genomes is an effective strategy for creating innovative crops.

Secondary contact zone and genetic introgression in closely related haplodiploid social spider mites

How frequently hybridisation and gene flow occur in the contact zones of diverging taxa is important for understanding the speciation process. Stigmaeopsis sabelisi and Stigmaeopsis miscanthi high-aggression form (hereafter, S. miscanthi HG) are haplodiploid, social spider mites that infest the Chinese silver grass, Miscanthus sinensis. These two species are closely related and parapatrically distributed in Japan. In mountainous areas, S. sabelisi and S. miscanthi HG are often found in the highlands and lowlands, respectively, suggesting that they are in contact at intermediate altitudes. It is estimated that they diverged from their common ancestors distributed in subtropical regions (south of Japan) during the last glacial period, expanded their distribution into the Japanese Archipelago, and came to have such a parapatric distribution (secondary contact). As their reproductive isolation is strong but incomplete, hybridisation and genetic introgression are expected at their distributional boundaries. In this study, we investigated their spatial distribution patterns along the elevation on Mt. Amagi using male morphological differences, and investigated their hybridisation status using single-nucleotide polymorphisms by MIG-seq. We found their contact zone at altitudes of 150-430 m, suggesting that their contact zone is prevalent in the parapatric area, which is in line with a previous study. Interspecific mating was predicted based on the sex ratio in the contact zone. No obvious hybrids were found, but genetic introgression was detected although it was extremely low.

Testing the applicability of environmental DNA metabarcoding to landscape genetics

Landscape genetics is a field dealing with local genetic differences and contributes to strategic conservation planning. Recently, environmental DNA (eDNA) metabarcoding has proven useful not only for detecting species but also for assessing genetic diversity and genetic structure on a large scale such as in phylogeography. However, it remains unclear whether eDNA analysis also has sufficient power to perform the landscape genetics, which focuses on a local scale. To reveal the applicability of eDNA to landscape genetics, we conducted an eDNA metabarcoding analysis of the mitochondrial DNA D-loop region of the fluvial sculpin Cottus nozawae in the upper Sorachi River in Japan and compared the results with inferences based on traditional tissue-based approaches by the same D-loop region and genome-wide SNP data. As a result, the spatial distribution of haplotypes was generally consistent between the eDNA- and tissue-based approaches. In addition, the genetic differentiation statistics calculated using eDNA and tissue samples were highly correlated when comparing both in the D-loop region. The removal of low-frequency reads or the conversion to semi-quantitative rankings of eDNA data did not alter the correlation of genetic diversity and differentiation statistics with tissue-based approaches much. Finally, we confirmed that analyses using eDNA data can reveal patterns such as isolation-by-distance shown in previous studies on this species, indicating the applicability of eDNA to basic landscape genetics. Even though some limitations remain, eDNA may have great potential for conducting basic landscape genetics.

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.

Evaluation of genetic consequences of stocking on the southern-margin populations of white-spotted charr

Coldwater-adapted freshwater fishes, especially their populations along warm-range margins, are most vulnerable to the climate oscillations associated with global warming. Stocking is a major strategy for avoiding the extinction of these species. However, while stocking can reverse the decline of isolated populations, it may also result in a loss of genetic diversity in the native local population due to the introgressive replacement of hatchery genes. To plan an adequate strategy for conserving locally adapted populations, the genetic impacts of stocking on native lineages should be evaluated from small river branches to wide-ranging drainage areas. We investigated the population genetic structure of white-spotted charr (Salvelinus leucomaenis) within its southern range (Lake Biwa basin, Japan). By applying genome-wide SNP analysis to the population's genetic structure, we assessed the extent of genetic introgression resulting from stocking. White-spotted charr in the Lake Biwa watershed constitutes a distinctive genetic group, within which apparent genetic differentiation was observed. The hatchery-reared fish line commonly used for supplementation stocking in the catchment was discernable from the native population, enabling us to analyze genetic introgression across the entire drainage area. Admixed individuals resulting from hatchery introgression were observed in most of the stocked sites that showed relatively high heterozygosity and nucleotide diversity. However, their genetic differentiation was much lower than that of native populations. The supplementation history as well as the road availability contributed substantially to the introgression of hatchery genes. Populations with the native genetic structure remained in the upstream regions of the tested rivers. However, their heterozygosity and nucleotide diversity were low when compared with that of the populations with hatchery supplementation. Our results shed light on the genetic impacts of stocking on isolated native populations and suggest that conventional supplementation methods cannot preserve a unique biodiversity in the distribution margin.

Phylogeography of Parnassius citrinarius Based on Mitochondrial DNA Reveals Large Differences in Genetic Structure between the Eastern and Western Japan

The Japanese Archipelago hosts a rich butterfly fauna, and elucidating the genetic structures of multiple species is necessary to clarify their formation processes. This study aimed to reveal the genetic structure and distribution formation process of Parnassius citrinarius, which is widely distributed across the Japanese Archipelago from Hokkaido to Shikoku, through phylogeographic analysis based on the mitochondrial cytochrome c oxidase subunit I (COI) gene sequence. Thirty haplotypes were revealed from 311 individuals from 47 sites, indicating significant differences in the genetic structures between the eastern and western parts of the Japanese Archipelago. In Eastern Japan, multiple genetic clusters were found, with some sites harboring two clusters. The divergence times among populations in Eastern Japan were relatively recent, and no genetic differentiation was observed between regions, including between Hokkaido and Honshu, which are separated by a narrow strait. In contrast, in Western Japan, including Shikoku, unique genetic clusters were observed in each region. The phylogenetic relationships among populations were regionally clustered, and the divergence times were relatively ancient. The distribution and genetic structure of P. citrinarius in the Japanese Archipelago have been significantly influenced by temperature fluctuations and the presence of geographical barriers during the Pleistocene glacial-interglacial cycles, including the potential formation of refugia in Western Japan.

Male mice are susceptible to brain dysfunction induced by early-life acephate exposure

Background

Acephate is a widely used organophosphate insecticide. Exposure to endocrine-disrupting chemicals, such as acephate, can interfere with neurodevelopment in childhood, increasing the risk of higher brain dysfunction later in life. Furthermore, brain dysfunction may be related to chemical exposure-related disturbances in the gut microbiota. However, the effects of early acephate exposure on the brains of adult males and females as well as on the adult gut environment remain poorly understood.

Methods

This study investigated the effects of perinatal acephate exposure on the central nervous system and gut microbiota of mice, including sex differences and environmentally relevant concentrations. C57BL/6 N pups were exposed to acephate (0, 0.3, 10, and 300 ppm) via the dam in their drinking water from embryonic day (E) 11.5 to postnatal day 14. We examined its effects on the central nervous system of adult males and females.

Results

In the male treatment group, impairments in learning and memory were detected. Immunohistochemical analysis revealed a decrease in SOX2-, NeuN-, DCX-, and GFAP-positive cells in the hippocampal dentate gyrus in males compared to the control group, whereas GFAP-positive cells were fewer in females. In addition, gut microbiota diversity was reduced in both sexes in the experimental group.

Conclusion

Our study demonstrates that the effects of early-life exposure to acephate are more pronounced in males than in females and can lead to a lasting impact on adult behavior, even at low doses, and that the gut microbiota may reflect the brain environment.

Fine-scale clonal structure of the lingonberry Vaccinium vitis-idaea under the nurse plant Pinus pumila vegetation in an alpine region, Mt. Norikura

The nurse effect is a positive interaction in which a nurse plant improves the abiotic environment for another species (beneficiary plant) and facilitates its establishment. The evergreen shrub Vaccinium vitis-idaea (a beneficiary plant) grows mainly under the dwarf shrub Pinus pumila (a nurse plant) in the alpine regions of central Japan. However, whether V. vitis-idaea shrubs under various P. pumila shrubs spread through clonal growth and/or seeds remains unclear. We investigated the clonal structure of V. vitis-idaea under the nurse plant P. pumila in Japanese alpine regions. MIG-seq analysis was conducted to clarify the clonal diversity of V. vitis-idaea in isolated and patchy P. pumila plots on a ridge (PATs), and in a plot covered by dense P. pumila on a slope adjacent to the ridge (MAT) on Mt. Norikura, Japan. We detected 28 multilocus genotypes in 319 ramets of V. vitis-idaea across 11 PATs and MAT. Three genets expanded to more than 10 m in the MAT. Some genets were shared among neighboring PATs or among PATs and MAT. These findings suggest that the clonal growth of V. vitis-idaea plays an important role in the sustainability of populations. The clonal diversity of V. vitis-idaea was positively related with the spatial size of PATs and was higher in MAT than in PATs at a small scale. Therefore, the spatial spread of the nurse plant P. pumila might facilitate the nurse effect on the genetic diversity of beneficiary plants, leading to the sustainability of beneficiary populations.

Hybridization and reticulate evolution in Diphasiastrum (flat-branched clubmosses, Lycopodiaceae) - New data from the island of Taiwan and Vietnam

In the species groups related to Diphasiastrum multispicatum and D. veitchii, hybridization was investigated in samples from northern and southern Vietnam and the island of Taiwan, including available herbarium specimens from southeast Asia. The accessions were analyzed using flow cytometry (living material only), Sanger sequencing and multiplexed inter-simple sequence repeat genotyping by sequencing. We detected two cases of ancient hybridization involving different combinations of parental species; both led via subsequent duplication to tetraploid taxa. A cross D. multispicatum × D. veitchii from Malaysia represents D. wightianum, a tetraploid taxon according to reported DNA content measurements of dried material (genome formulas MM, VV and MMVV, respectively). The second case involves D. veitchii and an unknown diploid parent (genome formula XX). Three hybridogenous taxa (genome formulas VVX, VVXX, VVVX) were discernable by a combination of flow cytometry and molecular data. Taxon I (VVX, three clones found on Taiwan island) is apparently triploid. Taxon II represents another genetically diverse and sexual tetraploid species (VVXX) and can be assigned to D. yueshanense, described from Taiwan island but occurring as well in mainland China and Vietnam. Taxon III is as well most likely tetraploid (VVVX) and represented by at least one, more likely two, clones from Taiwan island. Taxa I and III are presumably asexual and new to science. Two independently inherited nuclear markers recombine only within, not between these hybrids, pointing towards reproductive isolation. We present an evolutionary scheme which explains the origin of the hybrids and the evolution of new and fully sexual species by hybridization and subsequent allopolyploidization in flat-branched clubmosses.

Traditional multilocus phylogeny fails to fully resolve Palearctic ground squirrels (Spermophilus) relationships but reveals a new species endemic to West Siberia

Previous efforts to reconstruct evolutionary history of Palearctic ground squirrels within the genus Spermophilus have primarily relied on a single mitochondrial marker for phylogenetic data. In this study, we present the first phylogeny with comprehensive taxon sampling of Spermophilus via a conventional multilocus approach utilizing five mitochondrial and five nuclear markers. Through application of the multispecies coalescent model, we constructed a species tree revealing four distinct clades that diverged during the Late Miocene. These clades are 1) S. alaschanicus and S. dauricus from East Asia; 2) S. musicus and S. pygmaeus from East Europe and northwestern Central Asia; 3) the subgenus Colobotis found across Central Asia and its adjacent regions and encompassing S. brevicauda, S. erythrogenys, S. fulvus, S. major, S. pallidicauda, S. ralli, S. relictus, S. selevini, and S. vorontsovi sp. nov.; and 4) a Central/Eastern Europe and Asia Minor clade comprising S. citellus, S. taurensis, S. xanthoprymnus, S. suslicus, and S. odessanus. The latter clade lacked strong support owing to uncertainty of taxonomic placement of S. odessanus and S. suslicus. Resolving relationships within the subgenus Colobotis, which radiated rapidly, remains challenging likely because of incomplete lineage sorting and introgressive hybridization. Most of modern Spermophilus species diversified during the Early-Middle Pleistocene (2.2-1.0 million years ago). We propose a revised taxonomic classification for the genus Spermophilus by recognizing 18 species including a newly identified one (S. vorontsovi sp. nov.), which is found only in a limited area in the southeast of West Siberia. Employing genome-wide single-nucleotide polymorphism genotyping, we substantiated the role of the Ob River as a major barrier ensuring robust isolation of this taxon from S. erythrogenys. Despite its inherent limitations, the traditional multilocus approach remains a valuable tool for resolving relationships and can provide important insights into otherwise poorly understood groups. It is imperative to recognize that additional efforts are needed to definitively determine phylogenetic relationships between certain species of Palearctic ground squirrels.

Degenerate oligonucleotide primer MIG-seq: an effective PCR-based method for high-throughput genotyping

Next-generation sequencing (NGS) library construction often involves using restriction enzymes to decrease genome complexity, enabling versatile polymorphism detection in plants. However, plant leaves frequently contain impurities, such as polyphenols, necessitating DNA purification before enzymatic reactions. To overcome this problem, we developed a PCR-based method for expeditious NGS library preparation, offering flexibility in number of detected polymorphisms. By substituting a segment of the simple sequence repeat sequence in the MIG-seq primer set (MIG-seq being a PCR method enabling library construction with low-quality DNA) with degenerate oligonucleotides, we introduced variability in detectable polymorphisms across various crops. This innovation, named degenerate oligonucleotide primer MIG-seq (dpMIG-seq), enabled a streamlined protocol for constructing dpMIG-seq libraries from unpurified DNA, which was implemented stably in several crop species, including fruit trees. Furthermore, dpMIG-seq facilitated efficient lineage selection in wheat and enabled linkage map construction and quantitative trait loci analysis in tomato, rice, and soybean without necessitating DNA concentration adjustments. These findings underscore the potential of the dpMIG-seq protocol for advancing genetic analyses across diverse plant species.

Contrasting population genetic structure of three semi-terrestrial brachyuran crabs on the coast of the Japanese archipelago

Anthropogenic activities have reduced ecotones between the ocean and land, which is likely to threaten the population of brackish-water brachyuran crabs. To assess the current status of these crabs, we examine the population genetic structures of three semi-terrestrial brachyuran crabs widely distributed along the coast of the Japan and to clarify factors determining their genetic structures. We collected 184 Orisarma dehaani, 252 Chiromantes haematocheir, and 151 Helice tridens crabs from 36 localities of the Japanese archipelago. Genome-wide SNP data from these crabs were analyzed using MIG-seq. Bayesian clustering of STRUCTURE and DAPC analysis were used to identify genetically disturbed populations and to visualize genetic differentiation between local populations. Genetic population structure showed clear differentiation between populations on the Pacific coast of the Tohoku region and on other Japanese coasts in O. dehaani, but not in C. haematocheir or H. tridens. The inbreeding coefficient of O. dehaani was significantly higher on the Pacific coast of the Tohoku region compared to other Japanese coasts. C. haematocheir and H. tridens had homogeneous genetic structures along the Japanese coast, but showed genetic differentiation of a local population at their range limits. Thus, O. dehaani showed little gene flow and clear genetic differentiation between populations in the Tohoku Pacific region and those on other Japanese coasts due to ocean currents. Although such a regional differentiation was not found in C. haematocheir and H. tridens, one population of C. haematocheir was genetically isolated at the edge of its distribution range and likely vulnerable to environmental changes.

Intraspecific divergence in a coastal plant, Euphorbia jolkinii, at a major biogeographic boundary in East Asia

PREMISE

Quaternary climatic fluctuations and long-distance seed dispersal across the sea are critical factors affecting the distribution of coastal plants, but the spatiotemporal nature of population expansion and distribution change of East Asian coastal plants during this period are rarely examined. To explore this process, we investigated the genome-wide phylogenetic patterns of Euphorbia jolkinii Boiss. (Euphorbiaceae), which grows widely on littoral areas of Japan, Korea, and Taiwan.

METHODS

We used plastome sequences and genome-wide single nucleotide polymorphisms in samples across the species range to reveal phylogeographic patterns and spatiotemporal distributional changes. We conducted ecological niche modeling for the present and the last glacial maximum (LGM).

RESULTS

Genetic differentiation was observed between the northern and southern populations of E. jolkinii, separated by the major biogeographic boundary, the Tokara Gap. These two groups of populations differentiated during the glacial period and subsequently intermingled in the intermorainic areas of the central Ryukyu Islands after the LGM. Ecological niche models suggested that the potential range of E. jolkinii was restricted to southern Kyushu; however, it was widespread in the southern Ryukyu Islands and Taiwan during the LGM.

CONCLUSIONS

This study provides evidence of genetic differentiation among coastal plant populations separated by the prominent biogeographical boundary. Although coastal plants are typically expected to maintain population connectivity through sea-drifted seed dispersal, our findings suggest that genetic differences may arise because of a combination of limited gene flow and changes in climate during the glacial period.

Wasabi Gone Wild? Origin and Characterization of the Complete Plastomes of Ulleung Island Wasabi (Eutrema japonicum; Brassicaceae) and Other Cultivars in Korea

Korean wasabi occurs naturally on the young oceanic, volcanic Ulleung Island off the east coast of the Korean Peninsula. Although the Ulleung Island wasabi is reported as Eutrema japonicum and has been suggested to be morphologically identical to cultivars in Korea, very little is known about its taxonomic identity and relationship with other cultivars. In this study, we sequenced the complete chloroplast DNA sequences of three naturally occurring Ulleung Island wasabi plants and six cultivars ('Daewang', 'Daruma', 'Micado', 'Orochi', 'Green Thumb', and 'Shogun') from continental Korea and determined the taxonomic identity of Korean wasabi on Ulleung Island. The size and organization of the complete chloroplast genomes of the nine accessions were nearly identical to those of previously reported wasabi cultivars. In addition, phylogenetic analysis based on the complete plastomes suggested that Ulleung Island wasabi most likely comprises various wasabi cultivars with three chlorotypes ('Shogun', 'Green Thumb', and a unique Chusan type). Based on the complete plastomes, we identified eight chlorotypes for the major wasabi cultivars and the Ulleung Island wasabi. Two major groups (1-'Mazuma' and 'Daruma', and 2-'Fujidaruma'/'Shimane No. 3'/Ulleung Island wasabi/five cultivars in Korea) were also identified based on mother line genealogical history. Furthermore, different types of variations (mutations, insertions/deletions (indels), mononucleotide repeats, and inversions) in plastomes were identified to distinguish different cultivar lines and five highly divergent hotspots. The nine newly obtained complete plastomes are valuable organelle genomic resources for species identification and infraspecific phylogeographic studies on wild and cultivated wasabi.

Study on the relationship between the dispersal of wild boar (Sus scrofa) and the associated variability of Cesium-137 concentrations in its muscle Post-Fukushima Daiichi Nuclear Power Plant accident

After the Fukushima Daiichi Nuclear Power Plant (FDNPP) accident in 2011, the wild boar (Sus scrofa) population within the Fukushima Evacuation Zone (FEZ) increased substantially in size and distribution. This growing population and their potential dispersal from the FEZ, where they are exposed to high levels of radionuclides, into the surrounding landscape underscores the need to better understand boar movement patterns in order to establish policies for managing shipping restrictions for boar meat and develop management strategies. In this study, we quantified the genetic population structure of boar in and around Fukushima prefecture using sequence data of the mitochondrial DNA control region and MIG-seq analysis using 348 boar samples to clarify boar dispersal patterns. Among boar samples, seven Asian haplotypes and one European haplotype were detected. The European haplotype originated from hybridization between domestic pigs and native boar in the evacuation zone after the accident and was detected in 15 samples across a broad geographic area. Our MIG-seq analysis revealed genetic structure of boar was significantly different between boar inhabiting the eastern (including FEZ. i.e., East clade) and western (i.e., West clade) regions in Fukushima prefecture. In addition, we investigated the relationships between boar dispersal and Cesium (Cs)-137 activity concentrations in boar muscle using MIG-seq genetic data in Nihonmatsu city, located in the central-northern region of Fukushima. High Cs-137 activity concentrations, exceeding 1000 Bq/kg, in boar muscle had a significantly high probability of belonging to the East clade within localized regions. Thus, our results provide evidence of the spatial scale of dispersal of individuals or offspring of boar from the FEZ. Results of this research also indicate that dispersal of individuals between areas with different Cs-137 contamination levels is one of the biggest factors contributing to variation in Cs-137 activity concentration in boar muscle within localized regions.

Population genetics of Sida fallax Walp. (Malvaceae) in the Hawaiian Islands

Introduction

Sida fallax (Malvaceae) is the most widespread and variable taxon of Malvaceae in the Hawaiian Islands, growing with a diversity of morphological forms in different habitats including Midway Atoll, Nihoa, and all the main islands. Morphological variation exists within and among populations. The study aimed to investigate the genetic variation within and among populations from various habitats and geographic locations throughout the Hawaiian range of S. fallax.

Methods

A total of 124 samples, with up to five samples per population where possible, were collected from 26 populations across six of the main Hawaiian Islands (Kaua'i, O'ahu, Maui, Moloka'i, Lāna'i, and Hawai'i) and Nihoa in the Northwestern Hawaiian Islands. The sampling strategy encompassed collecting populations from different habitats and geographic locations, including coastal and mountain ecotypes, with many intermediate morphological forms. Multiplexed ISSR genotyping by sequencing (MIG-seq) was used to detect single nucleotide polymorphisms (SNP) and genetic differences among individuals and populations were evaluated using PCO analyses.

Results

The relationship of FST with the geographical distance between the populations was assessed using the Mantel test. The results showed that populations on a single island were more closely related to each other and to populations on islands within their respective groups than they were to populations on other islands.

Discussion

The overall genetic relationships among islands were, to a large extent, predictive based on island position within the chain and, to a lesser extent, within island topography.

Natural habitat connectivity and organic management modulate pest dispersal, gene flow, and natural enemy communities

The simplification and fragmentation of agricultural landscapes generate effects on insects at multiple spatial scales. As each functional group perceives and uses the habitat differently, the response of pest insects and their associated natural enemies to environmental changes varies. Therefore, landscape structure may have consequences on gene flow among pest populations in space. This study aimed to evaluate the effects of local and landscape factors, at multiple scales, on the local infestation, gene flow and broad dispersion dynamics of the pest insect Bemisia tabaci (Genn.) Middle East-Asia Minor 1 (MEAM-1, former biotype B) (Hemiptera: Aleyrodidae) and its associated natural enemies in a tropical agroecosystem. We evaluated the abundance of B. tabaci populations and their natural enemy community in 20 tomato farms in Brazil and the gene flow between farms from 2019 to 2021. Landscapes dominated by agriculture resulted in larger B. tabaci populations and higher gene flow, especially in conventional farms. A higher density of native vegetation patches disfavored pest populations, regardless of the management system. The results revealed that whitefly responds to intermediate spatial scales and that landscape factors interact with management systems to modulate whitefly populations on focal farms. Conversely, whitefly natural enemies benefited from higher amounts of natural vegetation at small spatial scales, while the connectivity between natural habitat patches was beneficial for natural enemies regardless of the distance from the focal farm. The resulting dispersion model predicts that the movement of whiteflies between farms increases as the amount of natural vegetation decreases. Our findings demonstrate that landscape features, notably landscape configuration, can mediate infestation episodes, as they affect pest insects and natural enemies in opposite ways. We also showed that landscape features interact with farm traits, which highlights the need for management strategies at multiple spatial scales. In conclusion, we demonstrated the importance of the conservation of natural areas as a key strategy for area-wide ecological pest management and the relevance of organic farming to benefit natural enemy communities in tropical agroecosystems.

MCPtaggR: R package for accurate genotype calling in reduced representation sequencing data by eliminating error-prone markers based on genome comparison

Reduced representation sequencing (RRS) offers cost-effective, high-throughput genotyping platforms such as genotyping-by-sequencing (GBS). RRS reads are typically mapped onto a reference genome. However, mapping reads harbouring mismatches against the reference can potentially result in mismapping and biased mapping, leading to the detection of error-prone markers that provide incorrect genotype information. We established a genotype-calling pipeline named mappable collinear polymorphic tag genotyping (MCPtagg) to achieve accurate genotyping by eliminating error-prone markers. MCPtagg was designed for the RRS-based genotyping of a population derived from a biparental cross. The MCPtagg pipeline filters out error-prone markers prior to genotype calling based on marker collinearity information obtained by comparing the genome sequences of the parents of a population to be genotyped. A performance evaluation on real GBS data from a rice F2 population confirmed its effectiveness. Furthermore, our performance test using a genome assembly that was obtained by genome sequence polishing on an available genome assembly suggests that our pipeline performs well with converted genomes, rather than necessitating de novo assembly. This demonstrates its flexibility and scalability. The R package, MCPtaggR, was developed to provide functions for the pipeline and is available at https://github.com/tomoyukif/MCPtaggR.

Habitat diversification associated with urban development has a little effect on genetic structure in the annual native plant Commelina communis in an East Asian megacity

Urban development greatly alters the natural and semi-natural habitats of native plants. Urbanisation results in a range of diverse habitats including remnant agricultural lands, urban parks, and roadside habitats. This habitat diversity often promotes trait divergence within urban areas. However, the mechanisms by which diverse urban habitats influence the population genetic structure of individual plant species remain poorly understood. We investigated the effects of urbanisation on genetic diversity and structure within 24 Commelina communis populations across diverse habitat types (rural agricultural land, urban agricultural land, urban park land, and urban roadsides) within the Kyoto-Osaka-Kobe megacity in Japan. We conducted multiplexed inter-simple sequence repeat genotyping to compare genetic diversity among populations in different habitats. We also examined the correlation between Nei's genetic distance and geographic and environmental distances and performed principal coordinate analysis (PCoA) to evaluate genetic differentiation among urban habitats. There were no significant differences in genetic diversity indices between urban and rural populations and among urban habitat types. Although we detected no isolation-by-distance structure in population pairs of the same habitat type and in those of different habitats, the difference in surrounding landscape facilitated genetic differentiation not only between urban and rural habitats but also between different urban habitats. PCoA revealed no clear genetic differentiation among rural and urban habitat populations. Our findings indicate that the establishment of diverse habitat types through urbanisation has no and little impact on genetic diversity and structure, respectively, in C. communis, likely due to its high selfing rate and ability to adapt to urban conditions.

Double migration of the endangered Tricyrtis formosana (Liliaceae) in Japan

The Ryukyu Islands of Japan are a biodiversity hotspot due to geographical and historical factors. Tricyrtis formosana is a perennial herbaceous plant that commonly found in Taiwan. But only a few populations have been identified in a limited habitat on Iriomote Island, while populations of unknown origin occur near human settlements in an area on the main island of Okinawa. To better understand these populations of the phylogenetic uniqueness and intrinsic vulnerability, we conducted comparative analyses including (1) phylogeny and population structure with MIG-seq data, (2) photosynthesis-related traits of plants grown under common conditions and (3) transcriptome analysis to detect deleterious variations. Results revealed that T. formosana was split into two clades by the congeners and that Iriomote and Okinawa populations independently derived from ancestral Taiwanese populations in each clade. Photosynthetic efficiency was lowest in the Iriomote population, followed by Okinawa and Taiwan. Transcriptome analysis showed that the Iriomote population accumulated more deleterious variations, suggesting intrinsic vulnerability. These results indicate that each T. formosana population in Japan is phylogenetically unique and has been independently dispersed from Taiwan, and that the Iriomote population presents a high conservation difficulty with a unique photosynthesis-related characteristic and a larger amount of deleterious variations.

The phylogenetic reconstruction of the Neotropical cycad genus Ceratozamia (Zamiaceae) reveals disparate patterns of niche evolution

The cycad genus Ceratozamia comprises 40 species from Mexico, Guatemala, Belize, and Honduras, where cycads occur throughout climatically varied montane habitats. Ceratozamia has the potential to reveal the history and processes of species diversification across diverse Neotropical habitats in this region. However, the species relationships within Ceratozamia and the ecological trends during its evolution remain unclear. Here, we aimed to clarify the phylogenetic relationships, the timing of clade and species divergences, and the niche evolution throughout the phylogenetic history of Ceratozamia. Genome-wide DNA sequences were obtained with MIG-seq, and multiple data-filtering steps were used to optimize the dataset used to construct an ultrametric species tree. Divergence times among branches and ancestral niches were estimated. The niche variation among species was evaluated, summarized into two principal components, and their ancestral states were reconstructed to test whether niche shifts among branches can be explained by random processes, under a Brownian Motion model. Ceratozamia comprises three main clades, and most species relationships within the clades were resolved. Ceratozamia has diversified since the Oligocene, with major branching events occurring during the Miocene. This timing is consistent with fossil evidence, the timing estimated for other Neotropical plant groups, and the major geological events that shaped the topographic and climatic variation in Mexico. Patterns of niche evolution in the genus do not accord with the Brownian Motion model. Rather, non-random evolution with shifts towards more seasonal environments at high latitudes, or shifts towards humid or dry environments at low latitudes explain the diversification of Ceratozamia. We present a comprehensive phylogenetic reconstruction for Ceratozamia and identify for the first time the environmental factors involved in clade and species diversification within the genus. This study alleviates the controversies regarding the species relationships in the genus and provides the first evidence that latitude-associated environmental factors may influence processes of niche evolution in cycads.

Phylogenetic, population structure, and population demographic analyses reveal that Vicia sepium in Japan is native and not introduced

Vicia sepium (bush vetch) is a perennial legume widely distributed throughout the Eurasian continent. However, its distribution in Japan is limited to Mt. Ibuki and small parts of central and southern Hokkaido. Therefore, each Japanese V. sepium lineage has been considered to have been introduced separately from Europe. Here, we examined whether the species was introduced or not on the basis of cpDNA sequences and genome-wide SNPs from Japanese and overseas samples. Both the cpDNA haplotype network and the nuclear DNA phylogenetic tree showed that Japanese V. sepium is monophyletic. Furthermore, although the nuclear DNA phylogenetic tree also showed that each lineage is clearly monophyletic, genetic admixture of the genetic cluster dominated in the Hokkaido lineage was also detected in the Mt. Ibuki lineage. Population divergence analysis showed that the two lineages diverged during the last glacial period. The Mt. Ibuki lineage showed a sudden population decline 300-400 years ago, indicating that some anthropogenic activity might be involved, while the Hokkaido lineage showed a gradual population decline from 5000 years ago. Consequently, these two lineages show low current genetic diversity compared with overseas lineages. These results show that the Japanese V. sepium is not introduced but is native.

Molecular phylogeny and taxonomy of Hosta (Asparagaceae) on Shikoku Island, Japan, including five new species, one new subspecies, and two new status assignments

Japan has 16 native species of the genus Hosta Tratt. (Asparagaceae). A recent study on Hosta based on field surveys and molecular phylogenetic analyses resulted in the discovery of six unknown taxa in Kochi Prefecture, Shikoku Island, southwestern Japan. We aimed to identify these unknown taxa. Therefore, we constructed a finely resolved phylogeny for 320 Hosta samples collected from the Honshu, Shikoku, and Kyushu Islands using multiplex inter-simple sequence repeat genotyping by sequencing (MIG-seq). Based on this phylogenetic analysis and related morphological observations, we describe five new species, H.longipedicellatasp. nov., H.minazukiflorasp. nov., H.polyneuronoidessp. nov., H.samukazemontanasp. nov., and H.takiminazukiflorasp. nov. and one new subspecies, H.takiminazukiflorasubsp.grandissubsp. nov. In addition, we propose two new status assignments, H.tardivasubsp.densinerviacomb. and stat. nov. and H.scabrinerviastat. nov. We also propose classifying H.kikutiivar.tosana as a species, H.tosana. Further studies that combine MIG-seq with careful morphological observations are needed for Hosta plants on all Japanese islands, which may result in the discovery of even more undescribed species.

Scrutinising an inscrutable bark-nesting ant: Exploring cryptic diversity in the Rhopalomastix javana (Hymenoptera: Formicidae) complex using DNA barcodes, genome-wide MIG-seq and geometric morphometrics

Overlooking cryptic species diversity has grave implications on assessments of climate change impacts on biodiversity, ecosystems and organismal populations. Discriminating between cryptic species has long been challenging even for seasoned taxonomists, as interspecies morphological differences are often indiscernible by visual observation. Multi-disciplinary methods involving genetic analyses in conjunction with quantitative morphological data, should therefore be used to investigate boundaries between cryptic species. We adopted an integrated approach combining analyses of mitochondrial COI barcodes, a genome-wide dataset obtained via multiplexed inter-simple sequence repeats (ISSRs) genotyping by sequencing (MIG-seq), and geometric morphometrics to investigate species divergences in the inscrutable Rhopalomastix javana species complex. Objective clustering of COI suggested five putative molecular species units divergent from each other by thresholds within 4.2-10.6% uncorrected pairwise distance. Phylogenetic analyses based on concatenated MIG-seq data also recovered and strongly supported the monophyly of five major lineages in agreement with COI clusters. Co-ancestry analyses based on MIG-seq data using fineRADstructure resolved variable patterns of admixture linked to geography, and potential genetic drift within some putative species. Geometric morphometric analyses of specimen images further detected statistically significant differences in at least one of three anatomical aspects (Head, Meso, Profile) between all pairs of putative species. Head shape (full-face view) was determined to be the most informative character for species diagnosis, with relatively high classification accuracy. Thin-plate spline deformation grids highlighted areas of high variation between species in each shape for deeper taxonomic scrutiny. The presence of species from multiple distinct lineages existing in near-sympatry firmly demonstrates that R. javana comprises more than one closely-related species, but exact species boundaries are difficult to ascertain. Differences in elevation and its associated abiotic effects on ant adaptations and reproductive phenology may contribute to restricting gene flow and maintaining species boundaries between sympatric populations of the R. javana complex. We further assess the advantages and limitations of geometric morphometrics as a taxonomic tool. Despite its drawbacks, our combined approach has helped draw important insights on cryptic diversity in R. javana, and also identified gaps of knowledge that await address. Results from this study will inform and prime future in-depth taxonomic investigation on the R. javana complex, including formal descriptions and establishment of the five putative species.

Impact of mycoheterotrophy on the growth of Gentiana zollingeri (Gentianaceae), as suggested by size variation, morphology, and 13C abundance of flowering shoots

Gentiana zollingeri is an annual photosynthetic plant that employs a mycoheterotrophic growth strategy during its underground seedling stage (initial mycoheterotrophy). Notably, the morphological characteristics of its flowering shoots, such as shoot size, leaf size, and leaf color, are highly variable, and it was hypothesized that these variations may be linked to nutritional mode. The morphological characteristics of G. zollingeri individuals were thus investigated alongside environmental factors, 13C abundance, and diversity of colonizing arbuscular mycorrhizal (AM) fungi. The majority of G. zollingeri flowering individuals were found to exhibit a high affinity for the specific AM fungi that exclusively colonize roots of the mycoheterotrophic seedlings, while other phylogenetically diverse AM fungi could also be detected. The leaves to shoot dry weight ratio (leaf ratio) was negatively correlated with the canopy openness in the habitat, suggesting that leaf development is impeded in sunny conditions. Furthermore, the shoot weight of G. zollingeri was positively correlated with leaf 13C abundance. Given that 13C enrichment can provide indirect evidence of mycoheterotrophy in AM plants, the results suggest that the utilization of carbon obtained through mycoheterotrophy, at least during the underground seedling stage, is crucial for G. zollingeri.

Host-Shift Speciation Proceeded with Gene Flow in Algae Covering Shells

AbstractHost shifts represent the advancement of a novel niche and often lead to speciation in symbionts. However, its mechanisms are not well understood. Here, we focused on the alga Pseudocladophora conchopheria growing on the shells of intertidal snails. Previous surveys have shown that the alga has host specificity-only attaching to the shell of Lunella correensis-but we discovered that the alga attaches to the shells of multiple sympatric snails. A genome-wide single-nucleotide polymorphism analysis (MIG-seq) was performed to determine whether host-associated speciation occurred in the algae. As a result, there was no gene flow or limited gene flow among the algae from different hosts, and some algae were genetically differentiated among hosts. In addition, the demographic estimate revealed that speciation with gene flow occurred between the algae from different hosts. Therefore, these results support the idea that host-shift speciation gradually proceeded with gene flow in the algae, providing insight into the early evolution of host shifts.

Population genetic structure and hybrid zone analyses for species delimitation in the Japanese toad (Bufo japonicus)

Hybridization following secondary contact may produce different outcomes depending on the extent to which genetic diversity and reproductive barriers have accumulated during isolation. The Japanese toad, Bufo japonicus, is distributed on the main islands of Japan. In the present study, we applied multiplexed inter-simple sequence repeat genotyping by sequencing to achieve the fine-scale resolution of the genetic cluster in B. j. japonicus and B. j. formosus. We also elucidated hybridization patterns and gene flow degrees across contact zones between the clusters identified. Using SNP data, we found four genetic clusters in B. j. japonicus and B. j. formosus and three contact zones of the cluster pairs among these four clusters. The two oldest diverged lineages, B. j. japonicus and B. j. formosus, formed a narrow contact zone consistent with species distinctiveness. Therefore, we recommend that these two subspecies be elevated to the species level. In contrast, the less diverged pairs of two clusters in B. j. japonicus and B. j. formosus, respectively, admixed over a hundred kilometers, suggesting that they have not yet developed strong reproductive isolation and need to be treated as conspecifics. These results will contribute to resolving taxonomic confusion in Japanese toads.

Discrimination of Camellia cultivars using iD-NA analysis

Recently, many new cultivars have been taken abroad illegally, which is now considered an international issue. Botanical evidence found at a crime scene provides valuable information about the origin of the sample. However, botanical resources for forensic evidence remain underutilized because molecular markers, such as microsatellites, are not available without a limited set of species. Multiplexed intersimple sequence repeat (ISSR) genotyping by sequencing (MIG-seq) and its analysis method, identification of not applicable (iD-NA), have been used to determine several genome-wide genetic markers, making them applicable to all plant species, including those with limited available genetic information. Camellia cultivars are popular worldwide and are often planted in many gardens and bred to make new cultivars. In this study, we aimed to analyze Camellia cultivars/species through MIG-seq. MIG-seq could discriminate similar samples, such as bud mutants and closely related samples that could not be distinguished based on morphological features. This discrimination was consistent with that of a previous study that classified cultivars based on short tandem repeat (STR) markers, indicating that MIG-seq has the same or higher discrimination ability as STR markers. Furthermore, we observed unknown phylogenetic relationships. Because MIG-seq can be applied to unlimited species and low-quality DNA, it may be useful in various scientific fields.

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.

Population genetic structure of wild Angelica acutiloba, A. acutiloba var. iwatensis, and their hybrids by atpF-atpA intergenic spacer in chloroplast DNA and genome-wide SNP analysis using MIG-seq

Sampling surveys of Angelica acutiloba and A. acutiloba var. iwatensis, which are medicinal plants endemic to Japan, were conducted in the Chubu region in the central area of the main island of Japan. A. acutiloba grows in riverbeds in mountainous areas, while A. acutiloba. var. iwatensis grows on slopes near mountain ridges at 1000 m above sea level or on constantly collapsing rocky slopes and bare fields on developed land along asphalt roads in valleys of mountainous areas. Specimens of two wild Angelica species collected in this region were examined for maternal lineage by DNA polymorphism analysis of the atpF-atpA region for chloroplast DNA using direct sequencing and genomic component analysis by genome-wide SNP using MIG-seq. In this study area, while all A. acutiloba populations were monophyletic in both maternal and ancestral lineages, A. acutiloba var. iwatensis were genetically heterogeneous due to being composed of three maternal and three ancestral lineages to various degrees. In addition, a natural hybrid population with maternal lineage presumed to be A. acutiloba and paternal lineage A. acutiloba var. iwatensis was also found. In the present study, we report that the combined method of atpF-atpA and MIG-seq analyses is a useful tool for determining the population genetic structure of two wild Angelica species and for identifying hybrids.

Taxonomic reassessment of salamanders (genus Hynobius) from Tsushima Islands, Japan, with a resurrection of Hynobius tagoi Dunn, 1923 (Amphibia: Caudata)

To clarify the taxonomic status of salamanders from Tsushima Islands, Japan, we examined two species of salamanders, Hynobius tsuensis and Hynobius sp. from Tsushima Islands, and compared them with H. nebulosus from Kyushu. We found that the three taxa differ from each other in nuclear DNA, adult morphology, and egg-sac shape, and consider them to be independent species. Based on the investigation of type specimens of the synonyms, we identified Hynobius sp. as H. tagoi Dunn, 1923 and redescribe it herein. Also, we designate a neotype for H. tsuensis and redescribe H. tsuensis. The distribution areas of the two species partly overlap but they were distinguishable by their body colorations: H. tsuensis has uniformly blackish tail sides and a distinct yellow stripe on the upper tail edge; whereas H. tagoi has brown tail sides with numerous dark stippling and without a distinct yellow stripe on the upper tail edge.

The genetic structure and demographic history of Zabelia tyaihyonii, endemic to Korean limestone karst forests, based on genome-wide SNP markers

Similar to the global phenomenon, many plant species endemic to Korean limestone karst forests are at risk of extinction due to human intervention. Zabelia tyaihyonii is a familiar shrub, called "Hardy abelia" and "Fragrant abelia" growing in the karst forests of Korea, where it is one of the most threatened species. We investigated the genetic structure and demographic history of Z. tyaihyonii, which allow us to develop appropriate conservation and management strategies. The genetic structure was evaluated using a total of 187 samples from 14 populations, covering the entire distribution of Z. tyaihyonii in South Korea. We utilized 254 and 1753 SNP loci obtained via MIG-seq (Multiplexed ISSR Genotyping by sequencing) for structure and demographic analyses, respectively. The population demographic modeling was performed with site frequency spectrum. To gain further historical insights, we also employed ENM (Ecological Niche Modeling). We found two distinct clusters (CLI and CLII) of ancient origin (ca. 490 ka). Despite CLII experiencing a more severe bottleneck, both clusters showed similar levels of genetic diversity, indicating mutual historical gene flow. Their historical distribution range seems to have changed very little. We proposed a historical distribution scenario for Z. tyaihyonii, taking into account its intrinsic factors, and emphasized a more complex response to Quaternary climate change beyond simple allopatric speciation models. These findings provide valuable insights for conservation and management strategies for Z. tyaihyonii.

Rapid and easy construction of a simplified amplicon sequencing (simplified AmpSeq) library for marker-assisted selection

Marker-assisted selection (MAS) is fundamental for plant breeding programs, as it can identify desirable seedlings at a young stage and reduce the cost, time and space needed for plant maintenance, especially for perennial crops. To facilitate the process of genotyping, which is time consuming and laborious, we developed a simplified amplicon sequencing (simplified AmpSeq) library construction method for next-generation sequencing that can be applied to MAS in breeding programs. The method is based on one-step PCR with a mixture of two primer sets: the first consisting of tailed target primers, the second of primers that contain flow-cell binding sites, indexes and tail sequences complementary to those in the first set. To demonstrate the process of MAS using s implified AmpSeq, we created databases of genotypes for important traits by using cultivar collections including triploid cultivars and segregating seedlings of Japanese pear (Pyrus pyrifolia Nakai), Japanese chestnut (Castanea crenata Sieb. et Zucc.) and apple (Malus domestica Borkh.). Simplified AmpSeq has the advantages of high repeatability, ability to estimate allele number in polyploid species and semi-automatic evaluation using target allele frequencies. Because this method provides high flexibility for designing primer sets and targeting any variant, it will be useful for plant breeding programs.

Reassessment of species delimitation using nuclear markers in three lentic-breeding salamanders from the Chugoku District of Japan (Amphibia: Caudata: Hynobiidae)

Hynobius akiensis sensu lato has recently been split into three species based on short sequence analyses of cyt-b gene of mtDNA and without data of nuclear DNA, and strange sympatric distribution in some areas has been indicated in two species. We analyzed nuclear DNA marker (SNPs) and complete sequence of cyt-b in H. akiensis sensu lato to reassess species delimitation and genetic introgression among species. As a result, we found two lineages with discordant mitochondrial and nuclear DNA in some areas. Of H. akiensis sensu lato, each of the two contains the type locality of two species recently reported (H. sumidai and H. geiyoensis), and the use of these names has been previously advocated. However, their sympatric distribution was rejected based on nuclear DNA data, which we consider is more reliable than mtDNA. We thus clarify geographic boundary of these two species and revise the species delimitations.

Phylogeography of the temperate grassland plant Tephroseris kirilowii (Asteraceae) inferred from multiplexed inter-simple sequence repeat genotyping by sequencing (MIG-seq) data

A group of temperate grassland plant species termed the "Mansen elements" occurs in Japan and is widely distributed in the grasslands of continental East Asia. It has been hypothesized that these species are continental grassland relicts in Japan that stretch back to a colder age, but their migration history has not been elucidated. To assess the migration history of the Mansen elements, we performed phylogeographic analyses of Tephroseris kirilowii, a member of this group, using single-nucleotide polymorphisms (SNPs) obtained from multiplexed inter-simple sequence repeat genotyping by sequencing (MIG-seq). It was estimated that the Japanese populations of T. kirilowii were divided from those of continental East Asia at 25.2 thousand years ago (ka) with 95% highest probability density interval (HPD) of 15.3-40.0 ka and that Japanese clades first diverged at 20.2 ka with 95% HPD of 10.4-30.1 ka. As the climatically suitable range during the last glacial maximum (LGM) estimated using ecological niche modeling (ENM) was limited in Japan and there was a slight genetic differentiation among Japanese populations, a post-glacial expansion of T. kirilowii in the Japanese Archipelago was indicated.

Phylogenomic analyses reveal incongruences between divergence times and fossil records of freshwater snails in East Asia

Fossils provide important insight into our understanding of phylogenetic history by serving as calibration points for divergence time estimation. However, uncertainties in the fossil record due to parallel evolution and convergent evolution can critically affect estimates of node ages. Here, we compare and contrast estimates of phylogenetic divergence with geologic and fossil history for two freshwater snail genera of the family Viviparidae in East Asia (Cipangopaludina and Margarya). Cipangopaludina species are commonly widely distributed species in East Asia, but extant Margarya species are endemic to the ancient lakes in Yunnan, China. According to some previous studies, parallel evolution or convergent evolution of shell morphology has occurred in the family several times which may affect divergence time estimation using fossil records. In this study, we used SNP data derived from ddRAD-seq loci to investigate population demographic history of both genera. Our results show a common pattern of lake endemic lineages diversifying from widely distributed lineages in the Miocene, and multiple colonization to a single ancient lake occurred in the Pleistocene. Our results indicate substantial incongruence among estimated phylogenomic divergence times, some fossil records, and formation ages of ancient lakes. These findings suggest some fossil records may be misidentified in these groups and highlight the need to carefully evaluate geological evidence and fossil records when using these for divergence time estimation.

Taxonomic study of Scytosiphon (Phaeophyceae) from temperate coasts of Argentina

Scytosiphon is a common intertidal genus widely distributed on temperate coasts worldwide. Recently, eight species have been delimited with molecular tools. Although S. lomentaria is the only species that predominates in the macroalgal literature of the Southwestern Atlantic Ocean (SwAO), unpublished molecular data obtained for a population study of S. lomentaria revealed hidden species diversity of Scytosiphon among the individuals collected from four localities at the SwAO. The aim of this study was to revise the identity and phylogenetic relationships of Scytosiphon from temperate coasts of the SwAO using DNA data. Thalli were collected from the Argentinean coast between 39° S and 43° S, from which cox1 and rbcL gene sequences were obtained. Phylogenies and haplotype networks were inferred and morphology of gametophytes was studied. Four species were recognized, S. lomentaria, S. promiscuus, S. shibazakiorum, and one species that belongs to a complex of species known as "Scytosiphon Atlantic complex." This complex was known to occur only in the North Atlantic, however, the results found in this study revealed that it has an extended distribution range that includes the southern hemisphere, where its populations have high genetic diversity and unique haplotypes. The morphological differences among the four species were subtle; denoting that previous Scytosiphon records from the SwAO attributed to the renowned S. lomentaria could represent different species. In addition, sex ratio and genome-wide single nucleotide polymorphisms (SNPs) analyses were done for populations of S. promiscuus presumably introduced to the SwAO, and the results indicated that they included female-dominant parthenogenetic populations, which were probably introduced from Japan.

A lack of population structure characterizes the invasive Lonicera japonica in West Virginia and across eastern North America

Invasive plant species cause massive ecosystem damage globally, yet represent powerful case studies in population genetics and rapid adaptation to new habitats. The availability of digitized herbarium collections data, and the ubiquity of invasive species across the landscape make them highly accessible for studies of invasion history and population dynamics associated with their introduction, establishment, spread, and ecological interactions. Here we focus on Lonicera japonica, one of the most damaging invasive vine species in North America. We leveraged digitized collections data and contemporary field collections to reconstruct the invasion history and characterize patterns of genomic variation in the eastern USA, using a straightforward method for generating nucleotide polymorphism data and a recently published, chromosome-level genome for the species. We found an overall lack of population structure among sites in northern West Virginia, USA, as well as across sites in the central and eastern USA. Heterozygosity and population differentiation were both low based on Fst, analysis of molecular variance, principal components analysis, and cluster-based analyses. We also found evidence of high inbreeding coefficients and significant linkage disequilibrium, in line with the ability of this otherwise outcrossing, perennial species to propagate vegetatively. Our findings corroborate earlier studies based on allozyme data, and suggest that intentional, human-assisted spread explains the lack of population structure, as this species was planted for erosion control and as an ornamental, escaping cultivation repeatedly across the USA. Finally, we discuss how plant invasion genomics can be incorporated into experiential undergraduate education as a way to integrate teaching and research.

GRAS-Di SNP-based molecular characterization and fingerprinting of a Turkish Corylus avellana core set provide insights into the cultivation and breeding of hazelnut in Turkey

Hazelnut (Corylus avellana L.) is an economically and socially important product for Turkey, the country that leads global production of this crop. The preservation of Turkish hazelnut genetic diversity and informed breeding of new cultivars are crucial for maintaining quality and crop yield stability. In this study, genotyping by random amplicon sequencing (GRAS-Di) was used to identify single-nucleotide polymorphisms (SNPs) in a panel of 96 individuals representing the Turkish national hazelnut collection. The resulting 7609 high-quality SNPs were physically mapped to the Tombul cultivar reference genome and used for population structure and diversity analyses. These analyses revealed that cultivars are not less diverse than wild accessions and that 44% of the panel had admixed ancestry. The results also indicated that recently released Turkish cultivars are highly similar to each other, suggesting that diversity analysis is an important tool that should be employed to prevent future genetic bottlenecks in this crop. A minimal marker algorithm was used to select a set of seven SNP markers that were capable of differentiating the panel accessions. These fingerprinting markers should be useful for the propagation of true-to-type elite cultivars that can be used to renew Turkey's aging hazelnut orchards.

Neutral and adaptive genetic diversity in plants: An overview

Genetic diversity is a prerequisite for evolutionary change in all kinds of organisms. It is generally acknowledged that populations lacking genetic variation are unable to evolve in response to new environmental conditions (e.g., climate change) and thus may face an increased risk of extinction. Although the importance of incorporating genetic diversity into the design of conservation measures is now well understood, less attention has been paid to the distinction between neutral (NGV) and adaptive (AGV) genetic variation. In this review, we first focus on the utility of NGV by examining the ways to quantify it, reviewing applications of NGV to infer ecological and evolutionary processes, and by exploring its utility in designing conservation measures for plant populations and species. Against this background, we then summarize the ways to identify and estimate AGV and discuss its potential use in plant conservation. After comparing NGV and AGV and considering their pros and cons in a conservation context, we conclude that there is an urgent need for a better understanding of AGV and its role in climate change adaptation. To date, however, there are only a few AGV studies on non-model plant species aimed at deciphering the genetic and genomic basis of complex trait variation. Therefore, conservation researchers and practitioners should keep utilizing NGV to develop relevant strategies for rare and endangered plant species until more estimates of AGV are available.

Multiplexed ISSR Genotyping by Sequencing (MIG-Seq)

Multiplexed inter-simple sequence repeat (ISSR) genotyping by sequencing (MIG-seq) is a simple, rapid, and inexpensive method for detecting single-nucleotide polymorphisms (SNPs) using next-generation sequencing (NGS). The advantages of MIG-seq include easy application to various species without prior genetic information. In addition, this method opens the door to genome-wide nucleotide sequence analyses of low-quality and trace-level deoxyribonucleic acid (DNA) samples, which have previously been difficult to analyze. Another advantage is that the procedure is simple, time-saving, and inexpensive. Recently, MIG-seq has been applied to wild and cultivated plants and has produced novel results. Using invisible DNA information, questions related to gene flow through pollination and seed dispersal, the genetic structure and diversity of populations, clonality, and the hybridization of wild and cultivated plants are being rapidly answered. In this chapter, I present the results of plant research based on MIG-seq and describe the procedure for this method as a user of MIG-seq.

Genetic population structures of common scavenging species near hydrothermal vents in the Okinawa Trough

Deep-sea mining of hydrothermal deposits off the coast of Japan is currently under consideration, and environmental baseline studies of the area are required to understand possible impacts. The aim of this study is to clarify population structures of dominant benthic megafaunal species near hydrothermal vent fields in the Okinawa Trough, using a population genetics approach. We examined dominant deep-sea scavenging species including eels, several amphipods, and a decapod and performed population genetic analyses based on the mitochondrial cytochrome c oxidase subunit I region. Several sites were sampled within Okinawa Trough to examine intra-population diversity while two other locations 1400-2400 km away were chosen for inter-population comparisons. For synaphobranchid eels Simenchelys parasitica and Synaphobranchus kaupii, our results showed significant intra-population diversity but no inter-population genetic differentiation, suggesting strong genetic connectivity and/or large population sizes. In addition, single nucleotide polymorphism analysis also confirmed strong genetic connectivity for Simenchelys parasitica. Among scavenging amphipods, we detected seven putative species using molecular phylogenetic analysis. We evaluated population structures of the most abundant species of amphipods and a decapod species (Nematocarcinus lanceopes). Our results provide basic information on the genetic population structures of benthic megafaunal species near hydrothermal vent fields, which can be used to select candidate species for future connectivity analysis with high-resolution genetic markers and aid understanding of the potential population impacts of environmental disturbances.

Genomics reveals a genetically isolated population of the Pacific white-sided dolphin (Lagenorhynchus obliquidens) distributed in the Sea of Japan

The processes by which animals become genetically isolated in an open environment such as the ocean have not yet been fully elucidated. Morphologically different populations of Pacific white-sided dolphin Lagenorhynchus obliquidens are observed sympatrically in the Sea of Japan. However, genetic studies that exclusively used limited mitochondrial loci or microsatellite DNA have failed to demonstrate the existence of genetically distinct populations. Here, to reveal the population structure, we analysed genome-wide population genetic data using single nucleotide polymorphisms (SNPs) gathered in 2018-2021 from all domestic captive individuals in aquaria, the majority of which originated from the wild, as well as from some stranded individuals, together covering a wide area of coastal water around Japan (n = 123). Multiplexed intersimple sequence repeat genotyping-by-sequencing analysis was performed to obtain the SNP data. Principal coordinate analysis and the clustering method structure indicated that two genetically-distinct populations exist, with little interpopulation gene flow revealed. In addition, the genotypic segregation was reflected in differences in external morphotype. Furthermore, a population demographic analysis based on the whole-genome sequences of an individual from each population indicated that sea-level changes during the Last Glacial Period probably led to allopatric divergence of this species in a limited area of the Sea of Japan, with that group subsequently sharing a distribution area with the other population. These findings yield insights into the formation of genetically isolated sympatric populations in the ocean.

The contrary conservation situations of two local critically endangered species, Vaccinium emarginatum (Ericaceae) and Elatostema platyphyllum (Urticaceae), growing on the eastern edge of the distribution

As biodiversity loss continues, there is an urgent need to develop efficient conservation measures to protect diversity with limited conservation resources. Conservation targets have generally been selected based on their population size, but more detailed assessments clarifying the phylogenetic genetic status, history, and phylogenetic uniqueness of rare species is crucial to set more appropriate and effective conservation measures. In Japan, the Act on Conservation of Endangered Species of Wild Fauna and Flora designated endangered plants with high conservation priority, but >40% of these species also grow overseas. We conducted comparative analyses based on ddRADseq and MIG-seq to evaluate the population conservation status and value of Vaccinium emarginatum and Elatostema platyphyllum which are growing across national borders at the eastern edge of their species distribution range. The analyses revealed contrasting conservation status between the two species; the Japanese population of V. emarginatum had lower genetic diversity at the individual level and phylogenetically differentiated from Taiwanese populations, while that of E. platyphyllum had higher diversity at the individual level and is a relatively recent migrant with little phylogenetical differentiation from Taiwanese populations. The two species, which share the common feature of being critically rare in Japan, showed contrasting genetic/phylogenetic characteristics. This study provided useful information for appropriate conservation measures based on species’ phylogenetic traits and genetic diversity.

Spatiotemporal process of long-distance seed dispersal in a pantropically distributed sea hibiscus group

Long-distance dispersal (LDD) of seeds plays an essential role in the migration of plants to a new habitat and maintaining gene flow among geographically isolated populations. Pantropical plants with sea-drifted seeds, which have one of the largest distributions in all flowering plants, have achieved their global distribution by LDD. However, the spatiotemporal processes to achieve the wide distribution and the role of LDD in it have not yet been elucidated. In this study, we conducted phylogenomic analyses on the plastome, genome-wide nuclear SNP, and low-copy gene data of Hibiscus tiliaceus and its relatives. The dated phylogeny suggested that global expansion started approximately 4 million years ago (Ma), and species diversification occurred 1 Ma. Plastome phylogeny confirmed the nonmonophyly of the haplotypes in the two widely distributed coastal species, H. tiliaceus and H. pernambucensis. In contrast, genome-wide nuclear SNP phylogenies demonstrated clear genetic segregation among species and/or geographical regions. Ancestral polymorphisms in chloroplast genomes shared among widely distributed species have remained below the range of rapid expansion and speciation of marginal populations. This study demonstrated that the LDD of sea-drifted seeds contributed to the rapid expansion and pantropical distribution of sea hibiscus in the last few million years, and adaptation to local environment or isolation by regional effect after LDD promoted speciation, suppressing gene flow.

Monotropastrum kirishimense (Ericaceae), a new mycoheterotrophic plant from Japan based on multifaceted evidence

Due to their reduced morphology, non-photosynthetic plants have been one of the most challenging groups to delimit to species level. The mycoheterotrophic genus Monotropastrum, with the monotypic species M. humile, has been a particularly taxonomically challenging group, owing to its highly reduced vegetative and root morphology. Using integrative species delimitation, we have focused on Japanese Monotropastrum, with a special focus on an unknown taxon with rosy pink petals and sepals. We investigated its flowering phenology, morphology, molecular identity, and associated fungi. Detailed morphological investigation has indicated that it can be distinguished from M. humile by its rosy pink tepals and sepals that are generally more numerous, elliptic, and constantly appressed to the petals throughout its flowering period, and by its obscure root balls that are unified with the surrounding soil, with root tips that hardly protrude. Based on genome-wide single-nucleotide polymorphisms, molecular data has provided clear genetic differentiation between this unknown taxon and M. humile. Monotropastrum humile and this taxon are associated with different Russula lineages, even when they are sympatric. Based on this multifaceted evidence, we describe this unknown taxon as the new species M. kirishimense. Assortative mating resulting from phenological differences has likely contributed to the persistent sympatry between these two species, with distinct mycorrhizal specificity.

Building consensus around the assessment and interpretation of Symbiodiniaceae diversity

Within microeukaryotes, genetic variation and functional variation sometimes accumulate more quickly than morphological differences. To understand the evolutionary history and ecology of such lineages, it is key to examine diversity at multiple levels of organization. In the dinoflagellate family Symbiodiniaceae, which can form endosymbioses with cnidarians (e.g., corals, octocorals, sea anemones, jellyfish), other marine invertebrates (e.g., sponges, molluscs, flatworms), and protists (e.g., foraminifera), molecular data have been used extensively over the past three decades to describe phenotypes and to make evolutionary and ecological inferences. Despite advances in Symbiodiniaceae genomics, a lack of consensus among researchers with respect to interpreting genetic data has slowed progress in the field and acted as a barrier to reconciling observations. Here, we identify key challenges regarding the assessment and interpretation of Symbiodiniaceae genetic diversity across three levels: species, populations, and communities. We summarize areas of agreement and highlight techniques and approaches that are broadly accepted. In areas where debate remains, we identify unresolved issues and discuss technologies and approaches that can help to fill knowledge gaps related to genetic and phenotypic diversity. We also discuss ways to stimulate progress, in particular by fostering a more inclusive and collaborative research community. We hope that this perspective will inspire and accelerate coral reef science by serving as a resource to those designing experiments, publishing research, and applying for funding related to Symbiodiniaceae and their symbiotic partnerships.

Population Genomic Evidence for the Diversification of Bellamya aeruginosa in Different River Systems in China

Clarifying the genetic structure can facilitate the understanding of a species evolution history. It is crucial for the management of germplasm resources and providing useful guidance for effective selective breeding. Bellamya is an economically and ecologically important freshwater snail for fish, birds and even humans. Population genetic structures of the Bellamya species, however, were unknown in previous studies. Population genomics approaches with tens to hundreds of thousands of single nucleotide polymorphisms (SNPs) make it possible to detect previously unidentified structures. The population genomic study of seven populations of B. aeruginosa across three river systems (Yellow River, Yangtze River and Pearl River) in China was conducted by SLAF-seq. SLAF-seq obtained a total of 4737 polymorphisms SLAF-tags and 25,999 high-consistency genome-wide SNPs. The population genetic structure showed a clear division among populations from the Yellow River basin (YH and WL) and the Pearl River basin (QSH and LB), as well as population YC from the Yangtze River basin using the SNPs data. However, there existed no distinct population structure using the mitochondrial DNA (mtDNA). The anthropogenic translocation from the Yangtze River basin to the Pearl River basin and the passive dispersion from the Yangtze River basin to the Yellow River basin by flooding have weakened the phylogeographic pattern of B. aeruginosa. The divergence of B. aeruginosa in the three river systems suggests that the anthropogenic dispersal for aquaculture and breeding requires serious consideration of the population structure for the preservation of genetic diversity and effective utilization of germplasm resources.